design: Up to five specimens of investigated species were reserved and macerated together with 360 mL of phosphate-buffered saline, from which 140 mL of the supernatant was reserved for total RNA extraction and quantification procedures. The extraction products were then standardized to a concentration of 0.2 ng/L, fragmented, and labeled with adapter sequences. Subsequently, the genomic libraries obtained were quantified and sequenced by the paired-end method (2 x 150) on the NextSeq 500 System platform. The obtained sequencing files were demultiplexed and subjected to quality assessment using Fastp v.0.23.4, configured to remove adapter sequences and reads with base quality Phred < 20 and length less than 50 nt. The contigs were assembled by the de novo method using MEGAHIT v.1.2.9 in its standard configuration (k-mers of 21, 29, 39, 59, 79, 99, 119, and 141 nt). The contig corresponding to the mitochondrial genome of the species was identified using DIAMOND v.2.1.9.163 (Blastx modality, e-value of 0.00001) and manually inspected using Geneious v.11.1.5.