[Bio] / FigTutorial / connect_map_to_genes.py Repository:
ViewVC logotype

Annotation of /FigTutorial/connect_map_to_genes.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (view) (download) (as text)

1 : disz 1.1 from FigKernelPackages import FIG2
2 :     import sys
3 :    
4 :     fig = FIG2.FIG()
5 :    
6 :     def die(msg):
7 :     print msg
8 :     sys.exit(0)
9 :    
10 :     if len(sys.argv) != 3:
11 :     die("usage: connect_map_to_genes Map Genome")
12 :    
13 :     map = sys.argv[1]
14 :     genome = sys.argv[2]
15 :     name = fig.map_name(map);
16 :    
17 :     print "%s\t%s" % (map, name)
18 :     roles = fig.map_to_ecs(map)
19 :     for role in roles:
20 :     print role
21 :     pegs = fig.seqs_with_role(role,"master",genome)
22 :     for peg in pegs:
23 :     func = fig.function_of(peg);
24 :     print " %s\t%s" % (peg, func[0][1])

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3