We consider the mixed Dirichlet-conormal problem for the heat equation on cylindrical domains with a bounded and Lipschitz base 12 subset of Rd and a time-dependent separation Lambda. Under certain mild regularity assumptions on Lambda, we show that for any q > 1 sufficiently close to 1, the mixed problem in Lq is solvable. In other words, for any given Dirichlet data in the parabolic Riesz potential space L1q and the Neumann data in Lq, there is a unique solution and the non-tangential maximal function of its gradient is in Lq on the lateral boundary of the domain. When q = 1, a similar result is shown when the data is in the Hardy space. Under the additional condition that the boundary of the domain 12 is Reifenberg-flat and the separation is locally sufficiently close to a Lipschitz function of m variables, where m = 0, ... , d - 2, with respect to the Hausdorff distance, we also prove the unique solvability result for any q is an element of (1, (m + 2)/(m + 1)). In particular, when m = 0, i.e., Lambda is Reifenberg-flat of co-dimension 2, we derive the Lq solvability in the optimal range q is an element of (1, 2). For the Laplace equation, such results were established in [25,24,5] and [13]. (c) 2022 Elsevier Inc. All rights reserved.