Question: Bwa Bwasw Md Tag
gravatar for Random
8.0 years ago by
Random160 wrote:

I mapped some reads to a reference genome using bwa's bwasw option and it seems that the MD tag is never produced. Is there a way to make bwa create this tag for each mapped read?

Alternatively is there way to enforce bwa to output a CIGAR string with "X" and "=" characters (which stand respectively for sequence mismatch and match), instead of the ambiguous M (which stands for both match and mismatch) ?

I'm trying to find out how much of each mapped read mismatched the reference.


bwa • 2.7k views
ADD COMMENTlink written 8.0 years ago by Random160
gravatar for Pierre Lindenbaum
8.0 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum127k wrote:

did you try to use samtools calmd ?

samtools calmd [-EeubSr] [-C capQcoef] <aln.bam> <ref.fasta>

generate the MD tag. If the MD tag is already present, this command will give a warning if the MD tag generated is different from the existing tag. Output SAM by default.

ADD COMMENTlink written 8.0 years ago by Pierre Lindenbaum127k

I hadn't. This seems to solve my problem. Thanks

ADD REPLYlink written 8.0 years ago by Random160
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: 2080 users visited in the last hour