A search for the production of Higgs boson in association with a single top quark is presented, based on data collected in 2016 with the CMS detector at the LHC at a center of mass energy of 13~\mathrm{TeV}, which corresponds to an integrated luminosity of 35.9~\mathrm{fb}^{-1}. The production cross section of this process is highly sensitive to the absolute values of the modifiers of the top quark-Higgs boson coupling, , and the coupling of vector bosons to the Higgs boson, , as well as their relative signs with respect to the standard model. Analyses using multilepton signatures due to , , and decay modes and that with a single lepton for the decay are combined with a reinterpretation of a measurement of to constrain . For the observed data favor positive values of and exclude values of below about .