Abstract Background The accumulation of mtDNA mutations in different tissues from various mouse models has been widely studied especially in the context of mtDNA mutation-driven ageing but has been confounded by the inherent limitations of the most widely used approaches. By implementing a method to sequence mtDNA without PCR amplification prior to library preparation, we map the full unbiased mtDNA mutation spectrum across six distinct brain regions from mice. Results We demonstrate that ageing-induced levels of mtDNA mutations (single nucleotide variants and deletions) reach stable levels at 50 weeks of age but can be further elevated specifically in the cortex, nucleus accumbens (NAc), and paraventricular thalamic nucleus (PVT) by expression of a proof-reading-deficient mitochondrial DNA polymerase, Polg^D181A. The increase in single nucleotide variants increases the fraction of shared SNVs as well as their frequency, while characteristics of deletions remain largely unaffected. In addition, Polg^D181A also induces an ageing-dependent accumulation of non-coding control-region multimers in NAc and PVT, a feature that appears almost non-existent in wild-type mice. Conclusions Our data provide a novel view of the spatio-temporal accumulation of mtDNA mutations using very limited tissue input. The differential response of brain regions to a state of replication instability provides insight into a possible heterogenic mitochondrial landscape across the brain that may be involved in the ageing phenotype and mitochondria-associated disorders.