A new mathematical model incorporating epidemiological features of the co-dynamics of tuberculosis (TB) and SARS-CoV-2 is analyzed. Local asymptotic stability of the disease-free and endemic equilibria are shown for the sub-models when the respective reproduction numbers are below unity. Bifurcation analysis is carried out for the TB only sub-model, where it was shown that the sub-model undergoes forward bifurcation. The model is fitted to the cumulative confirmed daily SARS-CoV-2 cases for Indonesia from February 11, 2021 to August 26, 2021. The fitting was carried out using the fmincon optimization toolbox in MATLAB. Relevant parameters in the model are estimated from the fitting. The necessary conditions for the existence of optimal control and the optimality system for the co-infection model is established through the application of Pontryagin’s Principle. Different control strategies: face-mask usage and SARS-CoV-2 vaccination, TB prevention as well as treatment controls for both diseases are considered. Simulations results show that: (1) the strategy against incident SARS-CoV-2 infection averts about 27,878,840 new TB cases; (2) also, TB prevention and treatment controls could avert 5,397,795 new SARS-CoV-2 cases. (3) In addition, either SARS-CoV-2 or TB only control strategy greatly mitigates a significant number of new co-infection cases.