Question: Mismatch Option In Bowtie2
3
gravatar for Varun Gupta
4.2 years ago by
Varun Gupta880
United States
Varun Gupta880 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 • 7.1k views
ADD COMMENTlink modified 20 months ago by poojayakshi10 • written 4.2 years ago by Varun Gupta880
4
gravatar for Ashutosh Pandey
4.2 years ago by
Philadelphia
Ashutosh Pandey11k 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 4.2 years ago by Ashutosh Pandey11k

thanks just what i thought too

Varun

ADD REPLYlink written 4.2 years ago by Varun Gupta880
0
gravatar for kashifalikhan007
20 months ago by
Cologne
kashifalikhan00740 wrote:

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 COMMENTlink modified 20 months ago • written 20 months ago by kashifalikhan00740
Please log in to add an answer.

Help
Access

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