Toxoplasma gondii consists of three genotypes, namely genotype I, II and III. Based on its virulence, T. gondii can be divided into virulent and avirulent strains. This study intends to evaluate an alternative method for predicting T. gondii virulence using hierarchical cluster analysis based on complete coding sequences (CDS) of sag1, gra7 and rop18 genes. Dendrogram was constructed using UPGMA with a Kimura 80 nucleotide distance measurement. The results showed that the prediction errors of T. gondii virulence using sag1, gra7 and rop18 were 7.41%, 6.89% and 9.1%, respectively. Analysis based on CDS of gra7 and rop18 was able to differentiate avirulent strains into genotypes II and III, whereas sag1 failed to differentiate.