Methods
We analyze a mtDNA dataset comprising 1,492 specimens from 598 locations worldwide. Molecular operational taxonomic units (MOTUs) ranged from 267 to 520, based on the Barcode Index Number (BIN), Bayesian Poisson tree processes (bPTP), multi-rate Poisson tree processes (mPTP), single-rate Poisson tree processes (sPTP), and generalized mixed Yule coalescent (sGMYC) approaches.