Phylogenetic relationships of the Blastocystis sp. subtypes identified in the present study and other reference subtypes. The phylogeny was inferred with both neighbor-joining (NJ) and maximum likelihood (ML) analysis. Evolutionary distances were computed using the maximum composite likelihood method with rate variation among sites modeled using a gamma distribution (shape parameter = 0.5). The percentage of trees that clustered together based on a bootstrap test (1000 replicates) and posterior probabilities (expressed as a percentage) are shown beside the branch nodes. An asterisk indicates a value of less than 50%; if both analyses produced values lower than 50%, no values are shown for that node. The genotypes detected in present study in bold are shown with filled triangles.