We report the spectroscopic confirmation and modeling of the quadruply imaged quasar GRAL 113100–441959, the first gravitational lens (GL) to be discovered from a machine learning technique that only relies on the relative positions and fluxes of the observed images without considering colour informations. Follow-up spectra obtained with Keck/LRIS reveal the lensing nature of this quadruply imaged quasar with redshift z_s = 1.090 ± 0.002, but show no evidence of the central lens galaxy. Using the image positions and G-band flux ratios provided by Gaia Data Release 2 as constraints, we modeled the system with a singular power-law elliptical mass distribution (SPEMD) plus external shear, to different levels of complexity. We show that relaxing the isothermal constraint of the SPEMD does not lead to statistically significant different results in terms of fitting the lensing data. We thus simplified the SPEMD to a singular isothermal ellipsoid to estimate the Einstein radius of the main lens galaxy θ_E = 0.″851, the intensity and position angle of the external shear (γ,θ_γ) = (0.044, 11.°5), and we predict the lensing galaxy position to be (θ_gal,1, θ_gal,2) = (−0.″424, −0.″744) with respect to image A. We provide time delay predictions for pairs of images, assuming a plausible range of lens redshift values z_l between 0.5 and 0.9. Finally, we examine the impact on time delays of the so-called source position transformation, a family of degeneracies existing between different mass density profiles that reproduce most of the lensing observables equally well. We show that this effect contributes significantly to the time delay error budget and cannot be ignored during the modeling. This has implications for robust cosmography applications of lensed systems. GRAL 113100–441959 is the first in a series of seven new spectroscopically confirmed GLs discovered from Gaia Data Release 2.