Exact computation of pattern probabilities in random sequences generated by Markov chains
Central Institute of Molecular Biology, Academy of Sciences of the GDR 1115 Berlin Buch, Robert Rössle Strasse 10, GDR
Observed patterns in macromolecular sequences are often considered as words and compared with their probabilities of occurring in random sequences. Calculation of these probabilities, however, often lacks rigour. We have developed an algorithm for exact computation of such probabilities for stochastic sequences that follow a Markov chain model. The method is applicable to the case that a random sequence contains one out of two given patterns P and Q, or both simultaneously. Another application yields the probability Junction P(x) that a sequence contains pattern P exactly x times. An application to patterns that include wild-card characters yields probabilities for homonucleotide clusters of a given length. We prove the probability of multiple runs of single nucleotides in the SV40 genome to be in accordance with the dinucleotide composition of the sequence, although it is in conflict with mononucleotide composition.
Received on January 10, 1990; accepted on April 23, 1990
This article has been cited by other articles:
![]() |
J. R. Stone and G. A. Wray Rapid Evolution of cis-Regulatory Sequences via Local Point Mutations Mol. Biol. Evol., September 1, 2001; 18(9): 1764 - 1770. [Abstract] [Full Text] [PDF] |
||||
