gsrap 0.8.2__py3-none-any.whl → 0.8.3__py3-none-any.whl

gsrap/.ipynb_checkpoints/__init__-checkpoint.py

@@ -75,6 +75,7 @@ def main():
     parsedb_parser.add_argument("-z", "--initialize", metavar='', type=str, default='-', help="Initialize the universe on the provided medium. By default, the first medium in --media is used. Use 'none' to avoid initialization.")
     parsedb_parser.add_argument("--precursors", action='store_true', help="Verify biosynthesis of biomass precursors and show blocked ones.")
     parsedb_parser.add_argument("--biosynth", action='store_true', help="Check biosynthesis of all metabolites and detect dead-ends.")
+    parsedb_parser.add_argument("-t", "--taxon", metavar='', type=str, default='-', help="High-level taxon of interest. If provided, it must follow the syntax '{level}:{name}', where {level} is 'kingdom' or 'phylum'.")
     parsedb_parser.add_argument("-e", "--eggnog", nargs='+', metavar='', type=str, default='-', help="Path to the optional eggnog-mapper annotation table(s).")
     parsedb_parser.add_argument("-k", "--keggorg", metavar='', type=str, default='-', help="A single KEGG Organism code. If provided, it takes precedence over --eggnog.")
     parsedb_parser.add_argument("--goodbefore", metavar='', type=str, default='-', help="Syntax is {pure_mid}-{rid1}-{rid2}. From top to bottom, build the universe until reaction {rid1}, transport {rid2} and metabolite {pure_mid} are reached.")

gsrap/__init__.py

@@ -75,6 +75,7 @@ def main():
     parsedb_parser.add_argument("-z", "--initialize", metavar='', type=str, default='-', help="Initialize the universe on the provided medium. By default, the first medium in --media is used. Use 'none' to avoid initialization.")
     parsedb_parser.add_argument("--precursors", action='store_true', help="Verify biosynthesis of biomass precursors and show blocked ones.")
     parsedb_parser.add_argument("--biosynth", action='store_true', help="Check biosynthesis of all metabolites and detect dead-ends.")
+    parsedb_parser.add_argument("-t", "--taxon", metavar='', type=str, default='-', help="High-level taxon of interest. If provided, it must follow the syntax '{level}:{name}', where {level} is 'kingdom' or 'phylum'.")
     parsedb_parser.add_argument("-e", "--eggnog", nargs='+', metavar='', type=str, default='-', help="Path to the optional eggnog-mapper annotation table(s).")
     parsedb_parser.add_argument("-k", "--keggorg", metavar='', type=str, default='-', help="A single KEGG Organism code. If provided, it takes precedence over --eggnog.")
     parsedb_parser.add_argument("--goodbefore", metavar='', type=str, default='-', help="Syntax is {pure_mid}-{rid1}-{rid2}. From top to bottom, build the universe until reaction {rid1}, transport {rid2} and metabolite {pure_mid} are reached.")

gsrap/commons/.ipynb_checkpoints/downloads-checkpoint.py

@@ -243,7 +243,38 @@ def format_expansion(logger, eggnog):
     
 
 
+def check_taxon(logger, taxon, idcollection_dict):
+    
+    
+    # verify presence of needed assets
+    if 'ko_to_taxa' not in idcollection_dict.keys():
+        logger.error(f"Asset 'ko_to_taxa' not found in 'gsrap.maps'. Please update 'gsrap.maps' with 'gsrap getmaps'.")
+        return 1
+    
+    
+    # extract level and name
+    try: level, name = taxon.split(':')
+    except: 
+        logger.error(f"Provided --taxon is not well formatted: '{taxon}'.")
+        return 1
+    
+    
+    # compute available levels and check
+    avail_levels = set(['kingdom', 'phylum'])            
+    if level not in avail_levels:
+        logger.error(f"Provided level is not acceptable: '{level}' (see --taxon). Acceptable levels are {avail_levels}.")
+        return 1
+    
+    
+    # compute available taxa at input level
+    avail_taxa_at_level = set()
+    ko_to_taxa = idcollection_dict['ko_to_taxa']
+    for ko in ko_to_taxa.keys():
+        for taxon_name in ko_to_taxa[ko][level]:
+            avail_taxa_at_level.add(taxon_name)
+    if name not in avail_taxa_at_level:
+        logger.error(f"Provided taxon name is not acceptable: '{name}' (see --taxon). Acceptable taxon names for level '{level}' are {avail_taxa_at_level}.")
+        return 1
 
-
-
-
+    
+    return 0

gsrap/commons/.ipynb_checkpoints/excelhub-checkpoint.py

