code: def gnomad_coverage_stats(mt_path):
      mt = hl.read_matrix_table(mt_path)
  
      def get_coverage_expr(mt):
          cov_arrays = hl.literal(
              {
                  x: [1, 1, 1, 1, 1, 1, 1, 1, 0]
                  if x >= 50
                  else [1, 1, 1, 1, 1, 1, 1, 0, 0]
                  if x >= 30
                  else ([1] * (i + 2)) + ([0] * (7 - i))
                  for i, x in enumerate(range(5, 100, 5))
              }
          )
  
          return hl.bind(
              lambda array_expr: hl.struct(
                  **{
                      f"over_{x}": hl.int32(array_expr[i])
                      for i, x in enumerate(
                          [1, 5, 10, 15, 20, 25, 30, 50, 100]
                      )
                  }
              ),
              hl.agg.array_sum(
                  hl.case()
                  .when(mt.x >= 100, [1, 1, 1, 1, 1, 1, 1, 1, 1])
                  .when(mt.x >= 5, cov_arrays[mt.x - (mt.x % 5)])
                  .when(mt.x >= 1, [1, 0, 0, 0, 0, 0, 0, 0, 0])
                  .default([0, 0, 0, 0, 0, 0, 0, 0, 0])
              ),
          )
  
      mt = mt.annotate_rows(
          mean=hl.agg.mean(mt.x),
          median=hl.median(hl.agg.collect(mt.x)),
          **get_coverage_expr(mt),
      )
      mt.rows()._force_count()
  
imports: from os import path
  from tempfile import TemporaryDirectory, NamedTemporaryFile
  import hail
  from .resources import *
  from .utils import benchmark
label_desc: Use image node_id 25 for calling a global function with example usage: gnomad_coverage_stats(mt_path) and returns: hl