Incompatible numerical manifold method (INMM) uses interpolation functions based on the concept of partition of unity, and considers the asymptotic solution and the discontinuity of displacement. This paper describes the application of INMM to bi-material interfacial crack. The two dimensional near-tip asymptotic displacement functions are added to the trial function approximation. This enables the domain to be modeled by manifold elements without explicitly meshing the crack surfaces. The crack-tip enrichment functions are chosen as those that span the asymptotic displacement fields for an interfacial crack. The INMM facilitates the incorporation of the oscillatory nature of the singularity within a conforming manifold element approximation. The complex stress intensity factors for bi-material interfacial cracks are numerically evaluated. Good agreement between the numerical results and the analytical solutions for benchmark interfacial crack problems is realized.