@@ -148,7 +148,7 @@ def write_excel_model(model, filepath, nofigs, memote_results_dict, df_E, df_B,
             else: df_T.append(row_dict)
             
     for g in model.genes:
-        row_dict = {'gid': g.id, 'involved_in': '; '.join([r.id for r in g.reactions])}
+        row_dict = {'gid': g.id, 'name': g.name, 'involved_in': '; '.join([r.id for r in g.reactions])}
         
         for db in g.annotation.keys():
             annots = g.annotation[db]
@@ -171,7 +171,7 @@ def write_excel_model(model, filepath, nofigs, memote_results_dict, df_E, df_B,
     df_R = df_R[df_R_first_cols + sorted([c for c in df_R.columns if c not in df_R_first_cols])]
     df_T = df_T[df_R_first_cols + sorted([c for c in df_T.columns if c not in df_R_first_cols])]
     df_A = df_A[df_R_first_cols + sorted([c for c in df_A.columns if c not in df_R_first_cols])]
-    df_G_first_cols = ['gid', 'involved_in']
+    df_G_first_cols = ['gid', 'name', 'involved_in']
     df_G = df_G[df_G_first_cols + sorted([c for c in df_G.columns if c not in df_G_first_cols])]
     
     

gsrap/commons/downloads.py

@@ -243,7 +243,38 @@ def format_expansion(logger, eggnog):
     
 
 
+def check_taxon(logger, taxon, idcollection_dict):
+    
+    
+    # verify presence of needed assets
+    if 'ko_to_taxa' not in idcollection_dict.keys():
+        logger.error(f"Asset 'ko_to_taxa' not found in 'gsrap.maps'. Please update 'gsrap.maps' with 'gsrap getmaps'.")
+        return 1
+    
+    
+    # extract level and name
+    try: level, name = taxon.split(':')
+    except: 
+        logger.error(f"Provided --taxon is not well formatted: '{taxon}'.")
+        return 1
+    
+    
+    # compute available levels and check
+    avail_levels = set(['kingdom', 'phylum'])            
+    if level not in avail_levels:
+        logger.error(f"Provided level is not acceptable: '{level}' (see --taxon). Acceptable levels are {avail_levels}.")
+        return 1
+    
+    
+    # compute available taxa at input level
+    avail_taxa_at_level = set()
+    ko_to_taxa = idcollection_dict['ko_to_taxa']
+    for ko in ko_to_taxa.keys():
+        for taxon_name in ko_to_taxa[ko][level]:
+            avail_taxa_at_level.add(taxon_name)
+    if name not in avail_taxa_at_level:
+        logger.error(f"Provided taxon name is not acceptable: '{name}' (see --taxon). Acceptable taxon names for level '{level}' are {avail_taxa_at_level}.")
+        return 1
 
-
-
-
+    
+    return 0

gsrap/commons/excelhub.py

@@ -148,7 +148,7 @@ def write_excel_model(model, filepath, nofigs, memote_results_dict, df_E, df_B,
             else: df_T.append(row_dict)
             
     for g in model.genes:
-        row_dict = {'gid': g.id, 'involved_in': '; '.join([r.id for r in g.reactions])}
+        row_dict = {'gid': g.id, 'name': g.name, 'involved_in': '; '.join([r.id for r in g.reactions])}
         
         for db in g.annotation.keys():
             annots = g.annotation[db]
@@ -171,7 +171,7 @@ def write_excel_model(model, filepath, nofigs, memote_results_dict, df_E, df_B,
     df_R = df_R[df_R_first_cols + sorted([c for c in df_R.columns if c not in df_R_first_cols])]
     df_T = df_T[df_R_first_cols + sorted([c for c in df_T.columns if c not in df_R_first_cols])]
     df_A = df_A[df_R_first_cols + sorted([c for c in df_A.columns if c not in df_R_first_cols])]
-    df_G_first_cols = ['gid', 'involved_in']
+    df_G_first_cols = ['gid', 'name', 'involved_in']
     df_G = df_G[df_G_first_cols + sorted([c for c in df_G.columns if c not in df_G_first_cols])]
     
     

gsrap/getmaps/.ipynb_checkpoints/getmaps-checkpoint.py

@@ -4,6 +4,7 @@ import pickle
 
 
 from .kdown import download_raw_txtfiles
+from .kdown import create_dict_keggorg
 from .kdown import create_dict_ko
 from .kdown import create_dict_c
 from .kdown import create_dict_r
@@ -20,13 +21,19 @@ def do_kdown(logger, outdir, usecache, keeptmp):
     logger.info(f"Respectfully retrieving metabolic information from KEGG. Raw data are being saved into '{outdir}/kdown/'. Be patient, could take a couple of days...")
     os.makedirs(f'{outdir}/kdown/', exist_ok=True)
     
+    
     response = download_raw_txtfiles(logger, outdir, usecache)
     if type(response) == int: return 1
     else: RELEASE_kegg = response
     
+        
     
     logger.info("Parsing downloaded KEGG information...")
-
+    
+    response = create_dict_keggorg(logger, outdir)
+    if type(response) == int: return 1
+    else: dict_keggorg = response
+    
     response = create_dict_ko(logger, outdir)
     if type(response) == int: return 1
     else: dict_ko = response
@@ -49,7 +56,7 @@ def do_kdown(logger, outdir, usecache, keeptmp):
     
     
     # create 'idcollection_dict' and 'summary_dict' dictionaries
-    idcollection_dict = create_idcollection_dict(dict_ko, dict_c, dict_r, dict_map, dict_md)
+    idcollection_dict = create_idcollection_dict(dict_keggorg, dict_ko, dict_c, dict_r, dict_map, dict_md)
     summary_dict = create_summary_dict(dict_c, dict_r, dict_map, dict_md)
     
         
@@ -57,7 +64,6 @@ def do_kdown(logger, outdir, usecache, keeptmp):
 
 
 
-
 def main(args, logger):
     
     
@@ -67,7 +73,7 @@ def main(args, logger):
     os.makedirs(f'{args.outdir}/', exist_ok=True)
     
     
-    # KEGG
+    # KEGG download
     response = do_kdown(logger, args.outdir, args.usecache, args.keeptmp)
     if type(response) == int: return 1
     else: RELEASE_kegg, idcollection_dict, summary_dict = response[0], response[1], response[2]
@@ -76,7 +82,9 @@ def main(args, logger):
     # create 'gsrap.maps':
     with open(f'{args.outdir}/gsrap.maps', 'wb') as wb_handler:
         pickle.dump({
-            'RELEASE_kegg': RELEASE_kegg, 'idcollection_dict': idcollection_dict, 'summary_dict': summary_dict,
+            'RELEASE_kegg': RELEASE_kegg, 
+            'idcollection_dict': idcollection_dict, 
+            'summary_dict': summary_dict, 
         }, wb_handler)
     logger.info(f"'{args.outdir}/gsrap.maps' created!")
     
@@ -87,4 +95,5 @@ def main(args, logger):
         logger.info(f"Temporary raw files deleted!")
     
     
+    
     return 0

gsrap/getmaps/.ipynb_checkpoints/kdown-checkpoint.py

@@ -34,6 +34,7 @@ def download_raw_txtfiles(logger, outdir, usecache):
         'orthology',
         'module', 
         'pathway',
+        'organism',
     ]
     for db in databases:
         time.sleep(0.5)
@@ -45,8 +46,9 @@ def download_raw_txtfiles(logger, outdir, usecache):
 
     # mix the items to download to be respectful/compliant
     items_to_download = []
-
     for db in databases:
+        if db == 'organism': 
+            continue   # here we just need the list
         with open(f"{outdir}/kdown/{db}.txt", 'r') as file:
             res_string = file.read()
             rows = res_string.split('\n')
@@ -54,7 +56,6 @@ def download_raw_txtfiles(logger, outdir, usecache):
                 item_id = row.split('\t', 1)[0]
                 if item_id == '': continue
                 items_to_download.append({'db': db, 'id': item_id})
-
     random.shuffle(items_to_download)
 
 
@@ -79,6 +80,51 @@ def download_raw_txtfiles(logger, outdir, usecache):
 
 
 
+def create_dict_keggorg(logger, outdir):
+    
+    organisms_raw = open(f'{outdir}/kdown/organism.txt', 'r').read()
+
+    # create a dataframe listing all organisms in KEGG;
+    # columns are [tnumber, name, domain, kingdom, phylum, classification]
+    df = []  # list fo dicts
+    for line in organisms_raw.strip().split("\n"):
+        fields = line.split("\t")
+        if len(fields) == 4:
+            tnumber, keggorg, name, classification = fields
+            levels = classification.split(";")
+            domain = levels[0]
+            kingdom = levels[1]
+            phylum = levels[2]
+            df.append({
+                'tnumber':tnumber, 
+                'keggorg': keggorg, 
+                'name': name, 
+                'domain': domain, 
+                'kingdom': kingdom,
+                'phylum': phylum,
+                'classification': classification
+            })
+        else: 
+            # never verified during tests!
+            logger.warning(f'Strange number of fields found in this line of "organism.txt": """{line}""".')
+    df = pnd.DataFrame.from_records(df)
+    df = df.set_index('keggorg', drop=True, verify_integrity=True)
+    
+    
+    # convert dataframe to dict
+    dict_keggorg = {}
+    for keggorg, row in df.iterrows(): 
+        dict_keggorg[keggorg] = {
+            'kingdom': row['kingdom'],
+            'phylum': row['phylum'],
+            #'name': row['name'],   # not strictly needed. Commented to save disk space.
+        }
+        
+    if logger != None: logger.info(f'Number of unique items (org): {len(dict_keggorg.keys())}.')
+    return dict_keggorg
+
+
+
 def create_dict_ko(logger, outdir):
     
     dict_ko = {}         # main output
@@ -98,6 +144,7 @@ def create_dict_ko(logger, outdir):
                 'ecs': set(),
                 'cogs': set(),
                 'gos': set(),
+                'keggorgs': set(),
             }   
         else:
             logger.error(f"{ko_id} already included!")
@@ -175,7 +222,13 @@ def create_dict_ko(logger, outdir):
                         gos = content[len('GO: '):].strip().split(' ')
                         for go in gos:
                             dict_ko[ko_id]['gos'].add(go)
-                        
+                            
+                            
+                # parse the organism-specific genes
+                if curr_header == 'GENES       ':
+                    keggorg = content.split(': ',1)[0]
+                    dict_ko[ko_id]['keggorgs'].add(keggorg.lower()) # organism.txt has IDs in lowercase
+                    
                         
                 # parse the reactions
                 if curr_header == 'REACTION    ':
@@ -547,7 +600,7 @@ def create_dict_md(logger, outdir):
 
 
 
-def create_idcollection_dict(dict_ko, dict_c, dict_r, dict_map, dict_md):
+def create_idcollection_dict(dict_keggorg, dict_ko, dict_c, dict_r, dict_map, dict_md):
 
     idcollection_dict = {}
     
@@ -620,6 +673,24 @@ def create_idcollection_dict(dict_ko, dict_c, dict_r, dict_map, dict_md):
         for go in dict_ko[ko_id]['gos']:
             idcollection_dict['ko_to_gos'][ko_id].add(go)
             
+    
+    # creation of 'ko_to_keggorgs' skipped as it takes too much disk space. Replaced with 'ko_to_taxa'. 
+    idcollection_dict['ko_to_taxa'] = {}
+    missing_keggorgs = set()
+    for ko_id in dict_ko.keys():
+        idcollection_dict['ko_to_taxa'][ko_id] = {'kingdom': set(), 'phylum': set()}
+        for keggorg in dict_ko[ko_id]['keggorgs']:
+            try: 
+                kingdom = dict_keggorg[keggorg]['kingdom']
+                phylum = dict_keggorg[keggorg]['phylum']
+            except: 
+                if keggorg not in missing_keggorgs:
+                    missing_keggorgs.add(keggorg)
+                    #print(f"Organism '{keggorg}' appears in 'orthology/' but not in 'organism.txt'.")
+                continue
+            idcollection_dict['ko_to_taxa'][ko_id]['kingdom'].add(kingdom)
+            idcollection_dict['ko_to_taxa'][ko_id]['phylum'].add(phylum)
+
             
     idcollection_dict['map_to_name'] = {}
     for map_id in dict_map.keys():

gsrap/getmaps/getmaps.py

@@ -4,6 +4,7 @@ import pickle
 
 
 from .kdown import download_raw_txtfiles
+from .kdown import create_dict_keggorg
 from .kdown import create_dict_ko
 from .kdown import create_dict_c
 from .kdown import create_dict_r
@@ -20,13 +21,19 @@ def do_kdown(logger, outdir, usecache, keeptmp):
     logger.info(f"Respectfully retrieving metabolic information from KEGG. Raw data are being saved into '{outdir}/kdown/'. Be patient, could take a couple of days...")
     os.makedirs(f'{outdir}/kdown/', exist_ok=True)
     
+    
     response = download_raw_txtfiles(logger, outdir, usecache)
     if type(response) == int: return 1
     else: RELEASE_kegg = response
     
+        
     
     logger.info("Parsing downloaded KEGG information...")
-
+    
+    response = create_dict_keggorg(logger, outdir)
+    if type(response) == int: return 1
+    else: dict_keggorg = response
+    
     response = create_dict_ko(logger, outdir)
     if type(response) == int: return 1
     else: dict_ko = response
@@ -49,7 +56,7 @@ def do_kdown(logger, outdir, usecache, keeptmp):
     
     
     # create 'idcollection_dict' and 'summary_dict' dictionaries
-    idcollection_dict = create_idcollection_dict(dict_ko, dict_c, dict_r, dict_map, dict_md)
+    idcollection_dict = create_idcollection_dict(dict_keggorg, dict_ko, dict_c, dict_r, dict_map, dict_md)
     summary_dict = create_summary_dict(dict_c, dict_r, dict_map, dict_md)
     
         
@@ -57,7 +64,6 @@ def do_kdown(logger, outdir, usecache, keeptmp):
 
 
 
-
 def main(args, logger):
     
     
@@ -67,7 +73,7 @@ def main(args, logger):
     os.makedirs(f'{args.outdir}/', exist_ok=True)
     
     
-    # KEGG
+    # KEGG download
     response = do_kdown(logger, args.outdir, args.usecache, args.keeptmp)
     if type(response) == int: return 1
     else: RELEASE_kegg, idcollection_dict, summary_dict = response[0], response[1], response[2]
@@ -76,7 +82,9 @@ def main(args, logger):
     # create 'gsrap.maps':
     with open(f'{args.outdir}/gsrap.maps', 'wb') as wb_handler:
         pickle.dump({
-            'RELEASE_kegg': RELEASE_kegg, 'idcollection_dict': idcollection_dict, 'summary_dict': summary_dict,
+            'RELEASE_kegg': RELEASE_kegg, 
+            'idcollection_dict': idcollection_dict, 
+            'summary_dict': summary_dict, 
         }, wb_handler)
     logger.info(f"'{args.outdir}/gsrap.maps' created!")
     
@@ -87,4 +95,5 @@ def main(args, logger):
         logger.info(f"Temporary raw files deleted!")
     
     
+    
     return 0

gsrap/getmaps/kdown.py

@@ -34,6 +34,7 @@ def download_raw_txtfiles(logger, outdir, usecache):
         'orthology',
         'module', 
         'pathway',
+        'organism',
     ]
     for db in databases:
         time.sleep(0.5)
@@ -45,8 +46,9 @@ def download_raw_txtfiles(logger, outdir, usecache):
 
     # mix the items to download to be respectful/compliant
     items_to_download = []
-
     for db in databases:
+        if db == 'organism': 
+            continue   # here we just need the list
         with open(f"{outdir}/kdown/{db}.txt", 'r') as file:
             res_string = file.read()
             rows = res_string.split('\n')
@@ -54,7 +56,6 @@ def download_raw_txtfiles(logger, outdir, usecache):
                 item_id = row.split('\t', 1)[0]
                 if item_id == '': continue
                 items_to_download.append({'db': db, 'id': item_id})
-
     random.shuffle(items_to_download)
 
 
@@ -79,6 +80,51 @@ def download_raw_txtfiles(logger, outdir, usecache):
 
 
 
+def create_dict_keggorg(logger, outdir):
+    
+    organisms_raw = open(f'{outdir}/kdown/organism.txt', 'r').read()
+
+    # create a dataframe listing all organisms in KEGG;
+    # columns are [tnumber, name, domain, kingdom, phylum, classification]
+    df = []  # list fo dicts
+    for line in organisms_raw.strip().split("\n"):
+        fields = line.split("\t")
+        if len(fields) == 4:
+            tnumber, keggorg, name, classification = fields
+            levels = classification.split(";")
+            domain = levels[0]
+            kingdom = levels[1]
+            phylum = levels[2]
+            df.append({
+                'tnumber':tnumber, 
+                'keggorg': keggorg, 
+                'name': name, 
+                'domain': domain, 
+                'kingdom': kingdom,
+                'phylum': phylum,
+                'classification': classification
+            })
+        else: 
+            # never verified during tests!
+            logger.warning(f'Strange number of fields found in this line of "organism.txt": """{line}""".')
+    df = pnd.DataFrame.from_records(df)
+    df = df.set_index('keggorg', drop=True, verify_integrity=True)
+    
+    
+    # convert dataframe to dict
+    dict_keggorg = {}
+    for keggorg, row in df.iterrows(): 
+        dict_keggorg[keggorg] = {
+            'kingdom': row['kingdom'],
+            'phylum': row['phylum'],
+            #'name': row['name'],   # not strictly needed. Commented to save disk space.
+        }
+        
+    if logger != None: logger.info(f'Number of unique items (org): {len(dict_keggorg.keys())}.')
+    return dict_keggorg
+
+
+
 def create_dict_ko(logger, outdir):
     
     dict_ko = {}         # main output
@@ -98,6 +144,7 @@ def create_dict_ko(logger, outdir):
                 'ecs': set(),
                 'cogs': set(),
                 'gos': set(),
+                'keggorgs': set(),
             }   
         else:
             logger.error(f"{ko_id} already included!")
