We describe the details of the calculation of the full colour NNLO QCD corrections to jet production observables at the LHC with antenna subtraction. All relevant matrix elements for the process pp → jj at NNLO in full colour are colour-decomposed and given in a Nc and nf expansion, making identification of leading and subleading colour contributions transparent. The colour-ordered antenna subtraction method has previously successfully been used to construct the NNLO subtraction terms for processes with up to five partons or in the leading colour approximation. However, it is challenged by the more involved subleading colour structure of the squared matrix elements in processes with six or more partons. Here, we describe the methods needed to successfully construct the NNLO subtraction terms for the subleading colour contributions to dijet production within the antenna subtraction formalism.