Neighbour-joining topology of the barcoding region of our dataset of the analysed species based on Kimura 2-parameter distances. Triangles show the relative number of individuals sampled (height) and sequence divergence (width). Pale blue colouration indicates species associated with BINs that have a maximum within-BIN distance > 2.2%. Dark blue colouration indicates species associated with BINs that have a higher within-BIN p-distance compared to NN p-distance. Orange colouration indicates species maximum intraspecific distances > 2.2% (cases of “BIN-sharing”). Numbers next to nodes represent non-parametric bootstrap values > 90% (1,000 replicates). Sequences shorter than 400 bp were excluded from this analysis. The aligned sequences and N-J topology can be reviewed in Suppl. materials 3 and 4.