Question: bwa mem mismatch option
gravatar for Varun Gupta
5.7 years ago by
Varun Gupta1.1k
United States
Varun Gupta1.1k wrote:

Hi Everyone,

I am using bwa mem to align my paired end read data. My reads are 150bp in length(paired end). I am using bwa mem as the manual says to use it for read length greater than 70 bp.
My question is I could not find an explicit option to control the sequencing errors i.e number of mismatches. I would like to have around 5 mismatches for 150 bp read I have.

I used this bwa mem command line

/apps1/bwa/bwa-0.7.5a/bwa mem -t 12 -R '@RG\tID:$ID\tPL:illumina\tPU:$PU\tLB:$LB\tSM:$SM' -v 1 -a -M /cork/S.cerevisiae/indexes/bwa/sacCer3.fa file_1.fastq file_2.fastq >

With the above command line I am getting some times more than 11 mismatches in my 150 bp read. Because of this, the read is getting mapped with mismatches at a place where deletion should be present. Any help how to tune in mismatch error.

Also Can some one explain me the parameters

-O INT     gap open penalty [6]

 -L INT     penalty for clipping [5]

I am using bwa-0.7.5a

Hope to hear from you soon.


bwa mem • 4.4k views
ADD COMMENTlink modified 5.7 years ago • written 5.7 years ago by Varun Gupta1.1k
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: 1411 users visited in the last hour