Context. Flat spectrum radio quasars (FSRQs) can suffer strong absorption above E = 25∕(1 + z) GeV, due to gamma–gamma interaction if the emitting region is at sub-parsec scale from the super-massive black hole (SMBH). Aims. Gamma-ray flares from these astrophysical sources can be used to investigate the location of the high-energy emission region and the physics of the radiating processes. Methods. We present an episode of remarkable gamma-ray flaring activity from FSRQ PKS 2023-07 during April 2016, as detected by both the AGILE and Fermi satellites. An intensive multiwavelength campaign, triggered by Swift, covered the entire duration of the flaring activity, including the peak gamma-ray activity. Results. We report the results of multiwavelength observations of the blazar. We found that during the peak emission, the most energetic photon had an energy of 44 GeV, putting strong constraints on the opacity of the gamma-ray dissipation region. The overall spectral energy distribution (SED) is interpreted in terms of leptonic models for blazar jets, with the emission site located beyond the broad line region (BLR).