This study proposes a mathematical model for examining the COVID-19 and tuberculosis (TB) co-dynamics thoroughly. First, the single infection dynamics: COVID-19 infection and TB infection models are taken into consideration and examined. Following that, the co-dynamics with TB and COVID-19 is also investigated. In order to comprehend the developed model dynamics, the basic system attributes including the region of definition, theory of nonnegativity and boundedness of solution are investigated. Further, a qualitative analysis of the equilibria of the formulated model equations is performed. The equilibria of both infection models are globally asymptotically stable if their respective basic reproductive number is smaller than one. As the associated reproductive number reaches unity, they experience the forward bifurcation phenomenon. Additionally, it is demonstrated that the formulated co-dynamics model would not experience backward bifurcation by applying the center manifold theory. Moreover, model fitting is done by using daily reported COVID-19 cumulative data in Ethiopia between March 13, 2020, and May 31, 2022. For instance, the non-linear least squares approach of fitting a function to data was performed in the fitting process using scipy.optimize.curve_fit from the Python. Finally, to corroborate the analytical findings of the model equation, numerical simulations were conducted.