Analytical results on generating phase-equivalent potentials by supersymmetry: Removal and addition of bound states