This work is devoted to the investigation of virus quasi-species evolution and diversification due to mutations, competition for host cells, and cross-reactive immune responses. The model consists of a nonlocal reaction-diffusion equation for the virus density depending on the genotype considered to be a continuous variable and on time. This equation contains two integral terms corresponding to the nonlocal effects of virus interaction with host cells and with immune cells. In the model, a virus strain is represented by a localized solution concentrated around some given genotype. Emergence of new strains corresponds to a periodic wave propagating in the space of genotypes. The conditions of appearance of such waves and their dynamics are described. © 2020 by the authors.