Question: Mismatch Option In Bowtie2
gravatar for Varun Gupta
6.9 years ago by
Varun Gupta1.1k
United States
Varun Gupta1.1k wrote:

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 • 11k views
ADD COMMENTlink modified 4.3 years ago by poojayakshi10 • written 6.9 years ago by Varun Gupta1.1k
gravatar for Ashutosh Pandey
6.9 years ago by
Ashutosh Pandey12k wrote:

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 COMMENTlink written 6.9 years ago by Ashutosh Pandey12k

thanks just what I thought too


ADD REPLYlink modified 4 months ago by RamRS27k • written 6.9 years ago by Varun Gupta1.1k
gravatar for kashiff007
4.3 years ago by
kashiff007130 wrote:


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 COMMENTlink modified 4 months ago by RamRS27k • written 4.3 years ago by kashiff007130
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1850 users visited in the last hour