docs: def _array(fn, cls, genelist, **kwargs):
      """
      Returns a "meta-feature" array, with len(genelist) rows and `bins`
      cols.  Each row contains the number of reads falling in each bin of
      that row's modified feature.
      """
      reader = cls(fn)
      _local_coverage_func = cls.local_coverage
      biglist = []
      if 'bins' in kwargs:
          if isinstance(kwargs['bins'], int):
              kwargs['bins'] = [kwargs['bins']]
  
      for gene in genelist:
          if not isinstance(gene, (list, tuple)):
              gene = [gene]
          coverage_x, coverage_y = _local_coverage_func(
              reader, gene, **kwargs)
          biglist.append(coverage_y)
      return biglist