mlbi-bitool 0.2.2__tar.gz → 0.2.3__tar.gz

{mlbi_bitool-0.2.2 → mlbi_bitool-0.2.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: mlbi-bitool
-Version: 0.2.2
+Version: 0.2.3
 Summary: Python interface to BioInfo tools for use in Google.Colab
 Author-email: MLBI Lab <syoon@dku.edu>
 License-Expression: MIT

{mlbi_bitool-0.2.2 → mlbi_bitool-0.2.3}/pyproject.toml

@@ -6,7 +6,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "mlbi-bitool"
-version = "0.2.2"
+version = "0.2.3"
 description = "Python interface to BioInfo tools for use in Google.Colab"
 readme = "README.md"
 authors = [{ name = "MLBI Lab", email = "syoon@dku.edu" }]

{mlbi_bitool-0.2.2 → mlbi_bitool-0.2.3}/src/mlbi_bitool/tasm.py

@@ -1,7 +1,7 @@
 import re, copy, collections
 import numpy as np
 import pandas as pd
-import os, time
+import os, time, glob
 from subprocess import Popen, PIPE
 import shlex, shutil, tempfile, sys
 from typing import Optional, Tuple, Set
@@ -701,11 +701,12 @@ def Stringtie_merge_n_compare( src_dir: str,
                                verbose = True ):
 
     cmd_lst = []
+    if src_dir[-1] == '/':
+        src_dir = src_dir[:-1]
     flst = os.listdir(src_dir)
-    flst = [src_dir + f for f in flst if f.split('.')[-1] == 'gff']
+    flst = [os.path.join(src_dir, f) for f in flst if f.split('.')[-1] == 'gff']
     flst.sort()
 
-    if src_dir[-1] == '/': src_dir = src_dir[:-1]
     if out_dir is None:
         out_dir = src_dir + '_merged_info'
     if not os.path.isdir(out_dir): os.mkdir(out_dir)
@@ -752,7 +753,23 @@ def Stringtie_merge_n_compare( src_dir: str,
     if verbose: print('done.')
 
     if verbose: print('Applying to merged GFF ... ', end = '')
-    df_tmap = pd.read_csv( gffcomp_prefix + '.merged.gff.tmap', sep = '\t' )
+    tmap_file = gffcomp_prefix + '.merged.gff.tmap'
+    if not os.path.isfile(tmap_file):
+        tmap_candidates = sorted(glob.glob(gffcomp_prefix + '*.tmap'))
+        if len(tmap_candidates) == 1:
+            tmap_file = tmap_candidates[0]
+        elif len(tmap_candidates) > 1:
+            tmap_candidates.sort(key=lambda p: ('.merged.gff.tmap' not in p, len(p), p))
+            tmap_file = tmap_candidates[0]
+        else:
+            existing = sorted(glob.glob(gffcomp_prefix + '*'))
+            msg = (
+                "gffcompare did not create a .tmap file. "
+                f"Expected {gffcomp_prefix}.merged.gff.tmap or {gffcomp_prefix}*.tmap. "
+                f"Exit code: {exit_code}. Existing output files: {existing}"
+            )
+            raise FileNotFoundError(msg)
+    df_tmap = pd.read_csv( tmap_file, sep = '\t' )
     df_gff = reflect_gffcompare_res_to_stringtie_res( df_gtf, df_gff, df_tmap, 
                                                       class_code_to_drop = class_code_to_drop, 
                                                       verbose = False )

{mlbi_bitool-0.2.2 → mlbi_bitool-0.2.3}/src/mlbi_bitool.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: mlbi-bitool
-Version: 0.2.2
+Version: 0.2.3
 Summary: Python interface to BioInfo tools for use in Google.Colab
 Author-email: MLBI Lab <syoon@dku.edu>
 License-Expression: MIT