We study numerically the statistical properties of high-density soliton gas, in the framework of the focusing nonlinear Schr\"{o}dinger equation (NLSE). Multi-soliton solutions containing of up to 128 solitons are generated with the new recursive numerical scheme for the Zakharov-Mikhailov variant of the dressing method. The statistical properties of these solutions are then studied with repsect to soliton gas density and velocity distribution. Special attention is paid to examination of rogue waves that appear as a result of multi-soliton collisions.