Infectious diseases have long been a shaping force in human history, necessitating a comprehensive understanding of their dynamics. This study introduces a co-evolution model that integrates both epidemiological and evolutionary dynamics. Utilizing a system of differential equations, the model represents the interactions among susceptible, infected, and recovered populations for both ancestral and evolved viral strains. Methodologically rigorous, the model’s existence and uniqueness have been verified, and it accommodates both deterministic and stochastic cases. A myriad of graphical techniques have been employed to elucidate the model’s dynamics. Beyond its theoretical contributions, this model serves as a critical instrument for public health strategy, particularly predicting future outbreaks in scenarios where viral mutations compromise existing interventions.