In this paper we develop a potential theory for strongly degenerate parabolic operators of the formL:=∇X⋅(A(X,Y,t)∇X)+X⋅∇Y−∂t, in unbounded domains of the formΩ={(X,Y,t)=(x,xm,y,ym,t)∈Rm−1×R×Rm−1×R×R|xm>ψ(x,y,ym,t)}, where ψ is assumed to satisfy a uniform Lipschitz condition adapted to the dilation structure and the (non-Euclidean) Lie group underlying the operator L. Concerning A=A(X,Y,t) we assume that A is bounded, measurable, symmetric and uniformly elliptic (as a matrix in Rm). Beyond the solvability of the Dirichlet problem and other fundamental properties our results include scale and translation invariant boundary comparison principles, boundary Harnack inequalities and doubling properties of associated parabolic measures. All of our estimates are translation- and scale-invariant with constants only depending on the constants defining the boundedness and ellipticity of A and the Lipschitz constant of ψ. Our results represent a version, for operators of Kolmogorov type with bounded, measurable coefficients, of the by now classical results of Fabes and Safonov, and several others, concerning boundary estimates for uniformly parabolic equations in (time-dependent) Lipschitz type domains.