Note that this has the same functional form as the position-space probability distribution. This is specific to the problem of the simple harmonic oscillator and arises due to the symmetry between x and p in the equations of motion.