Solution space heterogeneity of the random $K$-satisfiability problem: Theory and simulations