Question: Mapping of metagenomic reads on genomes - length criteria ?
gravatar for lagartija
13 months ago by
lagartija80 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 • 284 views
ADD COMMENTlink written 13 months ago by lagartija80

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 13 months ago • written 13 months ago by Sej Modha4.7k

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 13 months ago by lagartija80
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: 1706 users visited in the last hour