We report on the first master-field simulations of QCD with 2+1 dynamical quark flavours using non-perturbatively improved stabilised Wilson fermions. Our simulations are performed at a lattice spacing of 0.094 fm with 96 and 192 points in each direction. On both lattices, the pion and kaon masses are equal to 270 and 450 MeV, respectively, and $m_{\pi} L$ thus reaches an unprecedented value of $25$ on the larger lattice. This setup matches a single point on a chiral trajectory with fixed trace of the quark mass matrix and allows for comparisons to standard large-scale simulations. We present our algorithmic setup and performance measures, and report about our experience in thermalising large master-field lattices with fermions.