Parametric analysis for ungapped Markov models of evolution

  • Authors:
  • David Fernández-Baca;Balaji Venkatachalam

  • Affiliations:
  • Iowa State University, Ames, IA;University of California, Davis, CA

  • Venue:
  • ACM Transactions on Algorithms (TALG)
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

Efficient sensitivity analysis algorithms are presented for two problems arising in the study of Markov models of sequence evolution: ancestral reconstruction in evolutionary trees and local ungapped alignment under log-odds scoring. The algorithms generate complete descriptions of the optimum solutions for all possible values of the evolutionary distance. The running time for the parametric ancestral reconstruction problem under the Kimura 2-parameter model is O(kn + kn2/3 log k), where n is the number of sequences and k is their length, assuming all edges have the same length. For the parametric gapless alignment problem under the Jukes-Cantor model, the running time is O(mn + mn2/3 log m), where m and n are the sequence lengths and n ≤ m.