docs: def as_bel(self) -> str:
          """Return this protein modification variant as a BEL string."""
          return 'pmod({}{})'.format(
              str(self[IDENTIFIER]),
              ''.join(', {}'.format(self[x]) for x in PMOD_ORDER[2:] if x in self)
          )