bowtie2 mapping qualities of 0
1
0
Entering edit mode
3.1 years ago
Egelbets ▴ 30

I have used bowtie2 to map paired-end reads together with unpaired reads to a reference genome (default parameters). I noticed when looking at my alignment in IGV, that at some locations there are a lot of reads with a mapping quality of 0. I know that in BWA a mapping quality of 0 indicates that the read had multiple equally good mapping locations and that BWA chose one alignment at random. My question is, does this also apply to bowtie2 alignments? The bowtie2 documentation does state that that the tool chooses a random alignment in case of multiple equally good alignments, but it doesn't say anything on the mapping quality it gives those reads.

bowtie2 • 1.7k views
ADD COMMENT
2
Entering edit mode
3.1 years ago

Yes, bowtie2 assigns a mapping quality of 0 (EDIT: 0 or 1 as per kalavattam's comment) to multi-mapping reads.

ADD COMMENT
0
Entering edit mode

Thank you!

ADD REPLY
0
Entering edit mode

This is not correct. In addition to MAPQ = 0, Bowtie2 can assign MAPQ = 1 to multimappers. Details here.

ADD REPLY

Login before adding your answer.

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