The correlation error and finite-size correction in an ungapped sequence alignment.
The correlation error and finite-size correction in an ungapped sequence alignment.
复制标题
无间隙序列比对中的相关误差和有限尺寸校正。
DOI：
10.1093/bioinformatics/18.9.1236
复制
发表时间：
2002
期刊：
影响因子：
--
通讯作者：
Spouge,JohnL
中科院分区：
文献类型：
--
作者：
Park,Yonil;Spouge,JohnL
Motivation: The BLAST program for comparing two sequences assumes independent sequences in its random model. The resulting random alignment matrices have correlations across their diagonals. Analytic formulas for the BLASTp-value essentially neglect these correlations and are equivalent to a random model with independent diagonals. Progress on the independent diagonals model has been surprisingly rapid, but the practical magnitude of the correlations it neglects remains unknown. In addition, BLAST uses a finite-size correction that is particularly important when either of the sequences being compared is short. Several formulas for the finite-size correction have now been given, but the corresponding errors in the BLASTp-values have not been quantified. As the lengths of compared sequences tend to infinity, it is also theoretically unknown whether the neglected correlations vanish faster than the finite-size correction.Results: Because we required certain analytic formulas, our study restricted its computer experiments to ungapped sequence alignment. We expect some of our conclusions to extend qualitatively to gapped sequence alignment, however. With this caveat, the finite-size correction appeared to vanish faster than the neglected correlations. Although the finite-size correction underestimated the BLASTp-value, it improved the approximation substantially for all but very short sequences. In practice, the Altschul–Gish finite-size correction was superior to Spouge's. The independent diagonals model was always within a factor of 2 of the true BLASTp-value, although fittingp-value parameters from it probably is unwise.Contact: spouge@ncbi.nlm.nih.gov
