Residuals obtained from the fit of the timing model for PSR J1012−4235. Top plot: Post-fit residuals showing ToAs from Parkes, GBT, MeerKAT, and Fermi-LAT data, with error bars representing 1σ uncertainty in ToAs. Bottom plot: Residuals of MeerKAT ToAs, excluding the fit for Shapiro delay effect but fixing all the other system and orbital parameters obtained from the full fit. The orbital phase is measured from periastron and the peak in the Shapiro delay signature occurs around the phase of superior conjunction of the pulsar.