In this article, mathematical modeling and stability analysis of Lyme disease and its transmission dynamics using Caputo fractional-order derivatives is presented. The compartmental model has been formulated to analyze the spread of Borrelia burgdorferi virus through tick vectors and mammalian hosts. The feasible region is established, and the boundedness of the model is verified. Analytically, the disease-free equilibrium and the basic reproduction number $\left({\Re }_0\right)$ has been determined to assess outbreak potential. By virtue of the fixed-point theory, the existence and uniqueness of solutions has been established. The numerical simulations are obtained via the Runge–Kutta 4 method, demonstrating the model’s ability to capture realistic disease progression. Finally, sensitivity analysis and control strategies (tick population reduction, host vaccination, public awareness, and early treatment) are evaluated, revealing that integrated control measures significantly reduce infection rates and enhance recovery.