Residue renumbering fails when the gap tracker > len(gap_letter).
We can extend gap_letter, to include two character identifiers. However, molecular visualization tools (PyMOL and Chimera X), do not recognize them, and some residues may be merged.
Proposed idea: include more identifiers in gap_letter, but print a warning if gap_tracker > A-Z+a-z+0-9.
Add one warning:
WARNING: Chain xxx contains a lot of gaps. Check alignment!
And an error message if the gap_tracker > len(gap_letter).
ERROR: Too many gaps, check your alignment!
Residue renumbering fails when the gap tracker > len(gap_letter).
We can extend gap_letter, to include two character identifiers. However, molecular visualization tools (PyMOL and Chimera X), do not recognize them, and some residues may be merged.
Proposed idea: include more identifiers in gap_letter, but print a warning if gap_tracker > A-Z+a-z+0-9.
Add one warning:
And an error message if the gap_tracker > len(gap_letter).