code: def gene_ids(self, contig=None, strand=None):
          
          return self._all_feature_values(
              column="gene_id",
              feature="gene",
              contig=contig,
              strand=strand)
docs: What are all the gene IDs
          (optionally restrict to a given chromosome/contig and/or strand)