Knowledge of the neutrino flux produced by the Neutrinos at the Main Injector (NuMI) beamline is essential to the neutrino oscillation and neutrino interaction measurements of the MINERvA, MINOS+, NOvA and MicroBooNE experiments at Fermi National Accelerator Laboratory. We have produced a flux prediction which uses all available and relevant hadron production data, incorporating measurements of particle production off of thin targets as well as measurements of particle yields from a spare NuMI target exposed to a 120 GeV proton beam. The result is the most precise flux prediction achieved for a neutrino beam in the one to tens of GeV energy region. We have also compared the prediction to in situ measurements of the neutrino flux and find good agreement.