In this study, we formulated and analyzed a deterministic mathematical model for the co-infection of COVID-19 and tuberculosis, to study the co-dynamics and impact of each disease in a given population. Using each disease's corresponding reproduction number, the existence and stability of the disease-free equilibrium were established. When the respective threshold quantities , and are below unity, the COVID-19 and TB-free equilibrium are said to be locally asymptotically stable. The impact of vaccine (i.e., efficacy and vaccinated proportion) and the condition required for COVID-19 eradication was examined. Furthermore, the presence of the endemic equilibria of the sub-models is analyzed and the criteria for the phenomenon of backward bifurcation of the COVID-19 sub-model are presented. To better understand how each disease condition impacts the dynamics behavior of the other, we investigate the invasion criterion of each disease by computing the threshold quantity known as the invasion reproduction number. We perform a numerical simulation to investigate the impact of threshold quantities with respect to their invasion reproduction number, co-infection transmission rate , and each disease transmission rate on disease dynamics. The outcomes established the necessity for the coexistence or elimination of both diseases from the communities. Overall, our findings imply that while COVID-19 incidence decreases with co-infection prevalence, the burden of tuberculosis on the human population increases.