New finite-size correction for local alignment score distributions

BMC Res Notes. 2012 Jun 12:5:286. doi: 10.1186/1756-0500-5-286.

Abstract

Background: Local alignment programs often calculate the probability that a match occurred by chance. The calculation of this probability may require a "finite-size" correction to the lengths of the sequences, as an alignment that starts near the end of either sequence may run out of sequence before achieving a significant score.

Findings: We present an improved finite-size correction that considers the distribution of sequence lengths rather than simply the corresponding means. This approach improves sensitivity and avoids substituting an ad hoc length for short sequences that can underestimate the significance of a match. We use a test set derived from ASTRAL to show improved ROC scores, especially for shorter sequences.

Conclusions: The new finite-size correction improves the calculation of probabilities for a local alignment. It is now used in the BLAST+ package and at the NCBI BLAST web site ( http://blast.ncbi.nlm.nih.gov).

Publication types

  • Research Support, N.I.H., Intramural

MeSH terms

  • Amino Acid Sequence*
  • Databases, Protein
  • Internet
  • Molecular Sequence Data
  • Probability
  • Research Design
  • Sequence Alignment / methods*
  • Sequence Alignment / statistics & numerical data
  • Software*