Low-energy electron holography (LEEH) is one of the few techniques capable of imaging large and complex three-dimensional molecules, such as proteins, on the single molecule level at subnanometer resolution. During the imaging process, the structural information about the object is recorded both in the amplitude and in the phase of the hologram. In low energy electron holography imaging of proteins, the object's amplitude distribution, which directly reveals molecular size and shape on the single-molecule level, can be retrieved via a onestep reconstruction process. However, such a one-step reconstruction routine cannot directly recover the phase information encoded in the hologram. In order to extract the full information about the imaged molecules, we thus implemented an iterative phase retrieval algorithm and applied it to experimentally acquired low-energy electron holograms, reconstructing the phase shift induced by the protein along with the amplitude data. We show that phase imaging can map the projected atomic density of the molecule given by the number of atoms in the electron path. This directly implies a correlation between reconstructed phase shift and projected mean inner potential of the molecule, and thus a sensitivity to local changes in potential, an interpretation that is further substantiated by the strong phase signatures induced by localized charges.