Maximum likelihood estimation of the phylogenetic relationships based on the mitochondrial COI barcode region (K2P model, discrete Gamma distribution for rate differences with 3 categories (+G, parameter = 0.0500)). Nucleotide positions with less than 98% site coverage were eliminated, resulting in 570 analysed positions. Numbers of major nodes indicate bootstrap values from 1000 pseudoreplicates from the NJ and ML method.