The paper is concerned with the asymptotic expansion of solutions to the following 2 × 2 Dirac type system: (Formula presented.) with a smooth matrix potential Q∈W1n0,1⊗C2×2 and b1 < 0 < b2. If b2 = −b1 = 1, this equation is equivalent to one dimensional Dirac equation. We apply these formulas to get the asymptotic expansion of the characteristic determinant of the boundary value problem associated with the above equation subject to the general two-point boundary conditions. This expansion directly yields new completeness result for the system of root functions of such BVP with nonregular boundary conditions. © The Author(s), under exclusive licence to Springer Nature Switzerland AG 2024.