The assembled contigs were firstly annotated by NCBI-BlastN based on the local database with an e-value < 1e−6. Then, the boundaries of each gene were confirmed by Mitofy and exported as Sequin formatted files. tRNA genes were further predicted by tRNAscan (http://lowelab.ucsc.edu/tRNAscan-SE/)(Supplementary Table 43).