We scrutinize the anomalies in diffusion observed in an extended long-range system of classical rotors, the HMF model. Under suitable preparation, the system falls into long-lived quasi-stationary states for which superdiffusion of rotor phases has been reported. In the present work, we investigate diffusive motion by monitoring the evolution of full distributions of unfolded phases. After a transient, numerical histograms can be fitted by the q -Gaussian form P(x) proportional to {1+(q-1)[x/beta]2}(1/(1-q)) , with parameter q increasing with time before reaching a steady value q approximately 32 (squared Lorentzian). From the analysis of the second moment of numerical distributions, we also discuss the relaxation to equilibrium and show that diffusive motion in quasistationary trajectories depends strongly on system size.