Evolutionary relationship of different Bacillus sp. on the basis of merA gene was inferred using the Neighbor-Joining method with the sum of branch length 0.53623375. The branch lengths and the evolutionary distances used to infer the phylogenetic tree are in the same units with the bootstrap test value 500. The analysis involved 17 nucleotide sequences and codon positions included were 1st+2nd+3rd+Noncoding. A total of 1279 positions were in the final dataset with no gaps and missing data. Evolutionary analyses were performed by MEGA 5 software (Tamura et al., 2004).