gsrap 0.8.2__tar.gz → 0.9.0__tar.gz

{gsrap-0.8.2 → gsrap-0.9.0}/PKG-INFO

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

{gsrap-0.8.2 → gsrap-0.9.0}/pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "gsrap"
-version = "0.8.2"
+version = "0.9.0"
 description = ""
 authors = ["Gioele Lazzari"]
 license = "GNU General Public License v3.0"

{gsrap-0.8.2 → gsrap-0.9.0}/src/gsrap/.ipynb_checkpoints/__init__-checkpoint.py

@@ -75,12 +75,14 @@ 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.")
     parsedb_parser.add_argument("--onlyauthor", metavar='', type=str, default='-', help="Build the universe by parsing contents of the specified author ID only. Contents affected by --goodbefore are parsed anyway.")
     parsedb_parser.add_argument("--nofigs", action='store_true', help="Do not generate figures.")
     parsedb_parser.add_argument("-j", "--justparse", action='store_true', help="Just parse the database without performing extra activities (saves time during universe expansion).")
+    parsedb_parser.add_argument("-d", "--keepdisconn", action='store_true', help="Do not remove disconnected metabolites.")
     
     
     

{gsrap-0.8.2 → gsrap-0.9.0}/src/gsrap/__init__.py

@@ -75,12 +75,14 @@ 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.")
     parsedb_parser.add_argument("--onlyauthor", metavar='', type=str, default='-', help="Build the universe by parsing contents of the specified author ID only. Contents affected by --goodbefore are parsed anyway.")
     parsedb_parser.add_argument("--nofigs", action='store_true', help="Do not generate figures.")
     parsedb_parser.add_argument("-j", "--justparse", action='store_true', help="Just parse the database without performing extra activities (saves time during universe expansion).")
+    parsedb_parser.add_argument("-d", "--keepdisconn", action='store_true', help="Do not remove disconnected metabolites.")
     
     
     

{gsrap-0.8.2 → gsrap-0.9.0}/src/gsrap/commons/.ipynb_checkpoints/downloads-checkpoint.py

@@ -243,7 +243,99 @@ 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
+    
+    
+    """
+    sorted(list(df.query("kingdom == 'Bacteria'")['phylum'].unique()))
+    ['Acidobacteriota',
+     'Actinomycetota',
+     'Alphaproteobacteria',
+     'Aquificota',
+     'Armatimonadota',
+     'Atribacterota',
+     'Bacilli',
+     'Bacteria incertae sedis',
+     'Bacteroidota',
+     'Balneolota',
+     'Bdellovibrionota',
+     'Betaproteobacteria',
+     'Caldisericota',
+     'Calditrichota',
+     'Campylobacterota',
+     'Chlamydiota',
+     'Chlorobiota',
+     'Chloroflexota',
+     'Chrysiogenota',
+     'Cloacimonadota',
+     'Clostridia',
+     'Coprothermobacterota',
+     'Cyanobacteriota',
+     'Deferribacterota',
+     'Deinococcota',
+     'Deltaproteobacteria',
+     'Dictyoglomota',
+     'Elusimicrobiota',
+     'Enterobacteria',
+     'Fibrobacterota',
+     'Fidelibacterota',
+     'Fusobacteriota',
+     'Gemmatimonadota',
+     'Ignavibacteriota',
+     'Kiritimatiellota',
+     'Lentisphaerota',
+     'Melainabacteria',
+     'Mycoplasmatota',
+     'Myxococcota',
+     'Nitrospinota',
+     'Nitrospirota',
+     'Omnitrophota',
+     'Planctomycetota',
+     'Rhodothermota',
+     'Spirochaetota',
+     'Synergistota',
+     'Thermodesulfobacteriota',
+     'Thermodesulfobiota',
+     'Thermomicrobiota',
+     'Thermosulfidibacterota',
+     'Thermotogota',
+     'Verrucomicrobiota',
+     'Vulcanimicrobiota',
+     'other Bacillota',
+     'other Gammaproteobacteria',
+     'other Pseudomonadota',
+     'unclassified Bacteria']
+    """
+    
+    return 0

gsrap-0.9.0/src/gsrap/commons/.ipynb_checkpoints/escherutils-checkpoint.py

