The objective of this study were to conduct association mapping for drought tolerance at the seedling stage and yield-related traits. 60 cowpea accessions were used in the study. Single-nucleotide polymorphisms (SNPs) discovered through genotyping by sequencing (GBS) were used for genotyping. Association mapping was conducted using single-marker regression (SMR) in Q Gene, and general linear model (GLM) and mixed linear model (MLM) built in TASSEL. The population of the cowpea accessions were analysed using STRUCTURE 2.3.4 and the peak of delta K in the greenhouse showed seven population types, whereas the peak of delta K in the glasshouse indicated the presence of six population types. One SNP marker, 14083649|F|0-9 was associated with NP with a p value <0.001. Fifty SNP markers were associated with PWT at p <0.001. Four SNP markers, 14074781|F|0-16, 100047392|F|0-36, 14083801|F|0-28 and 100051488|F|0-49 were associated with AVSPD at p <0.001. SNP markers, 14074781|F|0-16, 14083801|F|0-28 and 100051488|F|0-49 were associated with PL at P <0.001. Five SNP markers, 100047392|F|0-36, 14083801|F|0-28, 100072738|F|0-34, 14076881|F|0-49 and 14076881|F|0-49 were associated with PWDTH at p <0.001. The 65 SNP markers identified can be used in cowpea molecular breeding to select for AVSPD, NP, PL, PWDTH, PWT, and RR through marker assisted selection (MAS).