Question: Qual Versus Bq Tag In Sam/Bam
3
gravatar for Fixee
8.5 years ago by
Fixee60
Colorado
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.5k views
ADD COMMENTlink modified 8.5 years ago by lh331k • written 8.5 years ago by Fixee60
1

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

ADD REPLYlink written 8.5 years ago by Jeremy Leipzig18k
3
gravatar for Istvan Albert
8.5 years ago by
Istvan Albert ♦♦ 82k
University Park, USA
Istvan Albert ♦♦ 82k 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 8.5 years ago by Istvan Albert ♦♦ 82k
1
gravatar for lh3
8.5 years ago by
lh331k
United States
lh331k wrote:
ADD COMMENTlink written 8.5 years ago by lh331k
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: 834 users visited in the last hour