Model-based prediction of sequence alignment quality.
BACKGROUND
Multiple sequence alignment (MSA) is an essential prerequisite for many sequence analysis methods and valuable tool itself for describing relationships between protein sequences. Since the success of the sequence analysis is highly dependent on the reliability of alignments, measures for assessing the quality of alignments are highly requisite.


RESULTS
We present a statistical model-based alignment quality score. Unlike other quality scores, it does not require several parallel alignments for the same set of sequences or additional structural information. Our quality score is based on measuring the conservation level of reference alignments in Homstrad. Reference sequences were realigned with the Mafft, Muscle and Probcons alignment programs, and a sum-of-pairs (SP) score was used to measure the quality of the realignments. Statistical modelling of the SP score as a function of conservation level and other alignment characteristics makes it possible to predict the SP score for any global MSA. The predicted SP scores are highly correlated with the correct SP scores, when tested on the Homstrad and SABmark databases. The results are comparable to that of multiple overlap score (MOS) and better than those of normalized mean distance (NorMD) and normalized iRMSD (NiRMSD) alignment quality criteria. Furthermore, the predicted SP score is able to detect alignments with badly aligned or unrelated sequences.


BACKGROUND
The method is freely available at http://www.mtt.fi/AlignmentQuality/.
Ahola
Virpi
V
Biotechnology and Food Research, MTT Agrifood Research Finland, FI-31600 Jokioinen, Finland. virpi.ahola@mtt.fi
Aittokallio
Tero
T
Vihinen
Mauno
M
Uusipaikka
Esa
E
D000465:Algorithms; D000595:Amino Acid Sequence; D019295:Computational Biology / Q000379:methods*; D030562:Databases, Protein; D015233:Models, Statistical; D016415:Sequence Alignment / Q000592:standards*