Ucsc Track Exonstarts/Exonends To Bed

Example from UCSC track:

    bin    name    chrom    strand    txStart    txEnd    cdsStart    cdsEnd    exonCount    exonStarts    exonEnds
    73    TCONS_l2_00002368    chr1    -    89294    237877    237877    237877    2    89294,236614,    90404,237877,
    585    TCONS_l2_00002369    chr1    -    89550    91105    91105    91105    2    89550,90286,    90050,91105,

How to combine exonStarts with corresponding exonEnds to get only exon's bed file?

   chr1    89294    90404
   chr1    236614    237877
   chr1    89550    90050
   chr1    90286    91105

bed ucsc

• 2,928 views • link

updated 13.3 years ago by Duarte Molha

• written 13.3 years ago by PoGibas

2 answers

Install bedtools...

instead of trying to do that from the track, download the track as a bed12 file and use the command bed12ToBed6 to get all exon coordinates:

bed12ToBed6 -i UCSC_track.bed > UCSC_discrete_features.bed

bedtools bed12tobed6 [OPTIONS] -i UCSC_track.bed > UCSC_discrete_features.bed

If this is not possible a relatively simple awk script should do the trick:

awk 'BEGIN{OFS="\t"}{n_exons = split($10,e_starts,",");split($11,e_ends,",");for(i=1;i<n_exons;i++){print $3,e_starts[i],e_ends[i]}}' t.track > exons.bed

• 8 views • link

written 13.3 years ago by Duarte Molha

I was using bedtools for some time and didn't know it had such a cool feature.

• 0 views • link

written 13.1 years ago by PoGibas

I specifically created a biostars account so that I could upvote this reply. Thank you very much!

• 0 views • link

written 10.2 years ago by rmbannen •