We design and analyze a hybrid high-order method on unfitted meshes to approximate the Stokes interface problem. The interface can cut through the mesh cells in a very general fashion. A cell-agglomeration procedure prevents the appearance of small cut cells. Our main results are inf-sup stability and a priori error estimates with optimal convergence rates in the energy norm. Numerical simulations corroborate these results
