This study investigates a symmetric fractional-order epidemic model with time delays and non-monotonic incidence rates, considering two viral strains. By confirming the existence, uniqueness, and boundedness of the system’s solutions, the research ensures the model’s well-posedness, guaranteeing its mathematical soundness and practical relevance. The study calculates and evaluates the equilibrium points and the basic reproduction numbers R01 and R02 to understand the dynamic behavior of the model under different parameter settings. Through the application of the Lyapunov method, the research examines the asymptotic global stability of the system, determining whether it will converge to a particular equilibrium state over time. Furthermore, Hopf bifurcation theory is employed to investigate potential periodic solutions and bifurcation scenarios, highlighting how the system might shift from stability to periodic oscillations under certain conditions. By utilizing the Adams-Bashforth-Moulton numerical simulation method, the theoretical results are validated, reinforcing the conclusions and demonstrating the model’s applicability in real-world contexts. It emphasizes the importance of fractional-order models in addressing epidemiological issues related to time delays (τ), individual heterogeneity (m, k), and memory effects (θ), offering greater accuracy compared with traditional integer-order models. In summary, this research provides a theoretical foundation and practical insights, enhancing the understanding and management of epidemic dynamics through fractional-order epidemic models.