The BLAST computer program presently uses computer simulations to give accurate estimates of statistical significance for sequence matches. Recent advances in the statistical theory of gapped sequence alignments indicates a need for further rigorous mathematical results in the underlying statistical theory. This project aims to provide those results. In particular, I have developed as mathematical proofs the limiting evolution of the composition of high scoring segments as a Brownian motion, a result that can be used to improve the accuracy of approximations to the statistics of gapped alignments. In addition, the edge effects present because real sequences have limited lengths appear as a correction term in an asymptotic expansion of the probability of sequence matching. This edge effect is likely to be more important in the statistics of matching with gaps than it was in the statistics of matching without gaps, because gapped matches tend to be longer, exhausting the sequences being matched more easily. As a byproduct, this project has produced some new analytic results that are likely to be useful in detecting non-random positioning of genomic markers.

Agency
National Institute of Health (NIH)
Institute
National Library of Medicine (NLM)
Type
Intramural Research (Z01)
Project #
1Z01LM000088-02
Application #
6432764
Study Section
(CBB)
Project Start
Project End
Budget Start
Budget End
Support Year
2
Fiscal Year
2000
Total Cost
Indirect Cost
Name
National Library of Medicine
Department
Type
DUNS #
City
State
Country
United States
Zip Code
Spouge, John L (2007) Inequalities on the Overshoot beyond a Boundary for Independent Summands with Differing Distributions. Stat Probab Lett 77:1486-1489
Sheetlin, Sergey; Park, Yonil; Spouge, John L (2005) The Gumbel pre-factor k for gapped local alignment can be estimated from simulations of global alignment. Nucleic Acids Res 33:4987-94
Frith, Martin C; Spouge, John L; Hansen, Ulla et al. (2002) Statistical significance of clusters of motifs represented by position specific scoring matrices in nucleotide sequences. Nucleic Acids Res 30:3214-24
Park, Yonil; Spouge, John L (2002) The correlation error and finite-size correction in an ungapped sequence alignment. Bioinformatics 18:1236-42
Makalowska, I; Ferlanti, E S; Baxevanis, A D et al. (1999) Histone Sequence Database: sequences, structures, post-translational modifications and genetic loci. Nucleic Acids Res 27:323-4
Wolfsberg, T G; Makalowska, I; Makalowski, W (1999) Genomes and evolution. Web alert. Curr Opin Genet Dev 9:619