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

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

d__Bacteria|p__Cyanobacteria|o__Pleurocapsales|g__Stanieria 
d__Bacteria|p__Proteobacteria|c__Gammaproteobacteria|s__Calyptogena_okutanii_thioautotrophic_gill_symbiont

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 3 months ago by medelin.ocejo0
1

Maybe Braken can help you.

ADD REPLYlink written 3 months ago by h.mon15k

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

ADD REPLYlink written 3 months ago by Sej Modha2.9k

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 3 months ago by medelin.ocejo0
1

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 3 months ago by Sej Modha2.9k

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

ADD REPLYlink written 10 weeks ago by medelin.ocejo0
Please log in to add an answer.

Help
Access

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