500 Haemulidae phylogenies

These are the sample of 500 trees from the posterior distribution generated by BEAST that we used throughout the analyses. The nucleotide dataset was 2909 bp long and consisted of partial sequences of three mtDNA genes (16S rRNA, COI, and cytb) and two nuclear genes (RAG2 and S7 ribosomal protein intron 1). Relative divergence times of the sampled haemulid species were estimated using an uncorrelated lognormal (UCLN) model of molecular evolutionary rate heterogeneity implemented in the computer program BEAST v. 1.6.1 (Drummond et al. 2006; Drummond and Rambaut 2007). Each gene was treated as a separate data partition and for the protein coding genes (COI, cytb, and RAG1) we applied three partitions that corresponded to three codon positions.