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.