Create a list of targets for the content of the community metabolic potential about metage2metabo HOT 3 CLOSED

Hi @cfrioux, Hi @ArnaudBelcour,

This would be great! When will this option become available?

Best regards,

Noor Alsmadi

from metage2metabo.

Creating a list of targets for all the reachable compounds of the community (community scope) is not implemented yet.
I do not have time to properly include it in m2m these days, but, I've got a snippet to do it, if you're interested.

import sys
import json
from metage2metabo.sbml_management import create_species_sbml

# inputs
jsonfile = sys.argv[1]
output_file_sbml = sys.argv[2]

# read scope file and retrieve targets
with open(jsonfile, "r") as f:
    raw = json.load(f)

targets = raw["com_scope"]

# create a SBML with the targets
create_species_sbml(targets, output_file_sbml)

You write the above script in a file (let's say scope_to_targets.py) and call it this way:

python scope_to_targets.py community_analysis/comm_scopes.json community_analysis/comm_scopes.sbml

from metage2metabo.

Added in #61.

from metage2metabo.