code: def heterogzygote_counts(paired):
      
      work_dir = utils.safe_makedir(os.path.join(dd.get_work_dir(paired.tumor_data), "structural", "counts"))
      key = "germline_het_pon"
      het_bed = tz.get_in(["genome_resources", "variation", key], paired.tumor_data)
      vr = bedutils.population_variant_regions([x for x in [paired.tumor_data, paired.normal_data] if x])
      cur_het_bed = bedutils.intersect_two(het_bed, vr, work_dir, paired.tumor_data)
      tumor_counts = _run_collect_allelic_counts(cur_het_bed, key, work_dir, paired.tumor_data)
      normal_counts = (_run_collect_allelic_counts(cur_het_bed, key, work_dir, paired.normal_data)
                       if paired.normal_data else None)
      if normal_counts:
          tumor_counts, normal_counts = _filter_by_normal(tumor_counts, normal_counts, paired.tumor_data)
      return tumor_counts, normal_counts
docs: Provide tumor/normal counts at population heterozyogte sites with CollectAllelicCounts.