We report a detailed analysis of the energy spectra, second- and high-order structure functions of velocity differences in superfluid counterflow turbulence, measured in a wide range of temperatures and heat fluxes. We show that the one-dimensional energy spectrum (averaged over the plane, parallel to the channel wall), directly measured as a function of the wall-normal wave vector , gives more detailed information on the energy distribution over scales than the corresponding second-order structure function . In particular, we discover two intervals of with different apparent exponents: for and for . Here denotes the wave number that separates scales with relatively strong (for ) and relatively weak (for ) coupling between the normal-fluid and superfluid velocity components. We interpret these ranges as cascade-dominated and mutual-friction-dominated intervals, respectively. The general behavior of the experimental spectra agrees well with the predicted spectra [L'vov and Pomyalov, Phys. Rev. B 97, 214513 (2018)2469-995010.1103/PhysRevB.97.214513]. Analysis of the -order structure functions statistics shows that in the energy-containing interval, the statistics of counterflow turbulence is close to Gaussian, similar to the classical hydrodynamic turbulence. In the cascade- and mutual-friction-dominated intervals, we found some modest enhancement of intermittency with respect to its level in classical turbulence. However, at small scales, the intermittency becomes much stronger than in the classical turbulence.