Mismatch Option In Bowtie2
2
4
Entering edit mode
8.9 years ago
Varun Gupta ★ 1.2k

Hi Everyone

I am using bowtie2. I want to set my max mismatches to 2 for my read sequence.

Can this be done in bowtie2. I know bowtie1 has option of -v but i cannot find it in bowtie2.

Any help is appreciated

Regards Varun

bowtie2 • 14k views
ADD COMMENT
7
Entering edit mode
8.9 years ago

You can specify number of mismatch in a seed alignment using -N. Your alignment can still have more than -N mismatches as there can be mismatches in non-seed part of the read. If you want you can specify the seed length to some number that is comparable to the read length. But this way you will miss loose lot of alignments as it will be too stringent. If I would be in your place, I will use the default alignment parameter and then filter out all the reads that have two or less mismatches from the output sam file. You can use XM tag from the SAM file to do the same.

XM:i:<n> The number of mismatches in the alignment. Only present if SAM record is for an aligned read

ADD COMMENT
0
Entering edit mode

thanks just what I thought too

Varun

ADD REPLY
3
Entering edit mode
6.3 years ago
kashiff007 ★ 1.8k

Ashutosh

Filtering through XM tag for single end read is OK, but in case of pair-end I have a confusion because both pairs might have different number of mismatches and in this way one can miss a pair in all such cases.

ADD COMMENT

Login before adding your answer.

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