The ML phylogenetic reconstruction was implemented using RAxML38. Nucleotide phylogenies were reconstructed using the GTR (General Time Reversible) nucleotide substitution model, while the core protein coding phylogeny was generated using Jones-Taylor-Thornton (JTT) amino acid substitution model. For both substitution models the rate of heterogeneity was estimated with 4 discrete rate categories. Internal branches of the phylogeny were estimated with 100 bootstraps.