We make a comparative study of the neutrinoless double beta decay constraints on heavy sterile neutrinos versus other direct and indirect constraints from both lepton number conserving and violating processes, as a sensitive probe of the extent of lepton number violation and possible interference effects in the sterile sector. To this effect, we introduce a phenomenological parametrisation of the simplified one-generation seesaw model with one active and two sterile neutrino states in terms of experimentally measurable quantities, such as active-sterile neutrino mixing angles, CP phases, masses and mass splittings. This simple parametrisation enables us to analytically derive a spectrum of possible scenarios between the canonical seesaw with purely Majorana heavy neutrinos and inverse seesaw with pseudo-Dirac ones. We then go on to constrain the simplified parameters of this model from various experiments in energy, intensity and cosmic frontiers. We emphasise that the constraints from lepton number violating processes strongly depend on the mass splitting between the two sterile states and the relative CP phase between them. This is particularly relevant for the neutrinoless double beta decay constraint, which could now get significantly weaker for small mass splitting and opposite CP parities between the sterile states. It is important to keep this in mind while comparing the neutrinoless double beta decay constraint with the direct search limits in the electron flavour from high-energy colliders.