p-values obtained after performing the likelihood analysis on our datasets. The logarithm of the p-value corresponding to the maximum likelihood obtained for the SBG catalogue (P_SBG) catalogue is plotted against that obtained for 2MRS (P_2MRS) for each dataset. The astrophysical model considered to build our datasets corresponds to model A, the baseline catalogue and a 1 nG EGMF. The values reported by Auger in Caccianiga et al. (2019) are shown with large black full circle. The individual distributions of the different quantities plotted are shown on top of the coordinate axis. The left panel shows the case of the Sun+Planck GMF model for two different values of the coherence length. The right panel shows the case of the JF12+Planck model. In this case we also show datasets produced after excluding galaxies from the Virgo association from the source catalogue (see legend).