Question: Kraken - incomplete taxonomy levels/ranks in assignment
gravatar for medelin.ocejo
16 days ago by
medelin.ocejo0 wrote:

Hi all! In the output of Kraken, I got many assignments with uncomplete taxonomic levels like these:


In the first case, Phylum jumps to Order and jumps to Genus. In the second, Class jumps directly to Species. So, when I try to import the abundance table to phyloseq, the taxas are not arranged correctly to its taxonomic levels. How could I fix this? I have been days browsing and could not find anything. :(

Thank you!

ADD COMMENTlink written 16 days ago by medelin.ocejo0

Have you used the scripts provided within the Kraken package to extract this info from Kraken output?

ADD REPLYlink written 16 days ago by Sej Modha2.7k

Hi @Sej Modha This is the script ran:

kraken --threads 20 --db /data/db_dir/krakken_db Sample1_R1.bam Sample1_R2.bam --paired > kraken_out/Sample1.kraken
kraken-mpa-report --db /data/db_dir/krakken_db Sample1.kraken > kraken_report/Sample1.krakenReport
ADD REPLYlink written 9 days ago by medelin.ocejo0

I believe that Kraken works with .fasta or .fastq files? Are you sure you are getting correct results with .bam files? Kraken relies on the NCBI taxonomy data so you might want to ensure that you are using the most recent version of the Kraken and corresponding taxonomy databases.

ADD REPLYlink written 9 days ago by Sej Modha2.7k

Maybe Braken can help you.

ADD REPLYlink written 16 days ago by h.mon12k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1507 users visited in the last hour