Genome quality evaluation
After the genome assembly process, several genome quality evaluation methods were used. We used BUSCO (v2.0) software (Simao et al. 2015) to estimate the percentage of conserved orthologs in the assembled genome. The conserved gene sets of Eukaryota and Metazoa were used as a database and the complete BUSCO genes, fragment genes, and missing genes were detected. BWA (Li& Durbin 2009) and BLAT (Kent 2002) were used for the mapping ratio of NGS data and de novo assembled transcripts, respectively.