The genome of a G. elata individual was sequenced using a whole-genome shotgun (WGS) approach (Supplementary Table 1). Through K-mer distribution analysis, the genome size was estimated to be 1.18 Gb (Supplementary Fig. 1). The assembly consisted of 3779 scaffolds, with a scaffold N50 of 4.9 Mb (total length = 1061.09 Mb) and contig N50 of 68.9 kb (total length = 1025.5 Mb) (Supplementary Table 2). Overall, 98.51% of the raw sequence reads could be mapped to the assembly, suggesting that our assembly results contained comprehensive genomic information (Supplementary Table 3). Gene region completeness was evaluated by RNA-Seq data (Supplementary Table 4): of the 80,646 transcripts assembled by Trinity, 98.66% could be mapped to our genome assembly, and 94.41% were considered as complete (more than 90% of the transcript could be aligned to one continuous scaffold). The completeness of gene regions was further assessed using CEGMA (conserved core eukaryotic gene mapping approach): 239 of 248 (96.37%) conserved core eukaryotic genes from CEGMA were captured in our assembly, and 217 (87.5%) of these were complete (Supplementary Table 5).