code: def gcn(args):
      
      p = OptionParser(gcn.__doc__)
      p.set_cpus()
      p.set_tmpdir(tmpdir="tmp")
      p.set_outfile()
      opts, args = p.parse_args(args)
  
      if len(args) < 2:
          sys.exit(not p.print_help())
  
      exonbed = args[0]
      canvasvcfs = args[1:]
      tsvfile = opts.outfile
      tmpdir = opts.tmpdir
  
      mkdir(tmpdir)
      set_tempdir(tmpdir)
  
      df = vcf_to_df(canvasvcfs, exonbed, opts.cpus)
      for suffix in (".avgcn", ".medcn"):
          df_to_tsv(df, tsvfile, suffix)
docs: %prog gcn gencode.v26.exonunion.bed data/*.vcf.gz
  
      Compile gene copy njumber based on CANVAS results.