We present, how a systematic study of $$B \rightarrow D\ell N$$B→DℓN ($$D \rightarrow K \ell N$$D→KℓN) decays with $$\ell =\mu ,\tau $$ℓ=μ,τ, at Belle II (BESIII) can provide unambiguous signature of a heavy neutrino N and/or constrain its mixing with active neutrinos $$\nu _\ell $$νℓ, which is parameterized by $$| U_{\ell N} |^2$$|UℓN|2. Our constraint on $$\vert U_{\mu N} \vert ^2$$|UμN|2 that can be achieved from the full Belle II data is comparable with what can be obtained from the much larger data set of the upgraded LHCb. Additionally, our method offers better constraint on $$\vert U_{\mu N} \vert ^2$$|UμN|2 for mass of sterile neutrino $$m_N < 2$$mN<2 GeV. We can also probe the Dirac and Majorana nature of N by observing the sequential decay of N, including suppression from observation of a displaced vertex as well as helicity flip, for Majorana N.