ASE: only consider genes not filtered

Author: e-sollier

pyjacker/ase.py

@@ -55,6 +55,7 @@ def compute_ase_matrix(samples,ase_dir,genes,genes_index=None,prior_coef=2.0,imp
 
             llr = llr_betabinom(df_sample.loc[x,"altCount"], df_sample.loc[x,"altCount"]+df_sample.loc[x,"refCount"])
             for g in genes_at_locus(genes_index,df_sample.loc[x,"contig"].lstrip("chr"),int(df_sample.loc[x,"position"])):
+                if not g.gene_id in genes: continue
                 if not g.gene_id in geneIDs2llrs: geneIDs2llrs[g.gene_id]=[]
                 geneIDs2llrs[g.gene_id].append(llr)
         for gene_id in geneIDs2llrs:

pyproject.toml

@@ -7,7 +7,7 @@ packages = ["pyjacker"]
 
 [project]
 name = 'pyjacker'
-version = "1.1.0"
+version = "1.1.1"
 description = 'Enhancer hijacking detection from WGS and RNAseq.'
 readme = 'README.md'
 authors = [