  This paper is concerned with the evolution of haploid organisms that
reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed
the quasispecies model in an attempt to understand such an evolutionary
process. Their work has impacted antiviral treatment and vaccine design
strategies. Yet, predictions of the quasispecies model are at best viewed as a
guideline, primarily because it assumes an infinite population size, whereas
realistic population sizes can be quite small. In this paper we consider a
population genetics-based model aimed at understanding the evolution of such
organisms with finite population sizes and present a rigorous study of the
convergence and computational issues that arise therein. Our first result is
structural and shows that, at any time during the evolution, as the population
size tends to infinity, the distribution of genomes predicted by our model
converges to that predicted by the quasispecies model. This justifies the
continued use of the quasispecies model to derive guidelines for intervention.
While the stationary state in the quasispecies model is readily obtained, due
to the explosion of the state space in our model, exact computations are
prohibitive. Our second set of results are computational in nature and address
this issue. We derive conditions on the parameters of evolution under which our
stochastic model mixes rapidly. Further, for a class of widely used fitness
landscapes we give a fast deterministic algorithm which computes the stationary
distribution of our model. These computational tools are expected to serve as a
framework for the modeling of strategies for the deployment of mutagenic drugs.
