In this research, the relationships among the 31 microsatellite markers with charcoal rot disease resistance related indices in 130 different soybean cultivars and lines were evaluated using association analysis based on the general linear model (GLM) and the mixed linear model (MLM) by the Structure and Tassel software. The results of microsatellite markers showed that the genetic structure of the studied population has three subpopulations (K=3) which the results of bar plat also confirmed it. In association analysis based on GLM and MLM models, 31 and 35 loci showed significant relationships with the evaluated traits, respectively, and confirmed considerable variation of the studied traits. The identified markers related to some of the studied traits were the same which can probably be due to pleiotropic effects or tight linkage among the genomic regions controlling these traits. Some of these relationships were including, the relationship between Sat_252 marker with amount of charcoal rot disease, Satt359, Satt190 and Sat_169 markers with number of microsclerota in stem, amount of charcoal rot disease and severity of charcoal rot disease, Sat_416 marker with number of microsclerota in stem and amount of charcoal rot disease and the Satt460 marker with number of microsclerota in stem and severity of charcoal rot disease. The results of this research and the linked microsatellite markers with the charcoal rot disease-related characteristics can be used to identify the suitable parents and to improve the soybean population in future breeding programs.