@@ -175,7 +222,13 @@ def create_dict_ko(logger, outdir):
                         gos = content[len('GO: '):].strip().split(' ')
                         for go in gos:
                             dict_ko[ko_id]['gos'].add(go)
-                        
+                            
+                            
+                # parse the organism-specific genes
+                if curr_header == 'GENES       ':
+                    keggorg = content.split(': ',1)[0]
+                    dict_ko[ko_id]['keggorgs'].add(keggorg.lower()) # organism.txt has IDs in lowercase
+                    
                         
                 # parse the reactions
                 if curr_header == 'REACTION    ':
@@ -547,7 +600,7 @@ def create_dict_md(logger, outdir):
 
 
 
-def create_idcollection_dict(dict_ko, dict_c, dict_r, dict_map, dict_md):
+def create_idcollection_dict(dict_keggorg, dict_ko, dict_c, dict_r, dict_map, dict_md):
 
     idcollection_dict = {}
     
@@ -620,6 +673,24 @@ def create_idcollection_dict(dict_ko, dict_c, dict_r, dict_map, dict_md):
         for go in dict_ko[ko_id]['gos']:
             idcollection_dict['ko_to_gos'][ko_id].add(go)
             
+    
+    # creation of 'ko_to_keggorgs' skipped as it takes too much disk space. Replaced with 'ko_to_taxa'. 
+    idcollection_dict['ko_to_taxa'] = {}
+    missing_keggorgs = set()
+    for ko_id in dict_ko.keys():
+        idcollection_dict['ko_to_taxa'][ko_id] = {'kingdom': set(), 'phylum': set()}
+        for keggorg in dict_ko[ko_id]['keggorgs']:
+            try: 
+                kingdom = dict_keggorg[keggorg]['kingdom']
+                phylum = dict_keggorg[keggorg]['phylum']
+            except: 
+                if keggorg not in missing_keggorgs:
+                    missing_keggorgs.add(keggorg)
+                    #print(f"Organism '{keggorg}' appears in 'orthology/' but not in 'organism.txt'.")
+                continue
+            idcollection_dict['ko_to_taxa'][ko_id]['kingdom'].add(kingdom)
+            idcollection_dict['ko_to_taxa'][ko_id]['phylum'].add(phylum)
+
             
     idcollection_dict['map_to_name'] = {}
     for map_id in dict_map.keys():

