Let Ω ⊂ R n be a bounded NTA-domain and let Ω T = Ω × (0, T ) for some T > 0. We study the boundary behaviour of non-negative solutions to the equationWe assume that A(x, t) = {a ij (x, t)} is measurable, real, symmetric and thatfor some constant β ≥ 1 and for some non-negative and real-valued function λ = λ(x) belonging to the Muckenhoupt class A 1+2/n (R n ). Our main results include the doubling property of the associated parabolic measure and the Hölder continuity up to the boundary of quotients of non-negative solutions which vanish continuously on a portion of the boundary. Our results generalize previous results of Fabes, Kenig, Jerison, Serapioni, see [18], [19], [20], to a parabolic setting.