code: def _run_genotype_gvcfs_genomicsdb(genomics_db, region, out_file, data):
      
      if not utils.file_exists(out_file):
          with file_transaction(data, out_file) as tx_out_file:
              broad_runner = broad.runner_from_config(data["config"])
              params = ["-T", "GenotypeGVCFs",
                        "--variant", "gendb://%s" % genomics_db,
                        "-R", dd.get_ref_file(data),
                        "--output", tx_out_file,
                        "-L", bamprep.region_to_gatk(region)]
              params += ["-ploidy", str(ploidy.get_ploidy([data], region))]
              # Avoid slow genotyping runtimes with improved quality score calculation in GATK4
              # https://gatkforums.broadinstitute.org/gatk/discussion/11471/performance-troubleshooting-tips-for-genotypegvcfs/p1
              params += ["--use-new-qual-calculator"]
              resources = config_utils.get_resources("gatk", data["config"])
              params += [str(x) for x in resources.get("options", [])]
              cores = dd.get_cores(data)
              memscale = {"magnitude": 0.9 * cores, "direction": "increase"} if cores > 1 else None
              broad_runner.run_gatk(params, memscale=memscale)
      return vcfutils.bgzip_and_index(out_file, data["config"])
docs: GenotypeGVCFs from a merged GenomicsDB input: GATK4.
              ropts += [str(x) for x in resources.get("options", [])]
  
      No core scaling -- not yet supported in GATK4.