cellink.tl.subset_gene

cellink.tl.subset_gene(gdata, gene_id, gene_snp_map='variant_annotation_vep')

Subset the gdata object to only include variants that are annotated to a specific gene.

It is required that vep annotations have already been added to gdata and that ‘gene_id’ is present in gdata.uns[“variant_annotation_vep”]

Parameters:

• gdata (Anndata object containing the genomic data.)

• gene_id (Gene identifier (e.g., Ensembl gene ID) to subset by.)

Returns:

Subsetted Anndata object containing only variants annotated to the specified gene.