Robust E-values for gapped local alignments.
We examine a Poisson heuristic for judging the significance of local sequence alignments with gaps. Model parameters are estimated directly from the sequences to be aligned, so that laborious prior simulation studies or database comparisons for the estimation of parameters describing the connection between score and E-value are unnecessary. Simulation studies give evidence that this method gives reasonable results even when the usual assumptions like the independence of sequence positions are violated.
Metzler
Dirk
D
Institut für Informatik, Johann Wolfgang Goethe-Universität, Frankfurt am Main, Germany. metzler@cs.uni-frankfurt.de
D000595:Amino Acid Sequence; D001483:Base Sequence; D019295:Computational Biology; D003627:Data Interpretation, Statistical; D016012:Poisson Distribution; D016415:Sequence Alignment / Q000706:statistics & numerical data*