Photons propagating in an external magnetic field may oscillate into axions or axion-like particles (ALPs). Such oscillations will lead to characteristic features in the energy spectrum of high-energy photons from astrophysical sources that can be used to probe the existence of ALPs. In this work, we revisit the signatures of these oscillations and stress the importance of a proper treatment of turbulent magnetic fields. We implement axions into ELMAG, a standard tool for modelling in a Monte Carlo framework the propagation of gamma-rays in the Universe, complementing thereby the usual description of photon-axion oscillations with a Monte Carlo treatment of high-energy photon propagation and interactions. We also propose an alternative method of detecting axions through the discrete power spectrum using as observable the energy dependence of wiggles in the photon spectra.