  The bispectrum of single-field inflationary trajectories in which the speed
of sound of the inflationary trajectories $c_s$ is constant but not equal to
the speed of light $c=1$ is explored. The trajectories are generated as random
realisations of the Hubble Slow-Roll (HSR) hierarchy and the bispectra are
calculated using numerical techniques that extends previous work. This method
allows for out-of-slow-roll models with non-trivial time dependence and
arbitrarily low $c_s$. The ensembles obtained using this method yield
distributions for the shape and scale-dependence of the bispectrum and their
relations with the standard inflationary parameters such as scalar spectral
tilt $n_s$ and tensor-to-scalar ratio $r$. The distributions demonstrate the
squeezed-limit consistency relations for arbitrary single-field inflationary
models.
