geney 1.3.58__py2.py3-none-any.whl → 1.3.60__py2.py3-none-any.whl

geney/SeqMats.py

@@ -300,6 +300,10 @@ class SeqMat:
     def raw_indices(self):
         return self.seqmat[self.ROW_INDS, :] + (self.seqmat[self.ROW_SUPERINDS, :] / 10)
 
+    @property
+    def filter(self):
+        return self.seqmat[:, self.seqmat[self.ROW_SEQ, :] != 5]
+
 
     def mutate(self, mut, inplace=False):
         """
@@ -357,15 +361,19 @@ class SeqMat:
         if tis_index not in self.indices:
             return SeqMat.from_seq({'seq': ''})
 
-        temp = self.seqmat[:, self._rel_index(tis_index):]
-        temp = temp[:, temp[0, :] != 5]
-        temp = SeqMat(temp)  # .drop_indices()
-        raw_seq = temp.seq  # Extract the raw sequence
+        raw_seq = SeqMat(self.seqmat[:, self._rel_index(tis_index):])
+        raw_seq = SeqMat(raw_seq.filter)
+
+        # temp = temp[:, temp[0, :] != 5]
+        # temp = SeqMat(temp)  # .drop_indices()
+        # raw_seq = temp.seq  # Extract the raw sequence
         pattern = re.compile(r"(?:[NACGT]{3})*?(TAA|TAG|TGA)")
-        # match = pattern.match(raw_seq)
-        match = [(m.group(), m.start()+1) for m in re.finditer(r"(?:[NACGT]{3})*?(TAA|TAG|TGA)", raw_seq) if
-                 (m.start()+1) % 3 == 0]
+        match = pattern.match(raw_seq.seq)
+        print(match)
+        # match = [(m.group(), m.start()+1) for m in re.finditer(r"(?:[NACGT]{3})*?(TAA|TAG|TGA)", raw_seq.seq)]# if
+        #          # (m.start()+1) % 3 == 0]
 
+        # print(match)
         if len(match) > 0:
             stop_index = match[0][-1] #.end()
         else:
@@ -373,7 +381,7 @@ class SeqMat:
 
         end_index = stop_index
         # assert (end_index+1) % 3 == 0, f"{end_index} is not a multiple of 3"
-        return SeqMat(temp.seqmat[:, :end_index])
+        return SeqMat(raw_seq.seqmat[:, :end_index])
 
     def translate(self, tis_index):
         """

{geney-1.3.58.dist-info → geney-1.3.60.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: geney
-Version: 1.3.58
+Version: 1.3.60
 Summary: A Python package for gene expression modeling.
 Home-page: https://github.com/nicolaslynn/geney
 Author: Nicolas Lynn

{geney-1.3.58.dist-info → geney-1.3.60.dist-info}/RECORD

@@ -1,6 +1,6 @@
 geney/Fasta_segment.py,sha256=99HxNGNh_MfdVW6hhtlb1vOn7eSmT7oFoEfHDFMxG8w,11275
 geney/Gene.py,sha256=nMWJjoQaiVFm2iRjoiq7ghZqnXtW0tJDcq2S0AyOIvY,6883
-geney/SeqMats.py,sha256=AKR1EhyDhOaEhavrkXSrVjYsohB6Qvsq11Udy7SGwvs,19238
+geney/SeqMats.py,sha256=DwEq_452FsF-zrr2nE1MMSTfKf54931yJERfXq0Ak7k,19460
 geney/Transcript.py,sha256=CpfxYkuCwFILozrtLuiWnlr1mRnMKn4o84HVJislgYs,14499
 geney/__init__.py,sha256=eBdDl42N6UhcYeZDjOnv199Z88fI5_8Y6xW8447OKXM,755
 geney/_mutation_utils.py,sha256=dHssUsnii_mf-wuRoMmF13UlD7k3ml_VwQMItTYnXpU,1132
@@ -25,7 +25,7 @@ geney/translation_initiation/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NM
 geney/translation_initiation/tis_utils.py,sha256=AF3siFjuQH-Rs44EV-80zHdbxRMvN4woLFSHroWIETc,4448
 geney/translation_initiation/resources/kozak_pssm.json,sha256=pcd0Olziutq-6H3mFWDCD9cujQ_AlZO-iiOvBl82hqE,1165
 geney/translation_initiation/resources/tis_regressor_model.joblib,sha256=IXb4DUDhJ5rBDKcqMk9zE3ECTZZcdj7Jixz3KpoZ7OA,2592025
-geney-1.3.58.dist-info/METADATA,sha256=oZ_vyEbmxOKk6EUGpM35weoldpehanRTYZ4bnQm64lA,990
-geney-1.3.58.dist-info/WHEEL,sha256=AHX6tWk3qWuce7vKLrj7lnulVHEdWoltgauo8bgCXgU,109
-geney-1.3.58.dist-info/top_level.txt,sha256=O-FuNUMb5fn9dhZ-dYCgF0aZtfi1EslMstnzhc5IIVo,6
-geney-1.3.58.dist-info/RECORD,,
+geney-1.3.60.dist-info/METADATA,sha256=9RdNxrkkxfNmokUDoatZq-ALXfGVcfyXXR5bEFI0_Ck,990
+geney-1.3.60.dist-info/WHEEL,sha256=AHX6tWk3qWuce7vKLrj7lnulVHEdWoltgauo8bgCXgU,109
+geney-1.3.60.dist-info/top_level.txt,sha256=O-FuNUMb5fn9dhZ-dYCgF0aZtfi1EslMstnzhc5IIVo,6
+geney-1.3.60.dist-info/RECORD,,