Dear frineds, i wanna ask how can i know the mismatch from bowtie2 output
dear friends:
i wanna know how to know the mismatch for the bowtie2 output??? when i run the bowtie, the sam output would do like this:(and tell you the base conversion of each reads, also the site)
J00142:71:HGL3TBBXX:8:1101:12540:1666 2:N:0:1 + chr17 8773052 GCACACGGCTGAACGCCAG <FJJJ7FAJF7<<FFAA-7 0 12:C>A
J00142:71:HGL3TBBXX:8:1101:2717:1666 2:N:0:1 - chr2 33141441 GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGAA JJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJFFAAAJAA<77-JJJJJJJJJJJJFJJJJFF7JJJJAA7A77<<-< 51:G>A,32:A>G,72:A>G
but bowtie2 generate like this: (how can i know the base conversion from it or can you help me explain the meaning of the tag MD:Z:21A41A6T12T0T0T0)
J00142:71:HGL3TBBXX:8:1101:1164:1683 16 chrY 10037829 3 49M1I37M * 0 0 TGTGAATTGCAGGACACATTGCTCATCGACACTTCGAACGCACTTGCGGCCCCGGGTTCCTCCCGGGGCTACGCCTGTCTGAGCGCC JAJ7A7AA7-7AFA7J-AF-F-JFFJF<-F-JJFJFFJFAJJF7FJ<FA7-JA<AFFA<FFJAJJJJFFJJ<FJJJJJFFJAF<AAA AS:i:-36 XN:i:0 XM:i:6 XO:i:1 XG:i:1 NM:i:7 MD:Z:21A41A6T12T0T0T0 YT:Z:UU
thank you so much!!!
• 2,459 views
•
link
updated
by
michael.ante
•
written
by
biologo •
•
written
by
biologo •
1 answer
Hi biologo,
you can find the specification of the MD tag here: SAMtags.pdf .
MD:Z:21A41A means, you have 21 matches followed by a mismatch with an A in the reference. The next 41 bases are matches again, followed by a mismatch with an A in the reference.
I hope this will help for the start.
Cheers,
Michael
• 0 views
•
link
written
by
michael.ante
Log in to answer this question.
If I understad the question right
from this value
more could be found here point 1.4 section 6
regarding the explanation
more info
I am very grateful for your kind help!!! thank you
I find the MD tags to be a pain to deal with. Fortunately, the current Sam spec allows you to generate sam files with cigar strings using X to indicate mismatch and = to indicate match. If you map your reads with BBMap, it will do this by default, then you don't have to look at the MD tag. But for your existing mapped reads, you can also convert them with BBMap's Reformat utility, like this:
That will change reads with a cigar string like "100M" to something like "20=1X79=" which clearly indicates there is a mismatch at position 21.