@@ -0,0 +1,108 @@
+import warnings
+import logging
+
+
+import cobra
+
+
+
+def print_json_tree(data, level=0, max_level=2):
+    # explore contents of a json object
+    
+    if level > max_level:
+        return
+    indent = '  ' * level
+    if isinstance(data, dict):
+        for key, value in data.items():
+            print(f"{indent}{key}")
+            print_tree(value, level + 1, max_level)
+    elif isinstance(data, list):
+        for i, item in enumerate(data):
+            print(f"{indent}[{i}]")
+            print_tree(item, level + 1, max_level)
+            
+            
+
+def count_undrawn_rids(logger, universe, lastmap, focus):
+    
+    
+    rids = set([r.id for r in universe.reactions])
+    
+    drawn_rids = set()
+    for key, value in lastmap['json'][1]['reactions'].items():
+        drawn_rids.add(value['bigg_id'])
+        
+        
+    remainings = rids - drawn_rids
+    filename = lastmap['filename']
+    logger.debug(f"Last universal map version detected: '{filename}'.")
+    if len(remainings) > 0:
+        logger.warning(f"Our universal map is {len(remainings)} reactions behind. Please draw!")
+        if focus == '-':
+            logger.warning(f"Drawing is eased when using '--focus'...")
+    else:
+        logger.info(f"Our universal map is {len(remainings)} reactions behind. Thank you ♥")
+        
+        
+        
+def count_undrawn_rids_focus(logger, universe, lastmap, focus, outdir):
+    
+    
+    # get modeled reads for this --focus:
+    rids = set()
+    try: gr = universe.groups.get_by_id(focus)
+    except: 
+        logger.warning(f"Group '{focus}' not found!")
+        return 
+    for r in gr.members:
+        rids.add(r.id)
+
+
+    # get rids on Escher: 
+    drawn_rids = set()
+    for key, value in lastmap['json'][1]['reactions'].items():
+        drawn_rids.add(value['bigg_id'])
+
+        
+    # get remaining rids for this map:
+    remainings = rids - drawn_rids
+    remainings_krs = set()
+    for rid in remainings: 
+        r = universe.reactions.get_by_id(rid)
+        krs = r.annotation['kegg.reaction']
+        for kr in krs: 
+            remainings_krs.add(kr)
+    
+    
+    if len(remainings) > 0:
+        if focus != 'transport':
+            logger.warning(f"Focusing on '{focus}', our universal map is {len(remainings)} reactions behind: {' '.join(list(remainings_krs))}.")
+        else:
+            logger.warning(f"Focusing on '{focus}', our universal map is {len(remainings)} reactions behind.")  # usually no kegg codes for tranport reactions
+        
+        
+        # subset the universe to ease the drawing: 
+        universe_focus = universe.copy()
+        to_remove = [r for r in universe_focus.reactions if r.id not in rids]
+        
+        
+        # trick to avoid the WARNING "cobra/core/group.py:147: UserWarning: need to pass in a list" 
+        # triggered when trying to remove reactions that are included in groups. 
+        with warnings.catch_warnings():  # temporarily suppress warnings for this block
+            warnings.simplefilter("ignore")  # ignore all warnings
+            cobra_logger = logging.getLogger("cobra.util.solver")
+            old_level = cobra_logger.level
+            cobra_logger.setLevel(logging.ERROR)   
+           
+            universe_focus.remove_reactions(to_remove,remove_orphans=True)
+
+            # restore original behaviour: 
+            cobra_logger.setLevel(old_level)  
+            
+            
+        # save the subset for drawing in Escher!
+        logger.info(f"Writing '{outdir}/{focus}.json' to ease your drawing workflow...")
+        cobra.io.save_json_model(universe_focus, f'{outdir}/{focus}.json')
+    else:
+        logger.info(f"Focusing on '{focus}', our universal map is {len(remainings)} reactions behind. Thank you ♥")
+        

{gsrap-0.8.2 → gsrap-0.9.0}/src/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-0.8.2 → gsrap-0.9.0}/src/gsrap/commons/downloads.py

@@ -243,7 +243,99 @@ 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
+    
+    
+    """
+    sorted(list(df.query("kingdom == 'Bacteria'")['phylum'].unique()))
+    ['Acidobacteriota',
+     'Actinomycetota',
+     'Alphaproteobacteria',
+     'Aquificota',
+     'Armatimonadota',
+     'Atribacterota',
+     'Bacilli',
+     'Bacteria incertae sedis',
+     'Bacteroidota',
+     'Balneolota',
+     'Bdellovibrionota',
+     'Betaproteobacteria',
+     'Caldisericota',
+     'Calditrichota',
+     'Campylobacterota',
+     'Chlamydiota',
+     'Chlorobiota',
+     'Chloroflexota',
+     'Chrysiogenota',
+     'Cloacimonadota',
+     'Clostridia',
+     'Coprothermobacterota',
+     'Cyanobacteriota',
+     'Deferribacterota',
+     'Deinococcota',
+     'Deltaproteobacteria',
+     'Dictyoglomota',
+     'Elusimicrobiota',
+     'Enterobacteria',
+     'Fibrobacterota',
+     'Fidelibacterota',
+     'Fusobacteriota',
+     'Gemmatimonadota',
+     'Ignavibacteriota',
+     'Kiritimatiellota',
+     'Lentisphaerota',
+     'Melainabacteria',
+     'Mycoplasmatota',
+     'Myxococcota',
+     'Nitrospinota',
+     'Nitrospirota',
+     'Omnitrophota',
+     'Planctomycetota',
+     'Rhodothermota',
+     'Spirochaetota',
+     'Synergistota',
+     'Thermodesulfobacteriota',
+     'Thermodesulfobiota',
+     'Thermomicrobiota',
+     'Thermosulfidibacterota',
+     'Thermotogota',
+     'Verrucomicrobiota',
+     'Vulcanimicrobiota',
+     'other Bacillota',
+     'other Gammaproteobacteria',
+     'other Pseudomonadota',
+     'unclassified Bacteria']
+    """
+    
+    return 0

