In the quantum theory of electromagnetic fields Maxwell’s equations play two roles.

One is to describe the behaviour of the actual field observables which measure the combined effects of all possible numbers of photons, and the other is that they are satisfied by something that is as close as possible to being the “wave function” for a single photon.

I say “as close as possible” and put “wave function” in scare quotes because it does not satisfy all the properties of a non-relativistic wave function. In relativistic theories, the concept of strict localization does not exist. It can only be approximated for massive particles in frames where they have low momentum, and cannot be done at all for photons. But nonetheless, (as discussed in this survey article (.pdf) by Iwo Bialynicki-Birula) with appropriate normalization, a function satisfying the complex form of Maxwell’s equations can be used to generate probabilities for detection of a single photon in various experimental contexts.

See also the answer by ‘Chiral Anomaly’ to this question at physics stack-exchange.