We study quantization of a minimally gauged massless Rarita-Schwinger field, by both Dirac bracket and functional integral methods. The Dirac bracket approach in covariant radiation gauge leads to an anticommutator that has a non-singular limit as gauge fields approach zero, is manifestly positive semidefinite, and is Lorentz invariant. The constraints also have the form needed to apply the Faddeev-Popov method for deriving a functional integral, using the same constrained Hamiltonian and inverse constraint matrix that appear in the Dirac bracket approach.