Due to the limited detector resolution and the challenges involving the tagging of individual quark flavors, it is difficult to measure the Yukawa coupling of bottom, charm and strange quarks from the final state.A precise prediction of the qTβ sprectrum for Higgs production via quark annihilation allows to access the Yukawa coupling for bottom and lighter quarks from the inital state as the qTβ spectra of the quark flavors show a different shape.Especially in the peak region this allows for a possible discrimination of bottom, charm, and strange initiated production once a full treatment of mass effects is included.I will present an N3LLβ²+ approximate N3LO prediction for Higgs production via quark annihilation where I will in particular focus on the differences to the related Drell-Yan process