Revamp single swap statistics

Author: Paulocracy

cosa_get_num_active_nad_nadp_reactions.py

@@ -31,7 +31,15 @@
             flux = results[growth_rate]["values"][reaction_id]
             if flux > 1e-6:
                 if reaction_id not in active_reactions:
-                    active_reactions.append(reaction_id)
+                    active_reactions.append(
+                        reaction_id
+                            .replace("_ORIGINAL_NADP_TCOSA", "")
+                            .replace("_ORIGINAL_NAD_TCOSA", "")
+                            .replace("_VARIANT_NAD_TCOSA", "")
+                            .replace("_VARIANT_NADP_TCOSA", "")
+                            .replace("_FWD", "")
+                            .replace("_REV", "")
+                    )
 print(active_reactions)
 active_reactions = list(set(active_reactions))
 print(f"Num  active reactions:", len(active_reactions))

cosa_swap_analysis_statistics.py

@@ -114,8 +114,8 @@
             id_optmdf_at_one = list(set(make_base(id_optmdf_at_one)))
             id_optsubmdf_at_one = list(set(make_base(id_optsubmdf_at_one)))
 
-            print(f"OptMDF at all (n={len(id_optmdf_at_all)}): {id_optmdf_at_all}")
-            print(f"OptSubMDF at all (n={len(id_optsubmdf_at_all)}): {id_optsubmdf_at_all}")
+            # print(f"OptMDF at all (n={len(id_optmdf_at_all)}): {id_optmdf_at_all}")
+            # print(f"OptSubMDF at all (n={len(id_optsubmdf_at_all)}): {id_optsubmdf_at_all}")
             print(f"OptMDF at >= 1 (n={len(id_optmdf_at_one)}): {id_optmdf_at_one}")
             print(f"OptSubMDF at >= 1 (n={len(id_optsubmdf_at_one)}): {id_optsubmdf_at_one}")