The neutron capture gamma-ray spectra for 8 nuclides, 89 Y, 93 Nb, 127 I, 133 Cs, 141 Pr, 197 Au, nat Tl, and 209 Bi, were calculated by using the Hauser-Feshbach statistical model, and their results were compared with the available experimental data. Two dominant ingredients to perform the statistical calculation were the level densities described by the Gilbert-Cameron approach with an improved systematics, and the gamma-ray transmission coefficients described by gamma-ray strength functions. Although various gamma-ray strength functions with a Lorentzian formula have been developed by using the photonuclear data or a microscopic analysis, they have failed to reproduce the occasional anomalous bumps observed near or below a neutron binding. In this work, we could reproduce the bumps well by adding a Lorentzian with an energy-temperature dependent width into a giant electric dipole resonance with an enhanced generalized Lorentzian. In addition, we introduced a correction function so as to compensate for the shortcomings due to missing levels or level-cuts.