PBH in single field inflation: the effect of shape dispersion and non-Gaussianities