Extraordinary variation has been found in mitochondrial (mt) genome inheritance, gene content and arrangement among bivalves. However, only few bivalve mt genomes have been comparatively analyzed to infer their evolutionary scenarios. In this study, the complete mt genome of the venerid Paphia euglypta (Bivalvia: Veneridae) was firstly studied and, secondly, it was comparatively analyzed with other venerids (e.g., Venerupis philippinarum and Meretrix petechialis) to better understand the mt genome evolution within a family. Though several common features such as the AT content, codon usage of protein-coding genes, and AT/GC skew are shared by the three venerids, a high level of variability is observed in genome size, gene content, gene order, arrangements and primary sequence of nucleotides or amino acids. Most of the gene rearrangement can be explained by the "tandem duplication and random loss" model. From the observed rearrangement patterns, we speculate that block interchange between adjacent genes may be common in the evolution of mt genomes in venerids. Furthermore, this study presents several new findings in mt genome annotation of V. philippinarum and M. petechialis, and hence we have reannotated the genome of these two species as: (1) the ORF of the formerly annotated cox2 gene in V. philippinarum is deduced by using a truncated "T" codon and a second cox2 gene is identified; (2) the trnS-AGN gene is identified and marked in the mt genome of both venerids. Thus, this study demonstrated a high variability of mt genomes in the Veneridae, and showed the importance of comparative mt genome analysis to interpret the evolution of the bivalve mt genome.