docs: def print_run_bidirectional_blast(reference, other_genome, dbtype, outdir):
      """Write torque submission files for running bidirectional blast on a server and print execution command.
  
      Args:
          reference (str): Path to "reference" genome, aka your "base strain"
          other_genome (str): Path to other genome which will be BLASTed to the reference
          dbtype (str): "nucl" or "prot" - what format your genome files are in
          outdir (str): Path to folder where Torque scripts should be placed
  
      """
      # TODO: add force_rerun option
  
      if dbtype == 'nucl':
          command = 'blastn'
      elif dbtype == 'prot':
          command = 'blastp'
      else:
          raise ValueError('dbtype must be "nucl" or "prot"')
  
      r_folder, r_name, r_ext = utils.split_folder_and_path(reference)
      g_folder, g_name, g_ext = utils.split_folder_and_path(other_genome)
  
      # Reference vs genome
      r_vs_g_name = r_name + '_vs_' + g_name
      r_vs_g = r_vs_g_name + '_blast.out'
      if op.exists(op.join(outdir, r_vs_g)) and os.stat(op.join(outdir, r_vs_g)).st_size != 0:
          log.debug('{} vs {} BLAST already run'.format(r_name, g_name))
      else:
          cmd = '{} -query {} -db {} -outfmt 6 -out {}'.format(command, reference, g_name, r_vs_g)
          utils.write_torque_script(command=cmd, err=r_vs_g_name, out=r_vs_g_name, name=r_vs_g_name,
                                    outfile=op.join(outdir, r_vs_g_name) + '.sh',
                                    walltime='00:15:00', queue='regular')
  
      # Genome vs reference
      g_vs_r_name = g_name + '_vs_' + r_name
      g_vs_r = g_vs_r_name + '_blast.out'
      if op.exists(op.join(outdir, g_vs_r)) and os.stat(op.join(outdir, g_vs_r)).st_size != 0:
          log.debug('{} vs {} BLAST already run'.format(g_name, r_name))
      else:
          cmd = '{} -query {} -db {} -outfmt 6 -out {}'.format(command, other_genome, r_name, g_vs_r)
          utils.write_torque_script(command=cmd, err=g_vs_r_name, out=g_vs_r_name, name=g_vs_r_name,
                                    outfile=op.join(outdir, g_vs_r_name) + '.sh',
                                    walltime='00:15:00', queue='regular')