We find the asymptotic representation of the solution of the variable-order fractional diffusion equation, which remains unsolved since it was proposed by Chechkin, Gorenflo, and Sokolov [J. Phys. A, 38, L679 (2005)JPHAC50305-447010.1088/0305-4470/38/42/L03]. We identify a new advection term that causes ultraslow spatial aggregation of subdiffusive particles due to dominance over the standard advection and diffusion terms in the long-time limit. This uncovers the anomalous mechanism by which nonuniform distributions can occur. We perform Monte Carlo simulations of the underlying anomalous random walk and find excellent agreement with the asymptotic solution.