view mytools/convertEnsembl.py @ 0:39217fa39ff2

Uploaded
author xuebing
date Tue, 13 Mar 2012 23:34:52 -0400
parents
children
line wrap: on
line source

'''
convert ensembl bed to ucsc
add chr to chromosome
1 = +
-1 = -
'''

import sys
f = open(sys.argv[1])
out = open(sys.argv[2],'w')
skip = int(sys.argv[3])

for i in range(skip):
    f.readline()

for line in f:
    flds = line.strip().split('\t')
    flds[0] = 'chr'+flds[0]
    if flds[5] == '1':
        flds[5] = '+'
    else:
        flds[5] = '-'
    out.write('\t'.join(flds)+'\n')
f.close()
out.close()