KrakenTools
KrakenTools copied to clipboard
kreport2krona.py not count the reads of bracken.Kreport
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.