Abstract
Multiple strains of the SARS-CoV-2 have arisen and jointly influence the trajectory of the coronavirus disease (COVID-19) pandemic. However, current models rarely account for this multi-strain dynamics and their different transmission rate and response to vaccines. We propose a new mathematical model that accounts for two virus variants and the deployment of a vaccination program. To demonstrate utility, we applied the model to determine the control reproduction number (Rc) and the per day infection, death and recovery rates of each strain in the US pandemic. The model dynamics predicted the rise of the alpha variant and shed light on potential impact of the delta variant in 2021. We obtained the minimum percentage of fully vaccinated individuals to reduce the spread of the variants in combination with other intervention strategies to deaccelerate the rise of a multi-strain pandemic.