That's great, what exactly do I have to check in the NM tag?
• 0 views
•
link
written
by
omariqbal2017 •
How to check if a read in a SAM file contains substitutions
How I do check whether a read contains substitutions or not. I'm not sure if I should check the flag, the cigar or the sequence.
• 2,864 views
•
link
updated
by
Ram
•
written
by
omariqbal2017 •
•
written
by
omariqbal2017 •
2 answers
Typically this is indicated in the NM auxiliary tag. Since you have C++ in the tags, I'll note that you can retrieve this in htslib via uint8_t *v = bam_aux_get(b, "NM") and then int64_t NM = bam_aux2i(v);.
• 0 views
•
link
written
by
Devon Ryan
I gave you two C commands that use htslib in my original reply.
• 0 views
•
link
written
by
Devon Ryan
True, one would also need to parse the CIGAR string and subtract appropriately.
• 0 views
•
link
written
by
Devon Ryan
If you are looking exclusively for base substitution for alignments that use M CIGAR operators (as opposed to X and =), then you either have to a) check each base against the reference and calculate mismatches yourself, or b) use the NM tag and adjust for any indels in the alignment.
To add to this answer, A while back a wrote a program, ExpandCigar, that expands the M operator with 'X' or '='. After that, it's easy to extract reads where the CIGAR string contains X.
Log in to answer this question.
More posts like this
-
gene strand determined by flag in read1 and read2 should be intepreted differently?
written by tomas4482 •
-
Determine fragment alignment overlaps by parsing CIGAR string in the MC tag
written by abcoxyzide •
-
Handling of Reverse Strand in GAM file
written by ved_vyas •
-
How to check if a read in a Bam file is paired or not?
written by omariqbal2017 •
-
How to convert from a cigar string to extended cigar string?
written by omariqbal2017 •
-
What's the difference between half mapped reads and single paired reads?
written by omariqbal2017 •
-
How to deal with undetected probes in microarray data?
written by Kristin Muench
-
How to check whether a genome is circular
written by utkarsh.sood •
-
How to check whether the mate read is fully mapped or not?
written by chongchu.cs •
-
Interpreting "*" In Extended Cigar Format
written by anand.sarwate •
