Abstract
We construct a mathematical model for the polymerase chain reaction and its mutations using the theory of branching processes. Under this model we study the number of mutations in a randomly chosen sequence after n PCR cycles. A method for estimating the mutation is proposed and the variance of this estimator is studied. We also study the distribution of the Hamming distance between two randomly chosen sequences and a method for estimating the mutation rate based on pairwise differences is proposed.