Sequence Alignment/Map (SAM) files

I was trying to get all the mismatch from a sam file and hey, the start position of reads written in a sam file is calculated after soft clipped nucleotides!

And ‘N’ in reference sequence seems to mess up mismatch tag

edited: samtools mpileup output a nice file for mismatch parsing.




Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License. If you liked this post, you can share it with your followers or follow me on Twitter!