We developed a TB-COVID-19 co-infection epidemic model using a non-linear dynamical system by subdividing the human population into seven compartments. The biological well-posedness of the formulated mathematical model was studied via proving properties like boundedness of solutions, no-negativity, and the solution's dependence on the initial data. We then computed the reproduction numbers separately for TB and COVID-19 sub-models. The criterion for stability conditions for stationary points was examined. The basic reproduction number of sub-models used to suggest the mitigation and persistence of the diseases. Qualitative analysis of the sub-models revealed that the disease-free stationary points are both locally and globally stable provided the respective reproduction numbers are smaller than unit. The endemic stationary points for each sub-models were globally stable if their respective basic reproduction numbers are greater than unit. In each sub-model, we performed an analysis of sensitive parameters concerning the corresponding reproduction numbers. Results from sensitivity indices of the parameters revealed that deceasing contact rate and increasing the transferring rates from the latent stage to an infected class of individuals leads to mitigating the two diseases and their co-infections. We have also studied the analytical behavior of the full co-infection model by deriving the equilibrium points and investigating the conditions of their stability. The numerical experiments of the proposed co-infection model agree with the findings in the analytical results.