We include finite top- and bottom-mass effects in the next-to-next-to-leading order parton shower (NNLOPS) event generator for inclusive Higgs boson production in gluon fusion based upon the POWHEG+MiNLO approach. Since fixed-order results for quark-mass effects only reach NLO accuracy, we add them to the NNLOPS generator at that accuracy. We explore uncertainties related to the unknown all-order logarithmic structure of bottom-mass effects by comparing the assumption of full exponentiation to no exponentiation at all. Phenomenological results showing the effects of finite quark-masses in the NNLOPS simulation are presented. These suggest that the aforementioned uncertainty is well contained within the envelope of plain renormalization and factorization scale uncertainties.