Each year, millions of people die from the airborne infectious illness tuberculosis (TB). Several drug-susceptible (DS) and drug-resistant (DR) forms of the causative agent, Mycobacterium tuberculosis (MTB), are currently common in the majority of affluent and developing nations, particularly in Bangladesh, and completely drug-resistant strains are beginning to arise. In this study, we developed fractional order tuberculosis model with fractal fractional operator by using generalized mittag-leffler kernel. By demonstrating characteristics such the boundedness of solutions, positivity, and the reliance of the solution on the original data, the biological well-posedness of the mathematical model formulation was investigated for TB cases of 2002 to 2017 in KPK Pakistan. Ulam-Hyres stability is also used to assess both local and global aspects of TB bacterial infection. Sensitivity analysis of the TB model with therapy was also examined. The advanced numerical technique is used to find the solution of the fractional-order system to check the impact of fractional parameters. Simulation highlights that all classes have converging qualities and retain established positions with time which shows the actual behavior of bacterial infection of TB. By examining the non-integer order with integer orders we find a better comparative result to maintain the position of non-integer order.