Pairwise genetic distances for COI and 16S were calculated in Geneious
(Biomatters Ltd). For 16S distances were calculated on the alignment
after trimming with Gblocks. The best-fitting model of evolution for
each gene was found using JModelTest 2.1.4 For all genes the GTRþI þGmodelwas
considered the best fit according to the Akaike Information Criterion,
but due to statistical concerns regarding the coestimation of the gamma and invariant-site parameters (discussed in the RAxML manual; theGTRþG model was chosen instead.