The Kitaev honeycomb lattice system is an important model of topological materials whose phase diagram exhibits both abelian and non-abelian topological phases. The latter, a so-called Ising phase, is related to topological superconductors. Its quasiparticle excitations, which are formed by Majorana fermions attached to vortices, show non-abelian fractional statistics and are known as Ising anyons. We investigate dislocation defects in the Ising phase of the Kitaev honeycomb model. After introducing them to the system, we accordingly generalize our solution of this model to the situation with the defects. The important part of this effort is developing an appropriate Jordan-Wigner fermionization procedure. It is expected that the presence of defects manifests itself by the formation of fermionic zero-energy modes around the defect end points. We numerically confirm this expectation and further investigate properties of these modes. The computational potential of our technique is demonstrated for both diagonalization and dynamical simulations. The latter focuses on the process of fusion of the vortex zero-energy modes with the Majorana fermions attached to the defect. This process simulates fusion of non-abelian Ising anyons.