gsrap-0.9.0/src/gsrap/commons/escherutils.py

@@ -0,0 +1,108 @@
+import warnings
+import logging
+
+
+import cobra
+
+
+
+def print_json_tree(data, level=0, max_level=2):
+    # explore contents of a json object
+    
+    if level > max_level:
+        return
+    indent = '  ' * level
+    if isinstance(data, dict):
+        for key, value in data.items():
+            print(f"{indent}{key}")
+            print_tree(value, level + 1, max_level)
+    elif isinstance(data, list):
+        for i, item in enumerate(data):
+            print(f"{indent}[{i}]")
+            print_tree(item, level + 1, max_level)
+            
+            
+
+def count_undrawn_rids(logger, universe, lastmap, focus):
+    
+    
+    rids = set([r.id for r in universe.reactions])
+    
+    drawn_rids = set()
+    for key, value in lastmap['json'][1]['reactions'].items():
+        drawn_rids.add(value['bigg_id'])
+        
+        
+    remainings = rids - drawn_rids
+    filename = lastmap['filename']
+    logger.debug(f"Last universal map version detected: '{filename}'.")
+    if len(remainings) > 0:
+        logger.warning(f"Our universal map is {len(remainings)} reactions behind. Please draw!")
+        if focus == '-':
+            logger.warning(f"Drawing is eased when using '--focus'...")
+    else:
+        logger.info(f"Our universal map is {len(remainings)} reactions behind. Thank you ♥")
+        
+        
+        
+def count_undrawn_rids_focus(logger, universe, lastmap, focus, outdir):
+    
+    
+    # get modeled reads for this --focus:
+    rids = set()
+    try: gr = universe.groups.get_by_id(focus)
+    except: 
+        logger.warning(f"Group '{focus}' not found!")
+        return 
+    for r in gr.members:
+        rids.add(r.id)
+
+
+    # get rids on Escher: 
+    drawn_rids = set()
+    for key, value in lastmap['json'][1]['reactions'].items():
+        drawn_rids.add(value['bigg_id'])
+
+        
+    # get remaining rids for this map:
+    remainings = rids - drawn_rids
+    remainings_krs = set()
+    for rid in remainings: 
+        r = universe.reactions.get_by_id(rid)
+        krs = r.annotation['kegg.reaction']
+        for kr in krs: 
+            remainings_krs.add(kr)
+    
+    
+    if len(remainings) > 0:
+        if focus != 'transport':
+            logger.warning(f"Focusing on '{focus}', our universal map is {len(remainings)} reactions behind: {' '.join(list(remainings_krs))}.")
+        else:
+            logger.warning(f"Focusing on '{focus}', our universal map is {len(remainings)} reactions behind.")  # usually no kegg codes for tranport reactions
+        
+        
+        # subset the universe to ease the drawing: 
+        universe_focus = universe.copy()
+        to_remove = [r for r in universe_focus.reactions if r.id not in rids]
+        
+        
+        # trick to avoid the WARNING "cobra/core/group.py:147: UserWarning: need to pass in a list" 
+        # triggered when trying to remove reactions that are included in groups. 
+        with warnings.catch_warnings():  # temporarily suppress warnings for this block
+            warnings.simplefilter("ignore")  # ignore all warnings
+            cobra_logger = logging.getLogger("cobra.util.solver")
+            old_level = cobra_logger.level
+            cobra_logger.setLevel(logging.ERROR)   
+           
+            universe_focus.remove_reactions(to_remove,remove_orphans=True)
+
+            # restore original behaviour: 
+            cobra_logger.setLevel(old_level)  
+            
+            
+        # save the subset for drawing in Escher!
+        logger.info(f"Writing '{outdir}/{focus}.json' to ease your drawing workflow...")
+        cobra.io.save_json_model(universe_focus, f'{outdir}/{focus}.json')
+    else:
+        logger.info(f"Focusing on '{focus}', our universal map is {len(remainings)} reactions behind. Thank you ♥")
+        

{gsrap-0.8.2 → gsrap-0.9.0}/src/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-0.8.2 → gsrap-0.9.0}/src/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