code: def _call_variants(example_dir, region_bed, data, out_file):
      
      tf_out_file = "%s-tfrecord.gz" % utils.splitext_plus(out_file)[0]
      if not utils.file_exists(tf_out_file):
          with file_transaction(data, tf_out_file) as tx_out_file:
              model = "wes" if strelka2.coverage_interval_from_bed(region_bed) == "targeted" else "wgs"
              cmd = ["dv_call_variants.py", "--cores", dd.get_num_cores(data),
                     "--outfile", tx_out_file, "--examples", example_dir,
                     "--sample", dd.get_sample_name(data), "--model", model]
              do.run(cmd, "DeepVariant call_variants %s" % dd.get_sample_name(data))
      return tf_out_file
docs: Call variants from prepared pileup examples, creating tensorflow record file.