As biologists began to use DNA sequences to infer phylogenies in the 1970s and 1980s, they developed explicit mathematical models describing how DNA sequences change over time. These models account for multiple changes at a given position in a DNA sequence. They also take into account different rates of change at different positions in a gene, at different positions in a codon, and among different nucleotides. For example, transitions (changes between two purines or between two pyrimidines) are usually more likely than are transversions (changes between a purine and pyrimidine).
Mathematical models can be used to compute how a tree might evolve given the observed data. A maximum likelihood method will identify the tree that most likely produced the observed data, given the assumed model of evolutionary change. Maximum likelihood methods can be used for any kind of characters, but they are most often used with molecular data, for which explicit mathematical models of evolutionary change are easier to develop. The principal advantages of maximum likelihood analyses are that they incorporate more information about evolutionary change than do parsimony methods, and they are easier to treat in a statistical framework. The principal disadvantages are that they are computationally intensive and require explicit models of evolutionary change (which may not be available for some kinds of character change).