It is widely accepted that the Feynman integral is one of the most promising methodologies for defining a generally covariant formulation of nonperturbative interacting quantum field theories (QFTs) without a fixed prearranged causal background. Recent literature suggests that if the spacetime metric is not fixed, e.g. because it is to be quantized along with the other fields, one may not be able to avoid considering the Feynman integral in the original Lorentz signature, without Wick rotation. Several mathematical phenomena are known, however, which are at some point showstoppers to a mathematically sound definition of Feynman integral in Lorentz signature. The Feynman integral formulation, however, is known to have a differential reformulation, called to be the master Dyson--Schwinger (MDS) equation for the field correlators. In this paper it is shown that a particular presentation of the MDS equation can be cast into a mathematically rigorously defined form: the involved function spaces and operators can be strictly defined and their properties can be established. Therefore, MDS equation can serve as a substitute for the Feynman integral, in a mathematically sound formulation of constructive QFT, in arbitrary signature, without a fixed background causal structure. It is also shown that even in such a generally covariant setting, there is a canonical way to define the Wilsonian regularization of the MDS equation. The main result of the paper is a necessary and sufficient condition for the regularized MDS solution space to be nonempty, for conformally invariant Lagrangians. This theorem also provides an iterative approximation algorithm for obtaining regularized MDS solutions, and is guaranteed to be convergent whenever the solution space is nonempty. The algorithm could eventually serve as a method for putting Lorentz signature QFTs onto lattice, in the original metric signature.