We exhibit the supersymmetric quantum mechanical structure of the full 3+1 dimensional Dirac equation considering ‘mass’ as a function of coordinates. Its usefulness in solving potential problems is discussed with specific examples. We also discuss the ‘physical’ significance of the supersymmetric states in this formalism.