Recent breakthroughs in the genetic manipulation of mitochondrial DNA (mtDNA) have enabled precise base substitutions and the efficient elimination of genomes carrying pathogenic mutations. However, reconstituting mtDNA deletions linked to mitochondrial myopathies remains challenging. Here, we engineered mtDNA deletions in human cells by co-expressing end-joining (EJ) machinery and targeted endonucleases. Using mitochondrial EJ (mito-EJ) and mito-ScaI, we generated a panel of clonal cell lines harboring a ∼3.5 kb mtDNA deletion across the full spectrum of heteroplasmy. Investigating these cells revealed a critical threshold of ∼75% deleted genomes, beyond which oxidative phosphorylation (OXPHOS) protein depletion, metabolic disruption, and impaired growth in galactose-containing media were observed. Single-cell multiomic profiling identified two distinct nuclear gene deregulation responses: one triggered at the deletion threshold and another progressively responding to heteroplasmy. Ultimately, we show that our method enables the modeling of disease-associated mtDNA deletions across cell types and could inform the development of targeted therapies.