We analyze the spatial and energy dependence of the local density of states in a SNS junction. We model our system as a one-dimensional tight-binding chain which we solve exactly by numerical diagonalization. We calculate the dependence of the Andreev bound states on position, phase difference, gate voltage, and coupling with the superconducting leads. Our results confirm the physics predicted by certain analytical approximations, but reveal a much richer set of phenomena beyond the grasp of these approximations, such as the metamorphosis of the discrete states of the normal link (the normal bound states) into Andreev bound states as the leads become superconducting.