Mapping of metagenomic reads on genomes - length criteria ?
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
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Filtering HMMSCAN Domtblout Output when using Metagenomic Sequence as Input
written by Jane123 •
-
filtering SAM/BAM to remove hits spanning short combined alignment lengths and low counts
written by Bryan •
-
Calculating read abundances across MAGs (question)
written by v.berriosfarias •
-
How to normalize gene abundance from metagenomic data?
written by v.berriosfarias •
-
Custom Database in Kraken2: how are plasmids in same fasta file as bacterial genome handled?
written by Raphaela •
-
small representative database
written by lagartija
-
Incomplete filtration of reads by length in .fastq and in .bam files
written by kerenyifarkas •
-
Resfams Database usage
written by ginna •
-
Is estimated coverage of Nonpareil reliable?
written by liuyifan2014 •
-
How to calculate organism relative abundance from Illumina metagenomic data?
written by nkuyfq
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
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.