kreport2krona.py not count the reads of bracken.Kreport
gracielad opened this issue · 2 comments
To obtain the krona visualization of Bracken reports, I run the kreport2krona.py using the file bracken_kreport (percentage of reads, Total number of reads, etc). However, I noted that the created output file does not present the majority of reads assigned in the bracken_Kreport. For example, in the bracken_kreport,
a vast amount of reads assigned of genus Ligilactobacillus and the created file the number is 0. Anyone, could you help me to solve this?
100.00 7095261 0 R 1 root
100.00 7095139 0 R1 131567 cellular organisms
99.69 7073316 0 D 2 Bacteria
88.16 6255049 0 D1 1783272 Terrabacteria group
86.95 6169667 0 P 1239 Firmicutes
86.60 6144312 0 C 91061 Bacilli
86.43 6132776 0 O 186826 Lactobacillales
86.04 6104444 0 F 33958 Lactobacillaceae
83.89 5952280 0 G 2767887 Ligilactobacillus
1.53 108246 0 G 2742598 Limosilactobacillus
0 k__Bacteria p__Firmicutes c__Bacilli o__Lactobacillales
0 k__Bacteria p__Firmicutes c__Bacilli o__Lactobacillales f__Lactobacillaceae
0 k__Bacteria p__Firmicutes c__Bacilli o__Lactobacillales f__Lactobacillaceae g__Ligilactobacillus
0 k__Bacteria p__Firmicutes c__Bacilli o__Lactobacillales f__Lactobacillaceae g__Limosilactobacillus
0 k__Bacteria p__Firmicutes c__Bacilli o__Lactobacillales f__Lactobacillaceae g__Lactobacillus
0 k__Bacteria p__Firmicutes c__Bacilli o__Lactobacillales f__Lactobacillaceae g__Leuconostoc
I'm aware this is a very old post, but I am getting a very similar problem. The krona plots being produced don't seem reflective of the bracken input.
I am facing the same problem with the krona output generated from bracken.k2report using kreport2krona.py
script. There are missing reads that have been previously assigned to relevant taxa in krona text file output.