A crossover from a non-Gaussian to Gaussian subdiffusion has been observed ubiquitously in various polymeric and molecular glassformers. We have developed a framework that generalizes the fractional Brownian motion model to incorporate non-Gaussian features by introducing a jump kernel. We illustrate that the non-Gaussian fractional Brownian motion model accurately characterizes the subdiffusion crossover. From the solutions of the non-Gaussian fractional Brownian motion model, we gain insights into the nature of van Hove self-correlation in non-Gaussian subdiffusive regime, which is found to exhibit exponential tails, providing first such experimental evidence in molecular glassformers. The validity of the model is supported by comparison with incoherent quasielastic neutron scattering data obtained from several molecular and polymeric glassformers.