Figure 4

Phylogenetic tree from Bayesian analysis of nucleotide sequences for 12 protein-coding genes with third codon positions removed (7,818 characters). The best-fit substitution model for each of 12 genes was estimated using the AIC criterion implemented in MrModeltest 2.3 [58]. 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 one-third of the trees was discarded as burn-in (see methods for analysis details).