I have re-stacked with more stringent quality criteria (FWHM and roundness). The noise is higher, but the stars look more round. Problem sorted I think! I’ve looked back at my determinations of transient locations, and the difference from the survey telescope position in this case isn’t unusual. Experimenting with the fit order doesn’t really change the measured position either.