A Monte Carlo study of the mixed phase region in the 3D Ashkin–Teller model