We present charmonium and bottomonium correlators and corresponding
reconstructed spectral functions from full QCD calculations in the pseudoscalar
channel. Correlators are obtained using a mixed-action approach, clover-improved
Wilson valence quarks on gauge field configurations generated with $N_f=2+1$ HISQ sea
quarks, with physical strange quark masses and light quark masses corresponding to $m_\pi=315$ MeV.
The charm and bottom quark masses are tuned to reproduce
the experimental mass spectrum of the spin averaged quarkonium vector mesons
from the particle data group. For the spectral
reconstruction, we use models based on perturbative spectral functions from different
frequency regions like resummed thermal contributions around the threshold from pNRQCD and vacuum contributions well above the threshold. We show preliminary results of the reconstructed spectral
function obtained for the first time in our study for full QCD.