Mean maps of the species for which the analysis is presented in this work for the first time (see Table 1). The mean maps were obtained by averaging the normalized maps of different transitions of the same species, and therefore they are in arbitrary units with peak close to unity. In dashed white lines are plotted the contour levels of the continuum map at 150, 60, and 20 times the value of rms = 0.8 mJy beam⁻¹ (from Mininni et al. 2020). The black dashed lines delimit the contour where the intensity of the mean map drops to 50%, 25%, and 5% of the peak value. The 2D Gaussian fit ellipse to the emission is coincident with the inner black-dashed contour (50% level) for all the molecular species. The cyan dashed line delimits the contour where the intensity of the mean map is equal to 10 times the rms of the map. The size of the 1.″2 beam is indicated in the lower left corner of the first panel (top left).