Domain switching pathways fundamentally control performance in ferroelectric thin film devices. In epitaxial bismuth ferrite (BiFeO 3 ) films, the domain morphology is known to influence the multiferroic orders. While both striped and mosaic domains have been observed, the origins of the latter have remained unclear. Here, it is shown that domain morphology is defined by the strain profile across the film-substrate interface. In samples with mosaic domains, X-ray diffraction analysis reveals strong strain gradients, while geometric phase analysis using scanning transmission electron microscopy finds that within 5 nm of the film-substrate interface, the out-of-plane strain shows an anomalous dip while the in-plane strain is constant. Conversely, if uniform strain is maintained across the interface with zero strain gradient, striped domains are formed. Critically, an ex situ thermal treatment, which eliminates the interfacial strain gradient, converts the domains from mosaic to striped. The antiferromagnetic state of the BiFeO 3 is also influenced by the domain structure, whereby the mosaic domains disrupt the long-range spin cycloid. This work demonstrates that atomic scale tuning of interfacial strain gradients is a powerful route to manipulate the global multiferroic orders in epitaxial films.