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

Annotation of /FigTutorial/contigs_for_genome.py

Parent Directory Parent Directory | Revision Log Revision Log


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

1 : disz 1.1 from FigKernelPackages import FIG
2 :     import sys
3 :    
4 :     fig = FIG.FIG()
5 :    
6 :     def die(msg):
7 :     print msg
8 :     sys.exit(0)
9 :    
10 :     if len(sys.argv) != 2:
11 :     die("usage: contigs_for_genome Genome")
12 :    
13 :     genome = sys.argv[1]
14 :    
15 :     contigs = fig.all_contigs(genome)
16 :    
17 :     for contig in contigs:
18 :     len = fig.contig_ln(genome,contig)
19 :     if len >= 50:
20 :     seq = fig.dna_seq(genome,contig+"_1_50")
21 :     print contig+"\t"+len+"\t"+seq
22 :    

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3