Best-fitting models resulting from the MCMC analysis as described in Sect. 4, using the light curves detrended by the classical (solid) and time-shift (dashed) procedures. For the latter only, we display the 1σ uncertainties (gray shaded area). These uncertainties have been evaluated using 10⁴ transit models generated using transit and phase modulation parameters randomly taken from the MCMC posteriors.