code: def align(args):
      
      p = OptionParser(align.__doc__)
      p.set_cpus()
      opts, args = p.parse_args(args)
  
      if len(args) != 1:
          sys.exit(not p.print_help())
  
      clustfile, = args
      parallel_musclewrap(clustfile, opts.cpus)
docs: %prog align clustfile
  
      Align clustfile to clustSfile. Useful for benchmarking aligners.