geney 1.3.65__py2.py3-none-any.whl → 1.3.66__py2.py3-none-any.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Potentially problematic release.
This version of geney might be problematic. Click here for more details.
geney/splicing_utils.py
CHANGED
|
@@ -166,7 +166,7 @@ def run_splicing_engine(seq, engine='spliceai'):
|
|
|
166
166
|
match engine:
|
|
167
167
|
case 'spliceai':
|
|
168
168
|
from .spliceai_utils import sai_predict_probs, sai_models
|
|
169
|
-
|
|
169
|
+
acceptor_probs, donor_probs = sai_predict_probs(seq, models=sai_models)
|
|
170
170
|
|
|
171
171
|
case 'pangolin':
|
|
172
172
|
from .pangolin_utils import pangolin_predict_probs, pang_models
|
|
@@ -214,6 +214,7 @@ def find_transcript_splicing(transcript, engine: str = 'spliceai') -> Tuple[Dict
|
|
|
214
214
|
# Create dictionaries and sort them by probability in descending order
|
|
215
215
|
donor_probs = dict(sorted(((i, p) for i, p in zip(ref_indices, ref_seq_donor_probs)),
|
|
216
216
|
key=lambda item: item[1], reverse=True))
|
|
217
|
+
|
|
217
218
|
acceptor_probs = dict(sorted(((i, p) for i, p in zip(ref_indices, ref_seq_acceptor_probs)),
|
|
218
219
|
key=lambda item: item[1], reverse=True))
|
|
219
220
|
|
|
@@ -332,8 +333,8 @@ def find_transcript_missplicing_seqs(ref_seq, var_seq, donors, acceptors, thresh
|
|
|
332
333
|
if ref_seq.seq == var_seq.seq:
|
|
333
334
|
return Missplicing({'missed_acceptors': {}, 'missed_donors': {}, 'discovered_acceptors': {}, 'discovered_donors': {}})
|
|
334
335
|
|
|
335
|
-
|
|
336
|
-
|
|
336
|
+
ref_seq_donor_probs, ref_seq_acceptor_probs = run_splicing_engine(ref_seq.seq, engine)
|
|
337
|
+
mut_seq_donor_probs, mut_seq_acceptor_probs = run_splicing_engine(var_seq.seq, engine)
|
|
337
338
|
ref_indices = ref_seq.indices[5000:-5000]
|
|
338
339
|
mut_indices = var_seq.indices[5000:-5000]
|
|
339
340
|
visible_donors = np.intersect1d(donors, ref_indices)
|
|
@@ -16,7 +16,7 @@ geney/pangolin_utils.py,sha256=9jdBXlOcRaUdfi-UpUxHA0AkTMZkUF-Lt7HVZ1nEm3s,2973
|
|
|
16
16
|
geney/power_utils.py,sha256=MehZFUdkJ2EFUot709yPEDxSkXmH5XevMebX2HD768A,7330
|
|
17
17
|
geney/seqmat_utils.py,sha256=wzb3PX5it5bpIFQvcxyzlxfhoJTbHHbsjg0rzh05iVs,19753
|
|
18
18
|
geney/spliceai_utils.py,sha256=tVY0T6F6l3fNoaktpn7Kq0oH5ZM0ThFYt9nPi_lfakw,3077
|
|
19
|
-
geney/splicing_utils.py,sha256=
|
|
19
|
+
geney/splicing_utils.py,sha256=_nXLCK41GhcrkXHXAqkhNV2IcwFltSxrR-rm8fUIrfE,40767
|
|
20
20
|
geney/survival_utils.py,sha256=KnAzEviMuXh6SnVXId9PgsFLSbgkduTvYoIthxN7FPA,6886
|
|
21
21
|
geney/tcga_utils.py,sha256=D_BNHm-D_K408dlcJm3hzH2c6QNFjQsKvUcOPiQRk7g,17612
|
|
22
22
|
geney/tis_utils.py,sha256=la0CZroaKe5RgAyFd4Bf_DqQncklWgAY2823xVst98o,7813
|
|
@@ -25,7 +25,7 @@ geney/translation_initiation/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NM
|
|
|
25
25
|
geney/translation_initiation/tis_utils.py,sha256=AF3siFjuQH-Rs44EV-80zHdbxRMvN4woLFSHroWIETc,4448
|
|
26
26
|
geney/translation_initiation/resources/kozak_pssm.json,sha256=pcd0Olziutq-6H3mFWDCD9cujQ_AlZO-iiOvBl82hqE,1165
|
|
27
27
|
geney/translation_initiation/resources/tis_regressor_model.joblib,sha256=IXb4DUDhJ5rBDKcqMk9zE3ECTZZcdj7Jixz3KpoZ7OA,2592025
|
|
28
|
-
geney-1.3.
|
|
29
|
-
geney-1.3.
|
|
30
|
-
geney-1.3.
|
|
31
|
-
geney-1.3.
|
|
28
|
+
geney-1.3.66.dist-info/METADATA,sha256=bl8lWCBcJsbfBPJmkoY8xG0n6G7z7X1C-6jA1bSevCk,990
|
|
29
|
+
geney-1.3.66.dist-info/WHEEL,sha256=AHX6tWk3qWuce7vKLrj7lnulVHEdWoltgauo8bgCXgU,109
|
|
30
|
+
geney-1.3.66.dist-info/top_level.txt,sha256=O-FuNUMb5fn9dhZ-dYCgF0aZtfi1EslMstnzhc5IIVo,6
|
|
31
|
+
geney-1.3.66.dist-info/RECORD,,
|
|
File without changes
|
|
File without changes
|