gsrap/parsedb/.ipynb_checkpoints/completeness-checkpoint.py

@@ -27,9 +27,8 @@ def parse_eggnog(model, eggnog, idcollection_dict):
             
             
     # PART 2. get reactions in the organism (even the GPR is not complete)
-    kr_to_kos = idcollection_dict['kr_to_kos']
     krs_org = set()
-    for kr, kos in kr_to_kos.items(): 
+    for kr, kos in idcollection_dict['kr_to_kos'].items(): 
         if any([ko in kos_org for ko in kos]):
             krs_org.add(kr)
     
@@ -49,9 +48,34 @@ def parse_keggorg(keggorg, outdir, idcollection_dict):
     
     
     # PART 2. get reactions in the organism (even the GPR is not complete)
-    kr_to_kos = idcollection_dict['kr_to_kos']
     krs_org = set()
-    for kr, kos in kr_to_kos.items(): 
+    for kr, kos in idcollection_dict['kr_to_kos'].items(): 
+        if any([ko in kos_org for ko in kos]):
+            krs_org.add(kr)
+    
+    
+    return krs_org
+
+
+
+def parse_taxon(taxon, idcollection_dict):
+    
+    
+    # formatting of --taxon was already verified at startup.
+    # also the presence of 'ko_to_taxa' in idcollection_dict was veryfied at startup.
+    level, name = taxon.split(':')
+           
+        
+    # PART 1. get KO codes available
+    kos_org = set()
+    for ko in idcollection_dict['ko_to_taxa'].keys():
+        if name in idcollection_dict['ko_to_taxa'][ko][level]:
+            kos_org.add(ko)
+            
+    
+    # PART 2. get reactions in the organism (even the GPR is not complete)
+    krs_org = set()
+    for kr, kos in idcollection_dict['kr_to_kos'].items(): 
         if any([ko in kos_org for ko in kos]):
             krs_org.add(kr)
     
