-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathrepeatMask.py
More file actions
36 lines (26 loc) · 788 Bytes
/
repeatMask.py
File metadata and controls
36 lines (26 loc) · 788 Bytes
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
import sys
for line in sys.stdin:
out=[]
line=line.rstrip()
ln=line.split("\t")
out.append(ln[0])
out.append(ln[1])
out.append(ln[2])
#rlen ref_span mapped_span
rlen = ln[3].split(":")[-1]
ref_span = int(ln[2]) - int(ln[1])
out.append( str(rlen))
out.append( str(ref_span) )
spann = int(ln[6]) - int(ln[5])
out.append( str(spann) )
if ln[15] != "-1":
#print(ln[11:13])
max_s = int(ln[15]) #max(int(ln[11]),int(ln[1]))
min_e = int(ln[16]) #min( int(ln[12]),int(ln[2]) )
else:
max_s = 0
min_e = 0
inter = min_e - max_s
out.append( str(inter))
sys.stdout.write( "\t".join(out) +"\n" )
## chr start end rlen ref_span mapped_span