A multinucleon neutrino interaction mode is implemented into the NEUT neutrino interaction generator. The model allows the simulation of fully differential event kinematics and uses the prediction of Nieves et al. for the total cross section and final-state charged lepton kinematics. The model is fit to T2K near-detector data using a Markov chain monte carlo method, and the goodness-of-fit is evaluated. The CCQE model in NEUT is fit to published data from the MiniBooNE experiment to inform the size of the prior systematic uncertainties associated with the charged-current quasielastic interaction mode. The priors were applied in the 2012 T2K analyses. The fit required development of monte carlo reweighting functions which allowed changes to the nuclear model to be simulated. Both the fits, and the development of the reweighting functions, are discussed in detail.