<p>The concatenated dataset comprising 278 amino acid sites of GAPDH, 280 amino acid sites of actin, 274 amino acid sites of EF-1α, and 1338 nucleotide sites of SSU rRNA gene sequences was analyzed in 28 parabasalian species. The tree was estimated in RAxML using separate models with the parameters and branch length optimized for each gene partitions individually. The supporting values (bootstrap in RAxML/Bayesian posterior probability) are indicated at the nodes. Values below 50% or 0.5 are indicated with hyphens. When the site-heterogeneous CAT model was used in each partition, the identical tree topology with similar bootstrap values was obtained (data not shown). Vertical bars to the right of the tree represent the parabasalian classes. The scale bar corresponds to 0.10 substitutions per site.</p
