CERN-LHC. The analysis searches for exotic Higgs boson decays to two light pseudoscalars, in the final state of two b quarks and two tau leptons. This signature is sensitive to the Higgs sector of the next to minimal supersymmetric standard model (NMSSM) as well as two Higgs doublet models extended with a scalar singlet (2HDM+S). The analysis is based on data from proton-proton collisions at a center-of-mass energy of 13 TeV collected with the CMS detector at the LHC, in 2016 corresponding to an integrated luminosity of 35.9 fb-1. No signal above the standard model background expectation is observed. The tables give upper limits on the cross section times branching fraction for Higgs bosons with masses between 15 and 60 GeV.