OBJECTIVE: To explore the association between the variant mutations within embB and ubiA, and the degree of ethambutol (EMB) resistance of Mycobacterium tuberculosis (M. tuberculosis) isolates. METHODS: A total of 146 M. tuberculosis isolates were used to determine the minimum inhibitory concentrations (MICs) of EMB with a 96-well microplate-based assay. The mutations within embB and ubiA among these isolates were identified with DNA sequencing. Moreover, a multivariate regression model and a computer model were established to assess the effects of mutations on EMB resistance. RESULTS: Our data showed that overall 100 isolates exhibited 28 mutated patterns within the sequenced embB and ubiA. Statistical analysis indicated that embB mutations Met306Val, Met306Ile, Gly406Ala, and Gln497Arg, were strongly associated with EMB resistance. Of these mutations, Met306Val and Gln497Arg were significantly associated with high-level EMB resistance. Almost all multiple mutations occurred in high-level EMB-resistant isolates. Although the mutation within ubiA accompanied with embB mutation presented exclusively in EMB-resistant isolates, four single ubiA mutations (Ala39Glu, Ser173Ala, Trp175Cys, and Val283Leu) leading to protein instability were observed in EMB-susceptible isolates. CONCLUSION: This study highlighted the complexity of EMB resistance. Some individual mutations and multiple mutations within embB and ubiA contributed to the different levels of EMB resistance.