MAPping Quality v.s. Alignment Score (AS) in SAM format

xiangwulu · 03-02-2016, 02:03 AM

MAPQ: MAPping Quality. It equals −10 log10 Pr{mapping position is wrong}, rounded to the nearest
integer. A value 255 indicates that the mapping quality is not available.

AS i Alignment score generated by aligner

What are the differences between these two? I am confused about this. If I want to filter the alignments and keep the 'good' alignment, which value should I look into?

Thanks.

dpryan · 03-02-2016, 02:27 AM

Alignment scores are internal to each aligner. MAPQ scores produced by the aligners typically involves the alignment score and other information. Filter by MAPQ.

xiangwulu · 03-02-2016, 03:58 AM

Thank you.

dariober · 03-02-2016, 04:05 AM

Just to expand a bit on Devon's answer... Alignment score is a metric that tells you how similar the read is to the reference. AS increases with the number of matches and decreases with the number of mismatches and gaps (rewards and penalties for matches and mismatches depend on the scoring matrix you use). MAPQ is a metric that tells you how confident you can be that the read comes from the reported position.

You can have high AS and low MAPQ if the read aligns perfectly at multiple positions, and you can have low AS and high MAPQ if the read aligns with mismatches but still the reported position is still much more probable than any other.

As Devon says, you probably want to filter for MAPQ, but "good" alignment may refer to AS if what you care is similarity between read and reference.

(Correct me if I'm wrong)

xiangwulu · 03-02-2016, 05:35 AM

Quote:

Originally Posted by dariober

Just to expand a bit on Devon's answer... Alignment score is a metric that tells you how similar the read is to the reference. AS increases with the number of matches and decreases with the number of mismatches and gaps (rewards and penalties for matches and mismatches depend on the scoring matrix you use). MAPQ is a metric that tells you how confident you can be that the read comes from the reported position.

You can have high AS and low MAPQ if the read aligns perfectly at multiple positions, and you can have low AS and high MAPQ if the read aligns with mismatches but still the reported position is still much more probable than any other.

As Devon says, you probably want to filter for MAPQ, but "good" alignment may refer to AS if what you care is similarity between read and reference.

(Correct me if I'm wrong)

Thank you for the explanation.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Detecting FASTQ Quality Score Format	bhan	Bioinformatics	2	01-02-2012 04:40 PM
Quality Score Format	vedbar	Bioinformatics	5	07-19-2011 04:41 AM
mapping quality score in tophat sam file	jjpurwar	Bioinformatics	5	06-22-2011 12:55 PM
FastX Quality Score format	kbushley	Illumina/Solexa	1	09-05-2010 02:56 AM
Different MAPPING QUALITY/PER-BASE QUALITY SCORE	m_elena_bioinfo	Bioinformatics	2	09-02-2010 10:00 AM

03-02-2016, 02:03 AM	#1
xiangwulu Member Location: ireland Join Date: Apr 2014 Posts: 18	MAPping Quality v.s. Alignment Score (AS) in SAM format MAPQ: MAPping Quality. It equals −10 log10 Pr{mapping position is wrong}, rounded to the nearest integer. A value 255 indicates that the mapping quality is not available. AS i Alignment score generated by aligner What are the differences between these two? I am confused about this. If I want to filter the alignments and keep the 'good' alignment, which value should I look into? Thanks.

03-02-2016, 02:27 AM	#2
dpryan Devon Ryan Location: Freiburg, Germany Join Date: Jul 2011 Posts: 3,480	Alignment scores are internal to each aligner. MAPQ scores produced by the aligners typically involves the alignment score and other information. Filter by MAPQ.

03-02-2016, 03:58 AM	#3
xiangwulu Member Location: ireland Join Date: Apr 2014 Posts: 18	Thank you.

03-02-2016, 04:05 AM	#4
dariober Senior Member Location: Cambridge, UK Join Date: May 2010 Posts: 311	Just to expand a bit on Devon's answer... Alignment score is a metric that tells you how similar the read is to the reference. AS increases with the number of matches and decreases with the number of mismatches and gaps (rewards and penalties for matches and mismatches depend on the scoring matrix you use). MAPQ is a metric that tells you how confident you can be that the read comes from the reported position. You can have high AS and low MAPQ if the read aligns perfectly at multiple positions, and you can have low AS and high MAPQ if the read aligns with mismatches but still the reported position is still much more probable than any other. As Devon says, you probably want to filter for MAPQ, but "good" alignment may refer to AS if what you care is similarity between read and reference. (Correct me if I'm wrong)