Best-fit models to the HST data from our atmosphere retrievals (top and upper-middle) and stellar contamination models (lower-middle). For the stellar contamination models, the solid lines depict the transit light source effect for the maximum spot filling factor, as defined in Rackham et al. (2018), while the dashed lines represent the mean. Transit depths from Diamond-Lowe et al. (2020) are also shown (bottom) and suggest the presence of solar spots with faculae.