Multijet production rates in neutral current deep inelastic scattering have been measured in the range of exchanged boson virtualities 10 < Q 2 < 5000 GeV 2 . The data were taken at the ep collider HERA with centre-of-mass energy √ s = 318 GeV using the ZEUS detector and correspond to an integrated luminosity of 82.2 pb −1 . Jets were identified in the Breit frame using the k T cluster algorithm in the longitudinally invariant inclusive mode. Measurements of differential dijet and trijet cross sections are presented as functions of jet transverse energy (E jet T,B ), pseudorapidity (η jet LAB ) and Q 2 with E jet T,B > 5 GeV and −1 < η jet LAB < 2.5. Nextto-leading-order QCD calculations describe the data well. The value of the strong coupling constant α s (M Z ), determined from the ratio of the trijet to dijet cross sections, is α s (M Z ) = 0.1179 ± 0.0013 (stat.) +0.0028 −0.0046 (exp.) +0.0064 −0.0046 (th.).