Question: Kraken - incomplete taxonomy levels/ranks in assignment
gravatar for medelin.ocejo
21 months 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 modified 14 months ago by Biostar ♦♦ 20 • written 21 months ago by medelin.ocejo0

Maybe Braken can help you.

ADD REPLYlink written 21 months ago by h.mon28k

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

ADD REPLYlink written 21 months ago by Sej Modha4.5k

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 21 months 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 21 months ago by Sej Modha4.5k

It seems that the problem is in the NCBI database. Thanks for your reply!

ADD REPLYlink written 20 months ago by medelin.ocejo0
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: 1405 users visited in the last hour