The disordered XXZ model is a prototype model of the many-body localization (MBL)
transition. Despite numerous studies of this model, the available
numerical evidence of multifractality of its eigenstates is not very
conclusive due to severe finite size effects. Moreover it is not clear
if similarly to the case of single-particle physics, multifractal
properties of the many-body eigenstates are related to anomalous transport,
which is observed in this model. In this work, using a state-of-the-art,
massively parallel, numerically exact method, we study systems of
up to 24 spins and show that a large fraction of the delocalized phase
flows towards ergodicity in the thermodynamic limit, while a region
immediately preceding the MBL transition appears to be multifractal
in this limit. We discuss the implication of our finding on the mechanism
of subdiffusive transport.