Mycobacterium tuberculosis transmission routes can be estimated from genotypic analysis of clinical isolates from patients. In Japan, still a middle-incidence country of TB, a unique genotype strain designated as 'M-strain' has been isolated nationwide recently. To ascertain the history of the wide spread of the strain, 10 clinical isolates from different areas were subjected to genome-wide analysis based on deep sequencers. Results show that all isolates possessed common mutations to those of referential strains. The greatest number of accumulated single nucleotide variants (SNVs) from the oldest coalescence was 13 nucleotides, indicating high clonality of these isolates. When an SNV common to the isolates was used as a surrogate marker of the clone, authentic clonal isolates with variation in a reliable subset of variable number of tandem repeat (VNTR) genotyping method can be selected successfully from clinical isolates populations of M. tuberculosis. When the authentic clones can also be assigned to sub-clonal groups by SNVs derived from the genomic comparison, they are classifiable into three sub-clonal groups with a bias of geographical origins. Feedback from genomic analysis of clinical isolates of M. tuberculosis to genotypic markers will be an efficient strategy for the big data in various settings for public health actions against TB.