Aligner version 3

Total genomic DNA was isolated using QIAamp DNA Mini Kit (Qiagen, Courtaboeuf, France) as per the manufacturer's instructions. The 5′ region of the mitochondrial cytochrome c oxidase subunit I (COI) gene was amplified with the primers TelF1 (5′-TCGACTAATCAYAAAGAYATYGGCAC-3′) and TelR1 (5′-ACTTCTGGGTGNCCAAARAATCARAA-3′) (Dettaï et al. 2011) . PCR reactions were performed in 20 μl, containing 1 ng of DNA, 1× CoralLoad PCR buffer, 3 mM MgCl2, 66 μM of each dNTP, 0.15 μM of each primer, and 0.5 units of Taq DNA polymerase (Qiagen). The amplification protocol was 4 min at 94 °C, followed by 40 cycles at 94 °C for 30 s, 48 °C for 40 sec, and 72 °C for 50 s, with a final extension at 72 °C for 7 min. PCR products were purified (Ampure XP Kit, Beckman Coulter) and sequenced in both directions on a 3730xl DNA Analyzer 96-capillary sequencer (Applied Biosystems, Foster City, USA). We used the CodonCode Aligner version 3.7.1 software (Codon Code Corporation, Dedham, MA, USA) to edit sequences, compared them to the GenBank database content with BLAST, and deposited them in GenBank under accession numbers MW484937-MW484939. Species identification was confirmed with the BOLD identification engine (Ratnasingham and Hebert 2007) (link).

Total genomic DNA was isolated using QIAamp DNA Mini Kit (Qiagen, Courtaboeuf, France)
as per the manufacturer's instructions. The 5 ʹ region of the mitochondrial cytochrome c oxidase subunit 1 (cox1) gene was amplified with the primers FishF1 (5 ʹ -TCA ACC AAC CAC AAA GAC ATT GGC AC-3 ʹ ) and FishR1 (5 ʹ -TAG ACT TCT GGG TGG CCA AAG AAT CA-3 ʹ ) (Ward et al., 2005) . PCR reactions were performed in a total volume of 20 µl, containing 1 ng of DNA, 1× CoralLoad PCR buffer, 3 mM MgCl 2 , 66 µM of each dNTP, 0.15 µM of each primer and 0.5 units of Taq DNA polymerase (Qiagen). The amplification protocol was 4 min at 94°C, followed by 40 cycles at 94°C for 30 s, 48°C for 40 s, and 72°C for 50 s, with a final extension step at 72°C for 7 min. PCR products were purified (Ampure XP Kit, Beckman Coulter, Brea, USA) and sequenced in both directions on a 3730xl DNA Analyzer 96-capillary sequencer (Applied Biosystems, Foster City, USA). We used CodonCode Aligner version 3.7.1 software (CodonCode Corporation, Dedham, MA, USA) to edit the sequence, which was 652 bp in length, compared it to the GenBank database content with BLAST, and deposited it in GenBank under the accession number MK275650. Species identification was confirmed with the BOLD identification engine (Ratnasingham & Hebert, 2007) .

Total genomic DNA was isolated using QIAamp DNA Mini Kit (Qiagen) as per the manufacturer's instructions. The 5′ region of the mitochondrial cytochrome c oxidase subunit I (COI) gene was amplified with the primers FishF1 (5′-TCAACCAACCACAAAGACATTGGCAC-3′) and FishR1 (5′-TAGACTTCTGGGTGGCCAAAGAATCA-3′) [19] . PCR reactions were performed in 20 μl, containing 1 ng of DNA, 1× CoralLoad PCR buffer, 3 mM MgCl2, 66 μM of each dNTP, 0.15 μM of each primer, and 0.5 units of Taq DNA polymerase (Qiagen). The amplification protocol was 4 min at 94°C, followed by 40 cycles at 94°C for 30 sec, 48°C for 40 sec, and 72°C for 50 sec, with a final extension at 72°C for 7 min.
PCR products were purified (Ampure XP Kit, Beckman Coulter) and sequenced in both directions on a 3730xl DNA Analyzer 96-capillary sequencer (Applied Biosystems). We used CodonCode Aligner version 3.7.1 software (CodonCode Corporation, Dedham, MA, USA) to edit sequences, which were 652 bp in length, compared them to the GenBank database content with BLAST, and deposited them in GenBank under accession numbers KX926437 -KX926442. Species identification was confirmed with the BOLD identification engine [20] .