Soybean is the most important legume cropped worldwide and can highly benefit from the biological nitrogen fixation (BNF) process. Brazil is recognized for its leadership in the use of inoculants and two strains,
Bradyrhizobium japonicum
CPAC 15 (=SEMIA 5079) and
Bradyrhizobium diazoefficiens
CPAC 7 (=SEMIA 5080) compose the majority of the 70 million doses of soybean inoculants commercialized yearly in the country. We studied a collection of natural variants of these two strains, differing in properties of competitiveness and efficiency of BNF. We sequenced the genomes of the parental strain SEMIA 566 of
B. japonicum
, of three natural variants of this strain (S 204, S 340 and S 370), and compared with another variant of this group, strain CPAC 15. We also sequenced the genome of the parental strain SEMIA 586 of
B. diazoefficiens
, of three natural variants of this strain (CPAC 390, CPAC 392 and CPAC 394) and compared with the genome of another natural variant, strain CPAC 7. As the main genes responsible for nodulation (nod, noe, nol) and BNF (nif, fix) in soybean
Bradyrhizobium
are located in symbiotic islands, our objective was to identify genetic variations located in this region, including single nucleotide polymorphisms (SNPs) and insertions and deletions (indels), that could be potentially related to their different symbiotic phenotypes. We detected 44 genetic variations in the
B. japonicum
strains and three in
B. diazoefficiens
. As the
B. japonicum
strains have gone through a longer period of adaptation to the soil, the higher number of genetic variations could be explained by survival strategies under the harsh environmental conditions of the Brazilian Cerrado biome. Genetic variations were detected in genes enconding proteins such as a dephospho-CoA kinase, related to the CoA biosynthesis; a glucosamine-fructose-6-phosphate aminotransferase, key regulator of the hexosamine biosynthetic pathway; a LysR family transcriptional regulator related to nodulation genes; and NifE and NifS proteins, directly related to the BNF process. We suggest potential genetic variations related to differences in the symbiotic phenotypes.