Hepatocellular carcinoma (HCC) is the third leading cause of cancer mortality worldwide. In this study, we had analysed the copy number variations and heteroplasmic mutations of mitochondria (MT) in 88 HCC individuals. The average copy number of MT genome in normal samples was significantly greater than that in tumor samples. Overall, the number of heteroplasmic mutations in 88 tumor and their matched normal samples were 241 and 173, respectively. There was higher positive ratio of heteroplasmic mutations in tumor samples (86%) than normal samples (73%). Worthwhile mention, ND1 gene harbored greater mutation frequency and more nonsynonymous mutations in tumor samples. Interestingly, 202 tumor-specific heteroplasmic mutations were detected. Moreover, ND1, ND3, ND4, ND5 and ND6 genes had higher ratio of nonsynonymous versus synonymous mutations in tumor-specific heteroplasmic mutations. It might suggest that the disorder of NADH dehydrogenase (complex I) resulted by heteroplasmic mutations may have close relation with tumorigenesis of hepatocellular carcinoma. This study provided theoretical basis for further understanding mechanism of tumorigenesis from the perspective of mitochondrial heteroplasmic mutations.