Fig. 4

Species delimitation results visualised as bars on an ultrametric Bayesian maximum clade credibility tree of A. torrentium COI gene. Next to the tree phylogroups obtained according to reconstruction of mtDNA and nuDNA are presented. Then follows two partitions of Automatic Barcode Gap Discovery (lumper and splitter) (ABGD [101]), General Mixed Yule Coalescent (GMYC [102]); Bayesian implementation of the Poisson Tree Processes (bPTP [140]); Multi-rate Poisson Tree Processes (mPTP [141]), and Templeton, Crandall and Sing method (TCS [129]). Also, in the last column phylogroups’ abbreviations are given: ZV Zeleni Vir, GK Gorski Kotar, ŽPB Žumberak, Plitvice and Bjelolasica, LD Lika and Dalmatia, BAN Banovina, SB southern Balkans, CSE central and southeastern Europe, APU Apuseni. As outgroups Austropotamobius pallipes, Astacus astacus and Pontastacus leptodactylus were used