We study gravitational lensing in strong-field limit by a static spherically symmetric black hole in quartic scalar field Horndeski gravity having additional hair parameter q, evading the no-hair theorem. We find an increase in the deflection angle alpha _D, photon sphere radius x_{ps}, and angular position theta _{infty } that increases more quickly while angular separation s more slowly, but the ratio of the flux of the first image to all other images r_{mag} decreases rapidly with increasing magnitude of the hair q. We also discuss the astrophysical consequences in the supermassive black holes at the centre of several galaxies and note that the black holes in Horndeski gravity can be quantitatively distinguished from the Schwarzschild black hole. Notably, we find that the deviation Delta theta _{infty } of black holes in Horndeski gravity from their general relativity (GR) counterpart, for supermassive black holes Sgr A* and M87*, for q=-0.2, respectively, can reach as much as 2.4227~mu as and 1.82026~mu as while Delta s is about 0.04650~mu as for Sgr A* and 0.03493~mu as for M87*. The ratio of the flux of the first image to all other images suggest that the Schwarzschild images are brighter than those of the black holes in Horndeski gravity, wherein the deviation |Delta r_{mag}| is as much as 0.70673. The results suggest that observational tests of hairy black holes in Horndeski gravity are indeed feasible. Taking the supermassive black holes Sgr A* and M87* as the lens, we also compare our hairy Horndeski black holes observable signatures with those of the neutral Horndeski black holes, Galileon black holes and charged Horndeski black holes. It turns out that although it is possible to detect some effects of the strong deflection lensing by the hairy Horndeski black holes and other black holes with the Event Horizon Telescope (EHT) observations, but it is unconvincing to discern these black holes as deviations are {mathcal {O}}(mu as). We also find that the shadow size is consistent with EHT observation if the deviation parameter q in (-0.281979,0)