Question: Bwa Bwasw Md Tag
1
gravatar for Random
7.0 years ago by
Random160
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.

Thanks.

bwa • 2.4k views
ADD COMMENTlink written 7.0 years ago by Random160
4
gravatar for Pierre Lindenbaum
7.0 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum118k 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 7.0 years ago by Pierre Lindenbaum118k

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

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