In this work, we implement Gaussian process regression to reconstruct the expansion history of the universe in a model-agnostic manner, using the Pantheon-Plus SN-Ia compilation in combination with two different BAO measurements (SDSS-IV and DESI DR1). In both the reconstructions, the ΛCDM model is always included in the 95\% confidence intervals. We find evidence that the DESI LRG data at zeff=0.51 is not an outlier within our model-independent framework. We study the Om-diagnostics and the evolution of the total equation of state (EoS) of our universe, which hint towards the possibility of a quintessence-like dark energy scenario with a very slowly varying EoS, and a phantom-crossing in higher z. The entire exercise is later complemented by considering two more SN-Ia compilations - DES-5YR and Union3 - in combination with DESI BAO. Reconstruction with the DESI BAO + DES-5YR SN data sets predicts that the ΛCDM model lies outside the 3σ confidence levels, whereas with DESI BAO + Union3 data, the ΛCDM model is always included within 1σ. We also report constraints on H0rd from our model-agnostic analysis, independent of the pre-recombination physics. Our results point towards an ≈ 2σ discrepancy between the DESI + Pantheon-Plus and DESI + DES-5YR data sets, which calls for further investigation.