pept_neighbors maps exact peptide sequence to the neighboring amino acids of the last protein in the database which contains that peptide:
This means that termini may be incorrect if the peptide maps to more than one protein. To fix this, the mapping would need to change from
sequence --> termini to
(sequence, protein_id) --> termini.
This "solution" would cause cascading failure in all output generation functions.