Question: Mapping of metagenomic reads on genomes - length criteria ?
gravatar for lagartija
5 weeks ago by
lagartija60 wrote:


I am mapping reads to genomes to estimate their abundance in several metagenomic samples. My problem is that when I filter for reads that are over 33 nt long I get a complete different result than when I keep all the reads. Which criterias should I apply to do the mapping ?

Thank you very much

mapping • 98 views
ADD COMMENTlink written 5 weeks ago by lagartija60

Doesn't that make sense? If you use different input reads you will get different results. It would be useful if you could explain the reason for deleting reads longer than 33 nt?

PS: Please refer to: Brief Reminder On How To Ask A Good Question and How To Ask Good Questions On Technical And Scientific Forums

ADD REPLYlink modified 5 weeks ago • written 5 weeks ago by Sej Modha4.2k

The reason that I thought about removing short reads is to get less false positives. I looked at the first alignments and saw that most of the matches were short reads, often with repeats. I thought that if only short reads are matching, it might be only false positives.

ADD REPLYlink written 5 weeks ago by lagartija60
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: 1141 users visited in the last hour