Get rid of multimapped and unmapped reads for snp calling improvement
0
0
Entering edit mode
3.0 years ago
Denis ▴ 310

I'm doing snp calling for several projects in which either BWA or Bowtie2 was used to map reads to the reference genome. I use only reads with mapping quality >= 20 in my workflow. Can i be completely sure by doing so that all the multimapped and unmapped reads are not used for snp calling or i have to use other approach to filter out these reads? I.e. may multimapped and unmapped reads have the mapping quality >=20 in the BAM file produced by BWA or Bowtie2 aligners?

genome mapping SNP Illumina • 718 views
ADD COMMENT
2
Entering edit mode

Multimappers have usually a MAPQ of zero and unmapped reads have no MAPQ at all.

ADD REPLY

Login before adding your answer.

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