The marginal can have several peaks. Fix one observed z=xt, then shrink the noise step to see why the local reverse conditional can still become Gaussian-like.
Fixed setup: Z=Y+√h ε, ε∼N(0,1), with a multimodal pY. The exact conditional is p(y|z) ∝ exp(−(z−y)²/(2h)) pY(y); the dashed curve below is only the zeroth-order approximation N(y; z, h).
Drag the red observation line, or use the slider. When the chart is focused, the left and right arrow keys also move it. The first two rows are peak-scaled to reveal their shapes; the two conditional curves in the last row share one density scale.