Question: How samtools calculate error rate
0
gravatar for jianzheng934963534
2.1 years ago by
jianzheng93496353410 wrote:

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.

sequencing samtools • 1.6k views
ADD COMMENTlink modified 2.1 years ago by Devon Ryan91k • written 2.1 years ago by jianzheng93496353410
0
gravatar for Devon Ryan
2.1 years ago by
Devon Ryan91k
Freiburg, Germany
Devon Ryan91k wrote:

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.

ADD COMMENTlink written 2.1 years ago by Devon Ryan91k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1664 users visited in the last hour