# HG changeset patch # User xuebing # Date 1333249093 14400 # Node ID 7c1a351ad7334c40fc0feedf5a0049efc562113a # Parent 4ac24aaed8d66e4b030880dd105515aa93104b10 Uploaded diff -r 4ac24aaed8d6 -r 7c1a351ad733 ensembl_to_ucsc.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/ensembl_to_ucsc.py Sat Mar 31 22:58:13 2012 -0400 @@ -0,0 +1,26 @@ +''' +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() +