Hello,
I'm looking for a script to extract from sam file the number of M, (match base in the CIGAR string) and the number mismatch (NM:i field).
The objective is to calculate a mismatch frequency in order to filter reads.
Thanks
Cheers
Cédric
I'm looking for a script to extract from sam file the number of M, (match base in the CIGAR string) and the number mismatch (NM:i field).
The objective is to calculate a mismatch frequency in order to filter reads.
Thanks
Cheers
Cédric
Comment