Disentangling spectra of young stars: Unveiling the accretion spectrum with Gaussian Processes