BedInverseExons

From genomewiki
Revision as of 12:09, 8 October 2007 by Max (talk | contribs)
Jump to navigationJump to search

When you export exons as bed with the table browser, they are numbered 0,1,2,3... by default, but always 5' to 3'. This is uncommon for biologists, they usually call the most upstream one exon 1. This scripts inverses the numbers.

#!/usr/bin/python 

# input: bed file with exons as exported from ucsc
# output: bed file with exons numbers reversed if strand="-"

from sys import *

f = open(argv[1], "r")
maxExon = {}
for l in f:
    fs = l.split()
    fs2 = fs[3].split("_")
    gene,exon = fs2[0],fs2[1] # might change this to fs2[0], fs[2]

    maxExon[gene] = int(exon)

f = open(argv[1], "r")

for l in f:
    fs = l.split()
    fs2 = fs[3].split("_")
    gene = fs2[0]
    exon = fs2[2]
    if fs[5]=="-":
        name = gene + "_exon_" + str(maxExon[gene] - int(exon))
    else:
        name = gene + "_exon_" + exon
    print "\t".join([fs[0], fs[1], fs[2],name,fs[4],fs[5]])