Many unclassified reads using Kraken2, is it normal?
0
0
Entering edit mode
21 months ago

Hi all, I ran kraken2 analysis on my dataset (8 samples from 2 different human gut micrbiome projects) and got above 50% no hits. I tried to ran those unclassified reads with Mataphlan2, and didn't get any classification either. Can you suggest a reason for that many unclassified reads, can I avoid that? Thanks in advance Michal

kraken2 kraken Mataphlan2 Mataphlan • 2.0k views
ADD COMMENT
0
Entering edit mode

How long are the reads?

ADD REPLY
0
Entering edit mode

We tried 250bp and 150bp reads

ADD REPLY
0
Entering edit mode

I am not familiar with kraken2 or even metaphlan but there can be many reasons. It also depends on your per-processing steps. Think the common reasons are

  • The DNA is just not sequenced yet
  • Your database is incomplete (most common)
  • The input data just contains a lot of non-coding DNA or no marker genes

Think if you want an answer from an expert you need to tell if you did any filtering before and which reference you are using.

ADD REPLY
0
Entering edit mode

First thanks for your help. I think its technical since the sample were prepared in different countries, so it is less likely they share the same contamination, I also tried to blast some sequence and got low alignment score, so I guess its not the db but i will try to run some more sequences to be extra sure it is not the problem. I wondered if it is common to have such low detection rate... Thanks again, Michal

ADD REPLY
0
Entering edit mode

You can blast and have a perfect hit but not all parts of the DNA can be used to determine a species.

ADD REPLY
0
Entering edit mode

I didn't filter the data and used the Kraken2 db.

ADD REPLY
0
Entering edit mode

Think it will always help do to at least some quality trimming

ADD REPLY
0
Entering edit mode

michal.poran : Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized. SUBMIT ANSWER is for new answers to original question.

ADD REPLY

Login before adding your answer.

Traffic: 1082 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6