The EXOD search for faint transients in XMM-Newton observations: Method and discovery of four extragalactic Type I X-ray bursters