Hi all,
I am running Kraken2 on some water eDNA samples. I am getting unusually high classification - only around 35% of my reads are classed as 'unclassified'. I know of other people who run similar analysis on water samples and their classification success is something like 0.9% of reads being classified.
Is there any reason why this might be?? Have I somehow set a low stringency for Kraken2 classification, or am I reading the output wrong? I am looking at the k2report file produced and going off that.
1 answer
Kraken 2 can generate information where, for every sequence, it lists what taxonomical hits it has seen and how the lowest common ancestor was determined (LCA)
https://github.com/DerrickWood/kraken2/wiki/Manual#confidence-scoring
that can give you a good sense of how accurate the classifications are
Log in to answer this question.
-
De novo transcriptome assembly (Trinity) with heavy fungal contamination in fruit/flower samples
written by soundslikedaisy •
-
Per base sequence content of Metagenomic samples
written by Savindu •
-
Low percentage of classified reads (~18%) after Kraken2 analysis – is this expected?
written by sumah.1si21bt017 •
-
Kraken2/bracken with 16S sequencing downstream analysis
written by Isidro •
-
extremely low classified reads from sequencing gut microbiome data, applying reads to another refer…
written by jobie1 •
-
Viral taxonomic classification question
written by v.berriosfarias •
-
Many unclassified reads using Kraken2, is it normal?
written by michal.poran •
-
Low Kraken mapping rate for unknown samples
written by LRStar
-
illumina sequencing reads with high quality scores have low percentage mapping to refseq and classi…
written by cecilio11 •
-
In simple words - what is k-mer??
written by alisa.kazarina21 •
Two things to check would be the database against you run your classification, and the
--confidencevalue (which I think is 0 by default).Thanks for your reply - do you mean checking the quality of the database sequences? I built it from only RefSeq sequences of a specific list of organisms, I don't see how this could have caused an issue :(
Sorry I wasn't clear. Since you were comparing your results to others, I suggested that you make sure your DB is comparable to the ones used in other studies. But probably the more important point is to check the confidence value. If I am not mistaken, by default even a single k-mer match is enough to classify a query sequence. So, as suggested by @Istvan , you should probably take a more careful look at the output and determine which taxa your sequences are classified to and with what confidence score.