In molecular epidemiological studies of Giardia intestinalis, an pathogenic intestinal flagellate, due to the presence of allelic sequence heterogeneity (ASH) on the tetraploid genome, the image of haplotype diversity in the field remains uncertain. Here we employed the nine assemblage B positive stool samples, which had previously reported from Kenyan children, for the clonal sequence analysis of multiple gene loci (glutamate dehydrogenase (GDH), triosephosphate isomerase (TPI), and beta-giardin (BG)). The diversified unique assemblage B haplotypes as GDH (n = 67), TPI (n = 84), and BG (n = 62), and the assemblage A haplotypes as GDH (n = 7), TPI (n = 14), and BG (n = 15), which were hidden in the previous direct-sequence results, were detected. Among the assemblage B haplotypes, Bayesian phylogeny revealed multiple statistically significant clusters (9, 7, and 7 clusters for GDH, TPI, and BG, respectively). A part of the clusters (2 for GDH and 1 for BG), which included >4 haplotypes from an individual sample, indicated the presence of co-transmission with multiple strains sharing a recent ancestor. Locus-dependent discrepancies, such as different compositions of derived samples in clusters and different genotyping results for the assemblages, were also observed and considered to be the traces of both intra- and inter-assemblage genetic recombination respectively. Our clonal sequence analysis for giardial population, which applied firstly in Kenya, could reveal the higher rates of ASH far beyond the levels reported in other areas and address the complex population structure. The clonal analysis is indispensable for the molecular field study of G. intestinalis.