annotate igblastmerge.py @ 6:f0bf6e978f77 draft

Uploaded
author davidvanzessen
date Tue, 10 Dec 2013 08:16:24 -0500
parents 39ff573fed64
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
1 import sys
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
2 # error
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
3 def stop_err( msg ):
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
4 sys.stderr.write( "%s\n" % msg )
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
5 sys.exit()
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
6
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
7 # main
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
8 def main():
6
f0bf6e978f77 Uploaded
davidvanzessen
parents: 3
diff changeset
9 args = sys.argv[1:-2]
3
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
10
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
11 try:
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
12 o = open(sys.argv[-1], 'w')
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
13 i = open(args[1], 'r')
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
14 separator = "\t"
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
15 newline = "\n"
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
16 line = i.readline()
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
17 #write the header
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
18 o.write(line[:line.rfind(newline)] + separator + "Sample" + separator + "Replicate" + newline)
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
19 i.close()
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
20
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
21 current = 1
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
22 sampleID = args[0]
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
23 count = 1
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
24
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
25 while True:
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
26 print str(o)
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
27 f = open(args[current], 'r')
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
28 line = f.readline()
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
29 line = f.readline()
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
30 while line:
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
31 o.write(line[:line.rfind(newline)] + separator + sampleID + separator + str(count) + newline)
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
32 line = f.readline()
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
33 f.close()
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
34
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
35 if current >= (len(args) - 1):
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
36 break
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
37 if args[current + 1].find("/") is -1:
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
38 sampleID = args[current + 1]
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
39 current += 1
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
40 count = 1
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
41 else:
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
42 count += 1
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
43 current += 1
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
44 o.close()
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
45
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
46 except Exception, ex:
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
47 stop_err('Error running new_column.py\n' + str(ex))
39ff573fed64 Uploaded
davidvanzessen
parents: 0
diff changeset
48 sys.exit(0)
0
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
49
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
50 if __name__ == "__main__":
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
51 print sys.argv
1f6de7db44b6 Uploaded
davidvanzessen
parents:
diff changeset
52 main()