The intention of this work is to study a mathematical model for fractal-fractional tuberculosis and COVID-19 co-infection under the Atangana-Baleanu fractal-fractional operator. Firstly, we formulate the tuberculosis and COVID-19 co-infection model by considering the tuberculosis recovery individuals, the COVID-19 recovery individuals, and both disease recovery compartment in the proposed model. The fixed point approach is utilized to explore the existence and uniqueness of the solution in the suggested model. The stability analysis related to solve the Ulam-Hyers stability is also investigated. This paper is based on Lagrange’s interpolation polynomial in the numerical scheme, which is validated through a specific case with a comparative numerical analysis for different values of the fractional and fractal orders.