The neutron capture gamma-ray spectra for 8 nuclides, (89)y, Nb-93, 1271, Cs-133, Pr-141, Au-197, (TI)-T-nat, and Bi-209, 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.