code: def coge(args):
      
      p = OptionParser(coge.__doc__)
      opts, args = p.parse_args(args)
  
      if len(args) < 1:
          sys.exit(not p.print_help())
  
      gffs = args
      for gff in gffs:
          atoms = op.basename(gff).split(".")
          gid = atoms[-2]
          assert gid.startswith("gid")
          gid = get_number(gid)
          genomefasta = "genome_{0}.faa.fasta".format(gid)
          species = "_".join(atoms[0].split("_")[:2])
          cdsfasta = species + ".cds.fasta"
          load([gff, genomefasta, "--id_attribute=Parent",
                "--outfile={0}".format(cdsfasta)])
docs: %prog coge *.gff
  
      Prepare coge datasets.