We study the eigenvalue problem for a self-adjoint 1D Dirac operator. It is known that, near an energy level where the square of the potential makes a simple well, the eigenvalues are approximated by a Bohr–Sommerfeld type quantization rule. A remarkable difference from the Schrödinger case appears in the Maslov correction term. This fact was recently found by Hirota [Real eigenvalues of a non-self-adjoint perturbation of the self-adjoint Zakharov–Shabat operator, J. Math. Phys. 58 (2017) 102–108] under the analyticity condition of the potential using a complex WKB method. In this paper, we approach this problem with a microlocal technique focusing on the asymptotic behavior of the eigenfunction along the characteristic set to generalize the result to [Formula: see text] potentials.