Spin chains with a symmetry-protected edge zero modes can be seen as prototypical systems for exploring topological signatures in quantum systems. However in an experimental realization of such a system, spurious interactions may cause the edge zero modes to delocalize. To combat this influence beyond simply increasing the bulk gap, it has been proposed to harness disorder which does not drive the system out of a topological phase. Equipped with numerical tools for constructing locally conserved operators that we introduce, we comprehensively explore the interplay of local interactions and disorder on localized edge modes in such systems. Contrary to established heuristic reasoning, we find that disorder has no effect on the edge mode localization length in the non-interacting regime. Moreover, disorder helps localize only a subset of edge modes in the truly interacting regime. We identify one edge mode operator that behaves as if subjected to a non-interacting perturbation, i.e., shows no disorder dependence. This implies that in finite systems, edge mode operators effectively delocalize at distinct interaction strengths despite the presence of disorder. In essence, our findings suggest that the ability to identify and control the best localized edge mode trumps any gains from introducing disorder.