Question: Qual Versus Bq Tag In Sam/Bam
gravatar for Fixee
9.5 years ago by
Fixee60 wrote:

I've noticed in many of the 1000genome BAMs there are qual strings (Phred), one per base, giving base quality; then there is a BQ tag at the end of the read giving base alignment quality.

What's the difference here? Is the first from the sequencing machine giving confidence on the base call and the second from the assembler saying how good the alignment is?

If true, it seems silly to try and give an alignment quality score for each base.

quality alignment bam • 2.8k views
ADD COMMENTlink modified 9.5 years ago by lh332k • written 9.5 years ago by Fixee60

please paste in a sample line of this BAM/SAM so people can see what you're talking about

ADD REPLYlink written 9.5 years ago by Jeremy Leipzig19k
gravatar for Istvan Albert
9.5 years ago by
Istvan Albert ♦♦ 86k
University Park, USA
Istvan Albert ♦♦ 86k wrote:

I have not needed to make use of this measure in my work but it intrigued me so I researched it a bit:

From the SAM format specification:

BQ = offset to base alignment quality (BAQ), of the same length as the read sequence. At the i-th read base, BAQi = Qi - (BQi - 64) where Qi is the i-th base quality.

The samtools manual states that:

Base Alignment Quality (BAQ) is a new concept deployed in samtools-0.1.9+. It aims to provide an efficient and effective way to rule out false SNPs caused by nearby INDELs.The BAQ is the Phred-scaled probability of a read base being misaligned.

ADD COMMENTlink written 9.5 years ago by Istvan Albert ♦♦ 86k
gravatar for lh3
9.5 years ago by
United States
lh332k wrote:
ADD COMMENTlink written 9.5 years ago by lh332k
Please log in to add an answer.


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