docs: def run_makeblastdb(infile, dbtype, outdir=''):
      """Make the BLAST database for a genome file.
  
      Args:
          infile (str): path to genome FASTA file
          dbtype (str): "nucl" or "prot" - what format your genome files are in
          outdir (str): path to directory to output database files (default is original folder)
  
      Returns:
          Paths to BLAST databases.
  
      """
      # TODO: add force_rerun option
      # TODO: rewrite using utils function command
  
      # Output location
      og_dir, name, ext = utils.split_folder_and_path(infile)
      if not outdir:
          outdir = og_dir
      outfile_basename = op.join(outdir, name)
  
      # Check if BLAST DB was already made
      if dbtype == 'nucl':
          outext = ['.nhr', '.nin', '.nsq']
      elif dbtype == 'prot':
          outext = ['.phr', '.pin', '.psq']
      else:
          raise ValueError('dbtype must be "nucl" or "prot"')
      outfile_all = [outfile_basename + x for x in outext]
      db_made = True
      for f in outfile_all:
          if not op.exists(f):
              db_made = False
  
      # Run makeblastdb if DB does not exist
      if db_made:
          log.debug('BLAST database already exists at {}'.format(outfile_basename))
          return outfile_all
      else:
          retval = subprocess.call('makeblastdb -in {} -dbtype {} -out {}'.format(infile, dbtype, outfile_basename), shell=True)
          if retval == 0:
              log.debug('Made BLAST database at {}'.format(outfile_basename))
              return outfile_all
          else:
              log.error('Error running makeblastdb, exit code {}'.format(retval))