Biases in amino acid replacement matrices and alignment scores due to rate heterogeneity.
Empirically derived amino acid replacement matrices are widely used in sequence comparison and database searches. We consider an extension of the usual Markov process model of protein evolution that admits site to site rate heterogeneity and demonstrates that rate heterogeneity can introduce a bias in estimated replacement probabilities and the corresponding alignment scores derived from these matrices. We suggest an approach to obtain unbiased estimates of replacement probabilities and alignment scores and derive the details for the case where rates are assumed to vary according to a gamma distribution.
Kelly
C
C
Department of Computer Science and Statistics, University of Rhode Island, Kingston 02881, USA.
Churchill
G A
GA
D000595:Amino Acid Sequence; D015982:Bias; D001699:Biometry; D019143:Evolution, Molecular; D008390:Markov Chains; D008962:Models, Theoretical; D011506:Proteins / Q000737:chemistry / Q000235:genetics; D016415:Sequence Alignment / Q000706:statistics & numerical data*
D011506:Proteins