We present a genome-wide comparative and comprehensive analysis of three different 12 sequencing methods (conventional next generation sequencing (NGS), tag-based single strand 13 sequencing (eg. SSCS), and Duplex Sequencing for investigating mitochondrial mutations in 14 human breast epithelial cells. Duplex Sequencing produces a single strand consensus sequence 15 (SSCS) and a duplex consensus sequence (DCS) analysis, respectively. Our study validates that 16 although high-frequency mutations are detectable by all the three sequencing methods with the 17 similar accuracy and reproducibility, rare (low-frequency) mutations are not accurately detectable 18 by NGS and SSCS. Even with conservative bioinformatical modification to overcome the high error 19 rate of NGS, the NGS frequency of rare mutations is 7.0x10 -4 . The frequency is reduced to 1.3x10 -4 20 with SSCS and is further reduced to 1.0x10 -5 using DCS. Rare mutation context spectra obtained 21 from NGS significantly vary across independent experiments, and it is not possible to identify a 22 dominant mutation context. In contrast, rare mutation context spectra are consistently similar in all 23 independent DCS experiments. We have systematically identified heat-induced artifactual 24 mutations and corrected the artifacts using Duplex Sequencing. All of these artifacts are 25 stochastically occurring rare mutations. C>A/G>T, a signature of oxidative damage, is the most 26 increased (170-fold) heat-induced artifactual mutation type. Our results strongly support the claim 27 that Duplex Sequencing accurately detects low-frequency mutations and identifies and corrects 28 artifactual mutations introduced by heating during DNA preparation. 29 Sequence (SSCS), Next-Generation Sequencing (NGS), sequencing error, rare mutations, oxidative 31 DNA damage, heat-induced mutations, mitochondrial DNA, human breast cells 32 33 1. Introduction 34 Next-generation sequencing (NGS) has rapidly transformed entire areas of basic research and 35 therapeutic applications by making large scale genomic studies feasible through reduced cost and 36faster turnaround time [1,2]. NGS has been extensively used to study clonal (high-frequency) 37 mutations, but not subclonal (low-frequency) mutations. A major impediment in investigating 38 subclonal (low-frequency) mutations is that conventional NGS methods have high error rates (10 -2 to 39 10 -3 ), which obscure true mutations that occur less frequently than errors [3,4]. These subclonal 40 mutations may account for the genetic heterogeneity of tumors and tumor recurrence, as well as 41 provide a reservoir for the rapid development of resistance to chemotherapy [5]. 42 Conventional sequencing technologies sequence only a single strand of DNA. In contrast, 43 Duplex Sequencing examines both strands of DNA and scores mutations only if they are present on 44 both strands of the same DNA molecule as complementary substitutions. This significantly reduces 45 sequencing error rates to < 5 x 10 -8 [6-9]. In the first report of Duplex Seq...