docs: def emboss_pepstats_on_fasta(infile, outfile='', outdir='', outext='.pepstats', force_rerun=False):
      """Run EMBOSS pepstats on a FASTA file.
  
      Args:
          infile: Path to FASTA file
          outfile: Name of output file without extension
          outdir: Path to output directory
          outext: Extension of results file, default is ".pepstats"
          force_rerun: Flag to rerun pepstats
  
      Returns:
          str: Path to output file.
  
      """
  
      # Create the output file name
      outfile = ssbio.utils.outfile_maker(inname=infile, outname=outfile, outdir=outdir, outext=outext)
  
      # Run pepstats
      program = 'pepstats'
      pepstats_args = '-sequence="{}" -outfile="{}"'.format(infile, outfile)
      cmd_string = '{} {}'.format(program, pepstats_args)
      ssbio.utils.command_runner(cmd_string, force_rerun_flag=force_rerun, outfile_checker=outfile, silent=True)
  
      return outfile