This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How samtools calculate error rate

I have aligned the fastq file to the reference genome and got a sam file. I know use "samtools stats" can get "error rate" of the sam file.

Is there any document shows how the error rate is calculated from the sam file? (since I want to seperately calculate the error rate of every read in the sam file instead of the average of the whole file)

I am really confused how they deal with softclips and hardclips when calculating the error rate.

Any answer is appreciated.

samtools sequencing
• 6,011 views • link
updated by Devon Ryan Devon Ryan • written by jianzheng934963534 • jianzheng934963534

1 answer

samtools adds up the number of mismatches (from the NM auxiliary tag) and divides that by the number of aligned bases. Soft and hard-clipped bases wouldn't be included, since they aren't aligned.

• 0 views • link
written by Devon Ryan Devon Ryan

Log in to answer this question.

More posts like this