docs: def run(self):
          """Run formula balance command"""
  
          # Create a set of excluded reactions
          exclude = set(self._args.exclude)
          count = 0
          unbalanced = 0
          unchecked = 0
          for reaction, result in formula_balance(self._model):
              count += 1
  
              if reaction.id in exclude or reaction.equation is None:
                  continue
  
              if result is None:
                  unchecked += 1
                  continue
  
              left_form, right_form = result
              if right_form != left_form:
                  unbalanced += 1
                  right_missing, left_missing = Formula.balance(
                      right_form, left_form)
  
                  print('{}\t{}\t{}\t{}\t{}'.format(
                      reaction.id, left_form, right_form,
                      left_missing, right_missing))
  
          logger.info('Unbalanced reactions: {}/{}'.format(unbalanced, count))
          logger.info('Unchecked reactions due to missing formula: {}/{}'.format(
              unchecked, count))
          logger.info('Reactions excluded from check: {}/{}'.format(
              len(exclude), count))