The minimum bias differential yield is reported for KS0,lambda,lambda(, )(sic)(-) , (sic)(+)omega(-),omega(+)and phi strange hadrons, obtained from Au - Au collisions at RHIC energies of root sNN = 11.5, 19.6, 27 and 39 GeV as a function of pT in kinematic range of 0 < pT < 5 GeV and pseudora-pidity region of |eta | < 0.5 using PYTHIA8 and HIJING2 event generators. The PYTHIA8 is tuned to different Colour Reconnection (CR) modes, and results are compared with STAR data at RHIC. PYTHIA8 predictions agree with STAR data at low pT values, while underestimate at high pT val-ues.The observed deviations in C R 0 are due to multiparton interactions, the deviations in C R 1 are due to QCD effects. At low pT, the non-pQCD effects are dominant, while at high pT, the pQCD effects are dominant. The deviations in C R 2 mode is due to the minimization of strings followed by different colour combinations, and the deviations in case of C R 3 mode is connected with strings as viewed diagonally in extended bags. HIJING2 predictions explain the STAR data very well over the entire pT region.