We establish two-sided heat kernel estimates for random conductance models with non-uniformly elliptic (possibly degenerate) stable-like jumps on graphs. These are long range counterparts of well known two-sided Gaussian heat kernel estimates by M.T. Barlow for nearest neighbor (short range) random walks on the supercritical percolation cluster. Unlike the cases for nearest neighbor conductance models, the idea through parabolic Harnack inequalities does not work, since even elliptic Harnack inequalities do not hold in the present setting. As an application, we establish the local limit theorem for the models.