Statistics of GMRES with embedded ensemble propagation for the efficient solution of parametric linear systems in uncertainty quantification of computational models