jenniferlu717 / KrakenTools

KrakenTools provides individual scripts to analyze Kraken/Kraken2/Bracken/KrakenUniq output files
GNU General Public License v3.0
297 stars 85 forks source link

kreport2krona.py not count the reads of bracken.Kreport #44

Open gracielad opened 2 years ago

gracielad commented 2 years ago

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 kBacteria pFirmicutes cBacilli oLactobacillales 0 kBacteria pFirmicutes cBacilli oLactobacillales fLactobacillaceae 0 kBacteria pFirmicutes cBacilli oLactobacillales fLactobacillaceae gLigilactobacillus 0 kBacteria pFirmicutes cBacilli oLactobacillales fLactobacillaceae gLimosilactobacillus 0 kBacteria pFirmicutes cBacilli oLactobacillales fLactobacillaceae gLactobacillus 0 kBacteria pFirmicutes cBacilli oLactobacillales fLactobacillaceae g__Leuconostoc

Somebodyatthdoor commented 1 week ago

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.