This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to get an exact match from the bam alignment file

How to get an exact match from the bam comparison file, it is well known that M in the cigar character is completely consistent or there is mismatch, but how do I get a completely consistent sequence?

bam bwa samtools
• 2,269 views • link
updated by yasokannan93 • yasokannan93 • written by xiaoguang xiaoguang

1 answer

You can use a combination of read length and the cigar character M from the alignment file. For example, if the read length is 100, you can try extracting reads with cigar score of 100M. This would mean the read is completely consistent or it is an exact match.

• 0 views • link
written by yasokannan93 • yasokannan93

That is not true, a M can also be a mismatch as xiaoguang correctly says. M only says that these was neither insertion nor deletion but a single base found at that position. A possible strategy could be no check the MD tag which should contain the position of mismatches. If it is 0 there should not be one, check the SAM documentation and other threads to confirm this https://samtools.github.io/hts-specs/SAMtags.pdf

• 0 views • link
written by ATpoint ATpoint

Thanks for the clarification. So for the example, we should expect it as MD:Z:100 ?

https://samtools.github.io/hts-specs/SAMtags.pdf says

The MD field aims to achieve SNP/indel calling without looking at the reference. For example, a string ‘10A5^AC6’ means from the leftmost reference base in the alignment, there are 10 matches followed by an A on the reference which is different from the aligned read base; the next 5 reference bases are matches followed by a 2bp deletion from the reference; the deleted sequence is AC; the last 6 bases are matches. The MD field ought to match the CIGAR string.

• 0 views • link
written by yasokannan93 • yasokannan93

Log in to answer this question.

More posts like this