Bayesian adaptive sequence alignment algorithms.
Bayesian adaptive sequence alignment algorithms.
复制标题
贝叶斯自适应序列比对算法。
DOI:
10.1093/bioinformatics/14.1.25
复制
发表时间:
1998
期刊:
影响因子:
--
通讯作者:
Lawrence,CE
中科院分区:
文献类型:
--
作者:
Zhu,J;Liu,JS;Lawrence,CE
The selection of a scoring matrix and gap penalty parameters continues to be an important problem in sequence alignment. We describe here an algorithm, the 'Bayes block aligner, which bypasses this requirement. Instead of requiring a fixed set of parameter settings, this algorithm returns the Bayesian posterior probability for the number of gaps and for the scoring matrices in any series of interest. Furthermore, instead of returning the single best alignment for the chosen parameter settings, this algorithm returns the posterior distribution of all alignments considering the full range of gapping and scoring matrices selected, weighing each in proportion to its probability based on the data. We compared the Bayes aligner with the popular Smith-Waterman algorithm with parameter settings from the literature which had been optimized for the identification of structural neighbors, and found that the Bayes aligner correctly identified more structural neighbors. In a detailed examination of the alignment of a pair of kinase and a pair of GTPase sequences, we illustrate the algorithm's potential to identify subsequences that are conserved to different degrees. In addition, this example shows that the Bayes aligner returns an alignment-free assessment of the distance between a pair of sequences.