This paper focused on the stability analysis of the Dirichlet-to-Neumann (DN) map for the fractional diffusion equation with a reaction coefficient $q$. The main result provided a Hölder-type stability estimate for the map, which was formulated in terms of the Dirichlet eigenvalues and normal derivatives of eigenfunctions of the operator $A_q:=-\mathrm{\Delta }+q$.