@@ -60,7 +84,7 @@ def parse_keggorg(keggorg, outdir, idcollection_dict):
 
     
     
-def check_completeness(logger, model, progress, module, focus, eggnog, keggorg, idcollection_dict, summary_dict, outdir): 
+def check_completeness(logger, model, progress, module, focus, taxon, eggnog, keggorg, idcollection_dict, summary_dict, outdir): 
     # check KEGG annotations in the universe model to get '%' of completeness per pathway/module.
     
             
@@ -69,6 +93,9 @@ def check_completeness(logger, model, progress, module, focus, eggnog, keggorg,
     if keggorg != '-':  # keggorg has precedence
         kr_uni = parse_keggorg(keggorg, outdir, idcollection_dict)
         kr_uni_label = f"organism code '{keggorg}'"
+    elif taxon != '-':
+        kr_uni = parse_taxon(taxon, idcollection_dict)
+        kr_uni_label = f"taxon '{taxon}'"
     elif eggnog != '-':
         for eggfile in eggnog:
             eggset = parse_eggnog(model, eggfile, idcollection_dict)

gsrap/parsedb/.ipynb_checkpoints/parsedb-checkpoint.py

@@ -17,6 +17,7 @@ from ..commons import show_contributions
 from ..commons import adjust_biomass_precursors
 from ..commons import count_undrawn_rids
 from ..commons import format_expansion
+from ..commons import check_taxon
 from ..commons import download_keggorg
 from ..commons import initialize_model
 from ..commons import get_memote_results_dict
@@ -46,6 +47,7 @@ from .cycles import verify_egc_all
 def main(args, logger):
     
     
+    ###### FORMAT ARGS NOT REQUIRING RESOURCES
     # adjust out folder path                      
     while args.outdir.endswith('/'):              
         args.outdir = args.outdir[:-1]
@@ -77,17 +79,8 @@ def main(args, logger):
     if args.onlyauthor == '-': args.onlyauthor = None
     
     
-    # format the --eggnog param
-    args.eggnog = format_expansion(logger, args.eggnog)  # now 'args.eggnog' could still be '-'
-    
-    # get the kegg organism if requested
-    if args.keggorg != '-':
-        response = download_keggorg(logger, args.keggorg, args.outdir)
-        if response == 1: return 1
-    
-    
-    
     
+    ###### LOAD LOCAL RESOURCES
     # check and extract the required 'gsrap.maps' file
     if os.path.exists(f'{args.inmaps}') == False:
         logger.error(f"File 'gsrap.maps' not found at {args.inmaps}.")
@@ -108,9 +101,27 @@ def main(args, logger):
             kegg_compound_to_others = pickle.load(handle)
     with resources.path("gsrap.assets", f"kegg_reaction_to_others.pickle") as asset_path: 
         with open(asset_path, 'rb') as handle:
-            kegg_reaction_to_others = pickle.load(handle)   
+            kegg_reaction_to_others = pickle.load(handle) 
+            
+            
+            
+    ###### FORMAT/CHECK FOCUSING ARGS
+    # format the --eggnog param
+    args.eggnog = format_expansion(logger, args.eggnog)  # now 'args.eggnog' could still be '-'
+    
+    # check the --taxon param
+    if args.taxon != '-':
+        response = check_taxon(logger, args.taxon, idcollection_dict)
+        if response == 1: return 1
+    
+    # get the kegg organism if requested
+    if args.keggorg != '-':
+        response = download_keggorg(logger, args.keggorg, args.outdir)
+        if response == 1: return 1
+    
     
     
+    # DOWNLOAD ONLINE RESOURCES
     # get dbuni and dbexp:
     logger.info("Downloading gsrap database...")
     response = get_databases(logger)
@@ -166,7 +177,7 @@ def main(args, logger):
     
     ###### CHECKS 1
     # check universe completness
-    df_C = check_completeness(logger, universe, args.progress, args.module, args.focus, args.eggnog, args.keggorg, idcollection_dict, summary_dict, args.outdir)
+    df_C = check_completeness(logger, universe, args.progress, args.module, args.focus, args.taxon, args.eggnog, args.keggorg, idcollection_dict, summary_dict, args.outdir)
     if type(df_C)==int: return 1
 
 

gsrap/parsedb/completeness.py

@@ -27,9 +27,8 @@ def parse_eggnog(model, eggnog, idcollection_dict):
             
             
     # PART 2. get reactions in the organism (even the GPR is not complete)
-    kr_to_kos = idcollection_dict['kr_to_kos']
     krs_org = set()
-    for kr, kos in kr_to_kos.items(): 
+    for kr, kos in idcollection_dict['kr_to_kos'].items(): 
         if any([ko in kos_org for ko in kos]):
             krs_org.add(kr)
     
@@ -49,9 +48,34 @@ def parse_keggorg(keggorg, outdir, idcollection_dict):
     
     
     # PART 2. get reactions in the organism (even the GPR is not complete)
-    kr_to_kos = idcollection_dict['kr_to_kos']
     krs_org = set()
-    for kr, kos in kr_to_kos.items(): 
+    for kr, kos in idcollection_dict['kr_to_kos'].items(): 
+        if any([ko in kos_org for ko in kos]):
+            krs_org.add(kr)
+    
+    
+    return krs_org
+
+
+
+def parse_taxon(taxon, idcollection_dict):
+    
+    
+    # formatting of --taxon was already verified at startup.
+    # also the presence of 'ko_to_taxa' in idcollection_dict was veryfied at startup.
+    level, name = taxon.split(':')
+           
+        
+    # PART 1. get KO codes available
+    kos_org = set()
+    for ko in idcollection_dict['ko_to_taxa'].keys():
+        if name in idcollection_dict['ko_to_taxa'][ko][level]:
+            kos_org.add(ko)
+            
+    
+    # PART 2. get reactions in the organism (even the GPR is not complete)
+    krs_org = set()
+    for kr, kos in idcollection_dict['kr_to_kos'].items(): 
         if any([ko in kos_org for ko in kos]):
             krs_org.add(kr)
     
@@ -60,7 +84,7 @@ def parse_keggorg(keggorg, outdir, idcollection_dict):
 
     
     
-def check_completeness(logger, model, progress, module, focus, eggnog, keggorg, idcollection_dict, summary_dict, outdir): 
+def check_completeness(logger, model, progress, module, focus, taxon, eggnog, keggorg, idcollection_dict, summary_dict, outdir): 
     # check KEGG annotations in the universe model to get '%' of completeness per pathway/module.
     
             
@@ -69,6 +93,9 @@ def check_completeness(logger, model, progress, module, focus, eggnog, keggorg,
     if keggorg != '-':  # keggorg has precedence
         kr_uni = parse_keggorg(keggorg, outdir, idcollection_dict)
         kr_uni_label = f"organism code '{keggorg}'"
+    elif taxon != '-':
+        kr_uni = parse_taxon(taxon, idcollection_dict)
+        kr_uni_label = f"taxon '{taxon}'"
     elif eggnog != '-':
         for eggfile in eggnog:
             eggset = parse_eggnog(model, eggfile, idcollection_dict)

gsrap/parsedb/parsedb.py

@@ -17,6 +17,7 @@ from ..commons import show_contributions
 from ..commons import adjust_biomass_precursors
 from ..commons import count_undrawn_rids
 from ..commons import format_expansion
+from ..commons import check_taxon
 from ..commons import download_keggorg
 from ..commons import initialize_model
 from ..commons import get_memote_results_dict
@@ -46,6 +47,7 @@ from .cycles import verify_egc_all
 def main(args, logger):
     
     
+    ###### FORMAT ARGS NOT REQUIRING RESOURCES
     # adjust out folder path                      
     while args.outdir.endswith('/'):              
         args.outdir = args.outdir[:-1]
@@ -77,17 +79,8 @@ def main(args, logger):
     if args.onlyauthor == '-': args.onlyauthor = None
     
     
-    # format the --eggnog param
-    args.eggnog = format_expansion(logger, args.eggnog)  # now 'args.eggnog' could still be '-'
-    
-    # get the kegg organism if requested
-    if args.keggorg != '-':
-        response = download_keggorg(logger, args.keggorg, args.outdir)
-        if response == 1: return 1
-    
-    
-    
     
+    ###### LOAD LOCAL RESOURCES
     # check and extract the required 'gsrap.maps' file
     if os.path.exists(f'{args.inmaps}') == False:
         logger.error(f"File 'gsrap.maps' not found at {args.inmaps}.")
@@ -108,9 +101,27 @@ def main(args, logger):
             kegg_compound_to_others = pickle.load(handle)
     with resources.path("gsrap.assets", f"kegg_reaction_to_others.pickle") as asset_path: 
         with open(asset_path, 'rb') as handle:
-            kegg_reaction_to_others = pickle.load(handle)   
+            kegg_reaction_to_others = pickle.load(handle) 
+            
+            
+            
+    ###### FORMAT/CHECK FOCUSING ARGS
+    # format the --eggnog param
+    args.eggnog = format_expansion(logger, args.eggnog)  # now 'args.eggnog' could still be '-'
+    
+    # check the --taxon param
+    if args.taxon != '-':
+        response = check_taxon(logger, args.taxon, idcollection_dict)
+        if response == 1: return 1
+    
+    # get the kegg organism if requested
+    if args.keggorg != '-':
+        response = download_keggorg(logger, args.keggorg, args.outdir)
+        if response == 1: return 1
+    
     
     
+    # DOWNLOAD ONLINE RESOURCES
     # get dbuni and dbexp:
     logger.info("Downloading gsrap database...")
     response = get_databases(logger)
@@ -166,7 +177,7 @@ def main(args, logger):
     
     ###### CHECKS 1
     # check universe completness
-    df_C = check_completeness(logger, universe, args.progress, args.module, args.focus, args.eggnog, args.keggorg, idcollection_dict, summary_dict, args.outdir)
+    df_C = check_completeness(logger, universe, args.progress, args.module, args.focus, args.taxon, args.eggnog, args.keggorg, idcollection_dict, summary_dict, args.outdir)
     if type(df_C)==int: return 1
 
 

{gsrap-0.8.2.dist-info → gsrap-0.8.3.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: gsrap
-Version: 0.8.2
+Version: 0.8.3
 Summary: 
 License: GNU General Public License v3.0
 Author: Gioele Lazzari

{gsrap-0.8.2.dist-info → gsrap-0.8.3.dist-info}/RECORD

@@ -1,5 +1,5 @@
-gsrap/.ipynb_checkpoints/__init__-checkpoint.py,sha256=o4xMVuraz_Qnf6ZOtxcYjNl-a7-zlRuZsP-GEw8lP4I,15240
-gsrap/__init__.py,sha256=o4xMVuraz_Qnf6ZOtxcYjNl-a7-zlRuZsP-GEw8lP4I,15240
+gsrap/.ipynb_checkpoints/__init__-checkpoint.py,sha256=Epw4X8B1O9kWnfC9v-X4jvDXE4X-G0XlRfkMdtNvMq0,15459
+gsrap/__init__.py,sha256=Epw4X8B1O9kWnfC9v-X4jvDXE4X-G0XlRfkMdtNvMq0,15459
 gsrap/assets/.ipynb_checkpoints/PM1-checkpoint.csv,sha256=0qjaMVG_t9aFxbHbxON6ecmEUnWPwN9nhmxc61QFeCU,8761
 gsrap/assets/.ipynb_checkpoints/PM2A-checkpoint.csv,sha256=rjYTdwe8lpRS552BYiUP3J71juG2ywVdR5Sux6fjZTY,8816
 gsrap/assets/.ipynb_checkpoints/PM3B-checkpoint.csv,sha256=42IGX_2O5bRYSiHoMuVKT-T-bzVj0cSRZBvGOrbnQMA,8130
@@ -14,9 +14,9 @@ gsrap/assets/kegg_reaction_to_others.pickle,sha256=AGW8CGN5hKeXZoYn3JRF4Xu832WyN
 gsrap/commons/.ipynb_checkpoints/__init__-checkpoint.py,sha256=9lrb0sBFSWEgV_e5FYzSgjTbam8b959rW_8VuxQHt1M,268
 gsrap/commons/.ipynb_checkpoints/biomass-checkpoint.py,sha256=4u7WBaUgo42tBoXDU1D0VUjICatb44e0jfswZrBeHYs,17987
 gsrap/commons/.ipynb_checkpoints/coeffs-checkpoint.py,sha256=qI3_GuqHkeA2KbK9pYdkqJaFwYemAVZJGLRR4QtHt6w,19182
-gsrap/commons/.ipynb_checkpoints/downloads-checkpoint.py,sha256=e-7ffMD4R07MWEgXyGcwjhScbWnG7A3L100YWbpNMk0,8461
+gsrap/commons/.ipynb_checkpoints/downloads-checkpoint.py,sha256=VnIFC6Y8QZ4oPyi04J_rtC2imzk5yOe0i5SpTP3MA6M,9713
 gsrap/commons/.ipynb_checkpoints/escherutils-checkpoint.py,sha256=_y0TgM0-Im0RT8W8z5rr4vlnGK55iRFds6DlDsjGD-8,1151
-gsrap/commons/.ipynb_checkpoints/excelhub-checkpoint.py,sha256=1KemKQ5oWziFlrT3__70-AKjXwdsqwEzVswEOSbNBDs,11433
+gsrap/commons/.ipynb_checkpoints/excelhub-checkpoint.py,sha256=_LtaWripY_D99f14Yk-3v9VIov2cUGzJ53AXgmNFpUk,11457
 gsrap/commons/.ipynb_checkpoints/figures-checkpoint.py,sha256=IRHSQXrCi4SQoISEfNB0rDhvUzbjcgsPi9zUSefsRto,4316
 gsrap/commons/.ipynb_checkpoints/fluxbal-checkpoint.py,sha256=jgC3-vI9Tbjvqohh2mJwFra4rl_pbUzHWrSa_QAxVO4,1262
 gsrap/commons/.ipynb_checkpoints/keggutils-checkpoint.py,sha256=M2nhHRiNH_xObHSxOIdt7ix59MrPdl9q3HNICC8X36M,4514
@@ -28,9 +28,9 @@ gsrap/commons/.ipynb_checkpoints/sbmlutils-checkpoint.py,sha256=gkY02qbGXrbYStn2
 gsrap/commons/__init__.py,sha256=9lrb0sBFSWEgV_e5FYzSgjTbam8b959rW_8VuxQHt1M,268
 gsrap/commons/biomass.py,sha256=4u7WBaUgo42tBoXDU1D0VUjICatb44e0jfswZrBeHYs,17987
 gsrap/commons/coeffs.py,sha256=qI3_GuqHkeA2KbK9pYdkqJaFwYemAVZJGLRR4QtHt6w,19182
-gsrap/commons/downloads.py,sha256=e-7ffMD4R07MWEgXyGcwjhScbWnG7A3L100YWbpNMk0,8461
+gsrap/commons/downloads.py,sha256=VnIFC6Y8QZ4oPyi04J_rtC2imzk5yOe0i5SpTP3MA6M,9713
 gsrap/commons/escherutils.py,sha256=_y0TgM0-Im0RT8W8z5rr4vlnGK55iRFds6DlDsjGD-8,1151
-gsrap/commons/excelhub.py,sha256=1KemKQ5oWziFlrT3__70-AKjXwdsqwEzVswEOSbNBDs,11433
+gsrap/commons/excelhub.py,sha256=_LtaWripY_D99f14Yk-3v9VIov2cUGzJ53AXgmNFpUk,11457
 gsrap/commons/figures.py,sha256=IRHSQXrCi4SQoISEfNB0rDhvUzbjcgsPi9zUSefsRto,4316
 gsrap/commons/fluxbal.py,sha256=jgC3-vI9Tbjvqohh2mJwFra4rl_pbUzHWrSa_QAxVO4,1262
 gsrap/commons/keggutils.py,sha256=M2nhHRiNH_xObHSxOIdt7ix59MrPdl9q3HNICC8X36M,4514
@@ -40,11 +40,11 @@ gsrap/commons/memoteutils.py,sha256=rulJFSVX4I-XGi4uHXloL0eGIkC5zhpuJYWJn9zCDbY,
 gsrap/commons/metrics.py,sha256=gvqF2c0e31m5qQWQ11JF4-eMqxtuONy_7lUiC7uaXX4,3291
 gsrap/commons/sbmlutils.py,sha256=gkY02qbGXrbYStn2F8J0KM0fmqati2Nbvi128EF7coI,365
 gsrap/getmaps/.ipynb_checkpoints/__init__-checkpoint.py,sha256=L4gLwk1vgnPlQuIP-FPnvy9uXnmGVQ4Rvjv3uyciDfk,92
-gsrap/getmaps/.ipynb_checkpoints/getmaps-checkpoint.py,sha256=BFAbBZFp-yaT35XO38YFHRqQr4mSAfRH6unOrPenQus,2703
-gsrap/getmaps/.ipynb_checkpoints/kdown-checkpoint.py,sha256=mv5wqP4DX9KhB1IIWx6cLrc-MjQdF5Nz-Gii8QSO_6I,24498
+gsrap/getmaps/.ipynb_checkpoints/getmaps-checkpoint.py,sha256=5Odf2PS1PmPZqQyKoLPYHYgNvYsp6IxrTy5HEQ5LNEM,2943
+gsrap/getmaps/.ipynb_checkpoints/kdown-checkpoint.py,sha256=OZGeGdnXMLzxvitwMTWKYZyL4dq3T8ZqbT4HVC3F3lI,27452
 gsrap/getmaps/__init__.py,sha256=L4gLwk1vgnPlQuIP-FPnvy9uXnmGVQ4Rvjv3uyciDfk,92
-gsrap/getmaps/getmaps.py,sha256=BFAbBZFp-yaT35XO38YFHRqQr4mSAfRH6unOrPenQus,2703
-gsrap/getmaps/kdown.py,sha256=mv5wqP4DX9KhB1IIWx6cLrc-MjQdF5Nz-Gii8QSO_6I,24498
+gsrap/getmaps/getmaps.py,sha256=5Odf2PS1PmPZqQyKoLPYHYgNvYsp6IxrTy5HEQ5LNEM,2943
+gsrap/getmaps/kdown.py,sha256=OZGeGdnXMLzxvitwMTWKYZyL4dq3T8ZqbT4HVC3F3lI,27452
 gsrap/mkmodel/.ipynb_checkpoints/__init__-checkpoint.py,sha256=PNze-26HMOwfdJ92KiXpr--VV1ftVfo3CAxBZgeokp8,92
 gsrap/mkmodel/.ipynb_checkpoints/biologcuration-checkpoint.py,sha256=Nn7z-js-mzzeO23kVM2L7sJ5PNle7AkCUeBcEAYjlFU,15378
 gsrap/mkmodel/.ipynb_checkpoints/gapfill-checkpoint.py,sha256=BPZw4sszlBhAYfHnV0pA7EpG0b2ePwS6kUfFt0Ww-ss,5159
@@ -61,19 +61,19 @@ gsrap/mkmodel/polishing.py,sha256=R1UdFPxN8N27Iu0jsYW2N_1BkWEbBHaMYW6NkCYZK_k,32
 gsrap/mkmodel/pruner.py,sha256=FAZid-0H6j66wR2dVKRAaMaDREVt1edflmZXbX7blXg,9836
 gsrap/parsedb/.ipynb_checkpoints/__init__-checkpoint.py,sha256=1k2K1gz4lIdXAwHEdJ0OhdkPu83woGv0Z4TpT1kGrTk,97
 gsrap/parsedb/.ipynb_checkpoints/annotation-checkpoint.py,sha256=Y02_zXJj_tS1GyBdfuLBy9YJjMgx3mjX6tqr1KhQ-9Q,4810
-gsrap/parsedb/.ipynb_checkpoints/completeness-checkpoint.py,sha256=yhFiEslK1qmMCk_GWZ7UZtX02FUqLU39UafG5886WsY,12016
+gsrap/parsedb/.ipynb_checkpoints/completeness-checkpoint.py,sha256=Op7VwmmwHmt1nhcl_0ISAejtLz-F9IkmnTrcJvO0BGc,12829
 gsrap/parsedb/.ipynb_checkpoints/cycles-checkpoint.py,sha256=HJ58LcHQseQ1eploysfXd5Y8Rip8n62qhje4pmL22VM,4761
 gsrap/parsedb/.ipynb_checkpoints/introduce-checkpoint.py,sha256=UuwGWGB2saG9VDMoboumeRBWhHOO68bs5_1r2RSkyVo,17145
 gsrap/parsedb/.ipynb_checkpoints/manual-checkpoint.py,sha256=qMKYshVftSGCRAjHC87E6n9-6kAiffFFCgHOUbqlyC0,3625
-gsrap/parsedb/.ipynb_checkpoints/parsedb-checkpoint.py,sha256=fTZlQvjV_wyRIxWJbW5TCVvydn0-CsYknsx9FYspdu8,8127
+gsrap/parsedb/.ipynb_checkpoints/parsedb-checkpoint.py,sha256=8mQgUTMOLpoeHK_X28s5jaW8adltKZ40nn_0uxmIXz8,8515
 gsrap/parsedb/.ipynb_checkpoints/repeating-checkpoint.py,sha256=WwPOzlZgsZWmJ-rvhFg21iOJ6gajgKFc2vCIHh6weBg,6103
 gsrap/parsedb/__init__.py,sha256=1k2K1gz4lIdXAwHEdJ0OhdkPu83woGv0Z4TpT1kGrTk,97
 gsrap/parsedb/annotation.py,sha256=Y02_zXJj_tS1GyBdfuLBy9YJjMgx3mjX6tqr1KhQ-9Q,4810
-gsrap/parsedb/completeness.py,sha256=yhFiEslK1qmMCk_GWZ7UZtX02FUqLU39UafG5886WsY,12016
+gsrap/parsedb/completeness.py,sha256=Op7VwmmwHmt1nhcl_0ISAejtLz-F9IkmnTrcJvO0BGc,12829
 gsrap/parsedb/cycles.py,sha256=HJ58LcHQseQ1eploysfXd5Y8Rip8n62qhje4pmL22VM,4761
 gsrap/parsedb/introduce.py,sha256=UuwGWGB2saG9VDMoboumeRBWhHOO68bs5_1r2RSkyVo,17145
 gsrap/parsedb/manual.py,sha256=qMKYshVftSGCRAjHC87E6n9-6kAiffFFCgHOUbqlyC0,3625
-gsrap/parsedb/parsedb.py,sha256=fTZlQvjV_wyRIxWJbW5TCVvydn0-CsYknsx9FYspdu8,8127
+gsrap/parsedb/parsedb.py,sha256=8mQgUTMOLpoeHK_X28s5jaW8adltKZ40nn_0uxmIXz8,8515
 gsrap/parsedb/repeating.py,sha256=WwPOzlZgsZWmJ-rvhFg21iOJ6gajgKFc2vCIHh6weBg,6103
 gsrap/runsims/.ipynb_checkpoints/__init__-checkpoint.py,sha256=6E6E1gWgH0V7ls4Omx4mxxC85gMJ_27YqhjugJzlZtY,97
 gsrap/runsims/.ipynb_checkpoints/biosynth-checkpoint.py,sha256=fUlHUo4CfB4rGX9Dth87B1p5E5sz7i6spR7ZoqDDGaI,2836
@@ -93,8 +93,8 @@ gsrap/runsims/precursors.py,sha256=1RNt_Rxs0L1lolDmYh4_CiZgiwHfU5B_AcomJO6vJ28,2
 gsrap/runsims/runsims.py,sha256=2FC5Gs8oSYyZTjHF3A7aXB_O6myVfcn3bCxQfLJlZTk,2842
 gsrap/runsims/simplegrowth.py,sha256=tCQHTMUqum1YwlBKRTNaQoag2co_yQlCaKmISOARAlE,2353
 gsrap/runsims/singleomission.py,sha256=jMuKAi0pINP8Jlrm-yI-tX7D110VzttR3YfTSnDRe4I,2847
-gsrap-0.8.2.dist-info/LICENSE.txt,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-gsrap-0.8.2.dist-info/METADATA,sha256=CEZpXejZ1b_KWKhxPbSXKGdW01ZU_nkKFqzfafDXXGE,898
-gsrap-0.8.2.dist-info/WHEEL,sha256=b4K_helf-jlQoXBBETfwnf4B04YC67LOev0jo4fX5m8,88
-gsrap-0.8.2.dist-info/entry_points.txt,sha256=S9MY0DjfnbKGlZbp5bV7W6dNFy3APoEV84u9x6MV1eI,36
-gsrap-0.8.2.dist-info/RECORD,,
+gsrap-0.8.3.dist-info/LICENSE.txt,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+gsrap-0.8.3.dist-info/METADATA,sha256=JADff6H-Y_SWY5PtR9qzEhmabesB6A5dLa-V0GTKqgc,898
+gsrap-0.8.3.dist-info/WHEEL,sha256=b4K_helf-jlQoXBBETfwnf4B04YC67LOev0jo4fX5m8,88
+gsrap-0.8.3.dist-info/entry_points.txt,sha256=S9MY0DjfnbKGlZbp5bV7W6dNFy3APoEV84u9x6MV1eI,36
+gsrap-0.8.3.dist-info/RECORD,,