Statistics of An $S$-adic characterization of minimal subshifts with first difference of complexity $1 p(n+1)-p(n)\le2$