As the SARS-CoV-2 pandemic progressed, distinct variants emerged and dominated in England. These variants, Wildtype, Alpha, Delta, and Omicron were characterized by variations in transmissibility and severity. We used a robust mathematical model and Bayesian inference framework to analyse epidemiological surveillance data from England. We quantified the impact of non-pharmaceutical interventions (NPIs), therapeutics, and vaccination on virus transmission and severity. Each successive variant had a higher intrinsic transmissibility. Omicron (BA.1) had the highest basic reproduction number at 8.1 (95% credible interval (CrI) 6.8-9.3). Varying levels of NPIs were crucial in controlling virus transmission until population immunity accumulated. Immune escape properties of Omicron decreased effective levels of protection in the population by a third. Furthermore, in contrast to previous studies, we found Alpha had the highest basic infection fatality ratio (2.8%, 95% CrI 2.3-3.2), followed by Delta (2.0%, 95% CrI 1.5-2.4), Wildtype (1.2%, 95% CrI 1.0-1.3), and Omicron (0.6%, 95% CrI 0.4-0.8). Our findings highlight the importance of continued surveillance. Long-term strategies for monitoring and maintaining effective immunity against SARS-CoV-2 are critical to inform the role of NPIs to effectively manage future variants with potentially higher intrinsic transmissibility and severe outcomes.