code: def statstable(args):
      
      p = OptionParser(statstable.__doc__)
      opts, args = p.parse_args(args)
  
      if len(args) < 1:
          sys.exit(not p.print_help())
  
      gff_files = args
      for metric in metrics:
          logging.debug("Parsing files in `{0}`..".format(metric))
  
          table = {}
          for x in gff_files:
              pf = op.basename(x).split(".")[0]
              numberfile = op.join(metric, pf + ".txt")
              ar = [int(x.strip()) for x in open(numberfile)]
              sum = SummaryStats(ar).todict().items()
              keys, vals = zip(*sum)
              keys = [(pf, x) for x in keys]
              table.update(dict(zip(keys, vals)))
  
          print(tabulate(table), file=sys.stderr)
docs: %prog statstable *.gff
  
      Print gene statistics table.