Multijet production rates in neutral current deep inelastic scattering (DIS) have been measured in the range of exchanged boson virtualities 10<Q2<5000 GeV2. The data were taken at the ep collider HERA with centre-of-mass energy [Formula: see text] using the ZEUS detector and correspond to an integrated luminosity of 82.2 pb-1. Jets were identified in the Breit frame using the kTcluster algorithm in the longitudinally invariant inclusive mode. Measurements of differential multijet cross sections are presented as functions of jet transverse energy [Formula: see text], pseudorapidity [Formula: see text] and Q2with [Formula: see text] and [Formula: see text]. Next-to-leading-order QCD calculations describe the data well. The value of the strong coupling constant αs(MZ), determined from the ratio of the trijet to dijet cross sections, is [Formula: see text].