Main Menu (Mobile)- Block

Main Menu - Block

A probabilistic model of local sequence alignment that simplifies statistical significance estimation.

Our Research / Publications
janelia7_blocks-janelia7_biblio_header | block
PLoS Computational Biology. 2008 May 30;4:e1000069. doi: 10.1371/journal.pcbi.1000069
View Article on
A probabilistic model of local sequence alignment that simplifies statistical significance estimation. Eddy/Rivas Lab
Eddy Sean
janelia7_blocks-janelia7_biblio_abstract | block

Abstract

Sequence database searches require accurate estimation of the statistical significance of scores. Optimal local sequence alignment scores follow Gumbel distributions, but determining an important parameter of the distribution (lambda) requires time-consuming computational simulation. Moreover, optimal alignment scores are less powerful than probabilistic scores that integrate over alignment uncertainty ("Forward" scores), but the expected distribution of Forward scores remains unknown. Here, I conjecture that both expected score distributions have simple, predictable forms when full probabilistic modeling methods are used. For a probabilistic model of local sequence alignment, optimal alignment bit scores ("Viterbi" scores) are Gumbel-distributed with constant lambda = log 2, and the high scoring tail of Forward scores is exponential with the same constant lambda. Simulation studies support these conjectures over a wide range of profile/sequence comparisons, using 9,318 profile-hidden Markov models from the Pfam database. This enables efficient and accurate determination of expectation values (E-values) for both Viterbi and Forward scores for probabilistic local alignments.

PMID: 18516236 [PubMed - indexed for MEDLINE]

node:body | entity_field
janelia7_blocks-janelia7_biblio_authors | block

Janelia Authors

janelia7_blocks-janelia7_biblio_tools | block