In this paper we consider the multispecies stirring process on the discrete torus. We prove a large deviation principle for the trajectory of the vector of densities of the different species. The technique of proof consists in extending the method of the foundational paper [15] based on the superexponential estimate to the multispecies setting. This requires a careful choice of the corresponding weakly asymmetric dynamics, which is parametrized by fields depending on the various species. We also prove the hydrodynamic limit of this weakly asymmetric dynamics, which is similar to the ABC model in [12, 2]. Using the appropriate asymmetric dynamics, we also obtain that the mobility matrix relating the drift currents to the fields coincides with the covariance matrix of the reversible multinomial distribution, which then further leads to the Einstein relation.