Data release for an extensive analysis of the GW190521 gravitational wave event, which is consistent with a high mass binary black hole system with a very short signal, detected during the O3a observation run of the Advanced LIGO and Virgo detectors. In this analysis we use the current generation of phenomenological waveform model families, IMRPhenomX and IMRPhenomT, as implemented in the LALSuite code (https://doi.org/10.7935/GT1W-FZ16). This data release corresponds to the paper "A detailed analysis of GW190521 with phenomenological waveform models" (https://arxiv.org/abs/2105.06360), which should be cited when using these data. It contains the posterior samples of the default runs using IMRPhenomXPHM and IMRPhenomTPHM, carried out using the Markov Chain Monte Carlo (MCMC) sampling implementation of the LALInference code and their corresponding command line.