Figure 4

Phylogenetic tree from Bayesian analysis of nucleotide sequences for 12 protein-coding genes with third codon positions included (11,784 characters). The best-fit substitution model for each of 12 genes was estimated using the AIC criterion implemented in MrModeltest 2.3. The resulting best-fit model for each of 12 genes was then used for Bayesian analysis. Bayesian posterior probability values (BPP), shown above the nodes, were estimated after the initial 200 trees (the first 2x10⁵ generations) was discarded as burn-in (see Methods for analysis details).