Mapping of metagenomic reads on genomes - length criteria ?
0
0
Entering edit mode
5.5 years ago
lagartija ▴ 160

Hi,

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 • 908 views
ADD COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY

Login before adding your answer.

Traffic: 1532 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