from FigKernelPackages import FIG2
import sys

fig = FIG2.FIG()

def die(msg):
	print msg
	sys.exit(0)

if len(sys.argv) != 3:
	die("usage: connect_subsystem_to_genes_the_right_way Subsystem Genome")

subsystem = sys.argv[1]
genome  = sys.argv[2]

print subsystem
roles = fig.subsystem_to_roles(subsystem)

for role in roles:
    print role
    pegs = fig.pegs_in_subsystem_cell(subsystem, genome, role)
    for peg in pegs: 
	func = fig.function_of(peg)
	print "    %s\t%s" % (peg, func[0][1])