We present the finite temperature spectra of both bottomonium and charmonium, obtained from a consistent lattice QCD based potential picture. Starting point is the complex in-medium potential extracted on full QCD lattices with dynamical u,d and s quarks, generated by the HotQCD collaboration. Using the generalized Gauss law approach, vetted in a previous study on quenched QCD, we fit Re[ V ] with a single temperature dependent parameter m D , the Debye screening mass, and confirm the up to now tentative values of Im[ V ]. The obtained analytic expression for the complex potential allows us to compute quarkonium spectral functions by solving an appropriate Schrödinger equation. These spectra exhibit thermal widths, which are free from the resolution artifacts that plague direct reconstructions from Euclidean correlators using Bayesian methods. In the present adiabatic setting, we find clear evidence for sequential melting and derive melting temperatures for the different bound states. Quarkonium is gradually weakened by both screening (Re[ V ]) and scattering (Im[ V ]) effects that in combination lead to a shift of their in-medium spectral features to smaller frequencies, contrary to the mass gain of elementary particles at finite temperature.