# HG changeset patch
# User francesco_lapi
# Date 1759332861 0
# Node ID 5956dcf9427766b50aa3481b4d60e89d0c6376f0
# Parent  ca98c149ec619552b03ce937df7066ef05aa8216
Uploaded

diff -r ca98c149ec61 -r 5956dcf94277 COBRAxy/ras_generator.py
--- a/COBRAxy/ras_generator.py	Wed Oct 01 14:21:26 2025 +0000
+++ b/COBRAxy/ras_generator.py	Wed Oct 01 15:34:21 2025 +0000
@@ -191,6 +191,7 @@
     for rule, reaction_ids in dict_rule_reactions.items():
         if len(rule) != 0:
             # there is one gene at least in the formula
+            rule = rule.replace("-","_")
             rule_split = rule.split()
             rule_split_elements = list(set(rule_split))                                # genes in formula
             
@@ -220,8 +221,6 @@
                             else:
                                 # complex expression (e.g. A or (B and C))
                                 data = count_df_filtered.loc[genes_in_count_matrix]  # dataframe of genes in the GPRs
-
-                                rule = rule.replace("-", "_")
                                 tree = ast.parse(rule, mode="eval").body
                                 values_by_cell = [dict(zip(data.index, data[col].values)) for col in data.columns]
                                 for j, values in enumerate(values_by_cell):