Mercurial > repos > xuebing > sharplab_seq_motif
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()