Abstract: Excited states in the mirror nuclei P-31 and S-31 were populated in the 1p and 1n exit channels of the reaction Ne-20 + C-12, at a beam energy of 33 MeV. The Ne-20 beam was delivered for the first time by the Piave-Alpi accelerator of the Laboratori Nazionali di Legnaro. Angular correlations of coincident gamma-rays and Doppler-shift attenuation lifetime measurements were performed using the multi-detector array GASP in conjunction with the EUCLIDES charged particle detector. In the observed B(E1) strengths, the isoscalar component, amounting to 24% of the isovector one, provides strong evidence for breaking of the isospin symmetry in the A = 31 mass region. Self-consistent beyond mean field calculations using Equation of Motion method based on a chiral potential and including two- and three-body forces reproduce well the experimental B(E1) strengths, reinforcing our conclusion. Coherent mixing from higher-lying states involving the Giant Isovector Monopole Resonance accounts well for the effect observed. The breaking of the isospin symmetry originates from the violation of the charge symmetry of the two- and three-body parts of the potential, only related to the Coulomb interaction.