Mercurial > repos > artbio > snvtocnv

annotate sanitize_snv.py @ 3:15b26e44b18d draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8-dirty"
author artbio
date Fri, 19 Feb 2021 15:22:05 +0000
parents 0e54da14f831
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
1 import sys
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
2
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
3 handle = open(sys.argv[1], 'r')
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
4 out = open(sys.argv[2], 'w')
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
5 for line in handle:
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
6 if line[0] == '#':
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
7 out.write(line)
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
8 continue
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
9 linelist = line.split('\t')
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
10 refcol = linelist[0].split('chr')
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
11 infocol = linelist[7].split('INDEL')
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
12 if len(infocol) > 1:
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
13 continue
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
14 if len(refcol) > 1:
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
15 refcol = refcol[1]
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
16 else:
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
17 refcol = refcol[0]
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
18 if refcol not in ['1', '2', '3', '4', '5', '6', '7', '8', '9', '10',
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
19 '11', '12', '13', '14', '15', '16', '17', '18', '19',
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
20 '20', '21', '22', 'X', 'Y']:
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
21 continue
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
22 else:
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
23 linelist[0] = refcol
0e54da14f831 "planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/snvtocnv commit fd049d0796cbfaa6c41a7d80e84d3a734b30acc8"
artbio
parents:
diff changeset
24 out.write('\t'.join(linelist))