Molecular diagnosis of complex dual genome mitochondrial disorders is a challenge. It requires the identification of deleterious mutations in one of the~1,500 nuclear genes and the mitochondrial genome. If the molecular defect is in the mitochondrial genome, quantification of degree of mutation load (heteroplasmy) in affected tissues is important. Due to the extreme clinical and genetic heterogeneity, conventional sequence analysis of the candidate genes one-byone is impractical, if not impossible. The newly developed massively parallel next generation sequencing (NGS) technique, that allows simultaneous sequence analysis of multiple target genes, when appropriately validated with deep coverage and proper quality controls, can be used as an effective comprehensive diagnostic approach in CLIA certified clinical laboratories.