annotate mytools/bedClean.py @ 0:39217fa39ff2

Uploaded
author xuebing
date Tue, 13 Mar 2012 23:34:52 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
1 import sys
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
2
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
3 def readChrSize(filename):
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
4 f = open(filename)
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
5 chrSize = {}
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
6 for line in f:
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
7 chrom,size = line.strip().split()
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
8 chrSize[chrom]=int(size)
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
9 f.close()
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
10 return chrSize
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
11
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
12 def cleanFile(filename,chrSize,outfile):
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
13 f = open(filename)
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
14 out = open(outfile,'w')
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
15 i = 0
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
16 for line in f:
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
17 i = i + 1
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
18 flds = line.strip().split('\t')
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
19 if len(flds) < 3:
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
20 print 'line',i,'incomplete line:\n',line
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
21 elif chrSize.has_key(flds[0]):
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
22 if int(flds[1]) > int(flds[2]):
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
23 tmp = flds[1]
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
24 flds[1] = flds[2]
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
25 flds[2] = tmp
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
26 if int( flds[1]) < 0 or int(flds[2]) <0:
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
27 print 'line',i,'negative coordinates:\n',line
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
28 elif int(flds[2]) > chrSize[flds[0]]:
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
29 print 'line',i,'end larger than chr size:\n',line
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
30 else:
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
31 out.write('\t'.join(flds)+'\n')
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
32 else:
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
33 print 'line',i,'chromosome',flds[0],'not found!\n',line
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
34 f.close()
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
35 out.close()
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
36
39217fa39ff2 Uploaded
xuebing
parents:
diff changeset
37 cleanFile(sys.argv[1],readChrSize(sys.argv[2]),sys.argv[3])