Morphological and molecular characterization of forty soybean (Glycine max (L.) Merr.) accessions were assessed using 7 morphological characters and 10 Randomly Amplified DNA Polymorphism (RAPD) primers, respectively. The field experiment was carried out at two locations using randomized complete block design with three replications. The aims of the research were to analyze the morphological and molecular organization of the existing diversity and to compare the genetic relatedness among forty soybean genotypes. The data were subjected to analysis of variance and correlation analysis to determine the extent of genetic variability and correlation coefficients among the characters, respectively. Principal Component Analysis (PCA) and Single Linkage Cluster Analysis (SLCA) were employed to group the accessions. The genetic relatedness among the accessions based on RAPD molecular markers was also presented in the form of a dendrogram generated by cluster analysis using the Unweighted Pair Group method with Arithmetic Mean (UPGMA). The relative effectiveness of the RAPD markers compared to botanical descriptors in assessing the diversity among the genotypes was investigated. The single linkage cluster technique classified the 40 accessions into seven clusters while the FASTCLUS technique revealed that the number of pods per plant, pod yield per plant, 100-seed weight and seed yield per plot contributed the largest proportion of morphological variation. Out of the 100 bands generated by the 10 primers, 31 were monomorphic and 69 polymorphic. The size of the fragment varied from 250 bp to 3000 bp. RAPDs markers were highly polymorphic and more discriminatory and informative as they were able to differentiate more pairs of genotypes than the botanical descriptors. The highest yield was recorded for TGx 1834-1E (477.60g/plot) and TGx 1910-2F (459.55).