Problem In Obtaining Gene Names From Ucsc
2
0
Entering edit mode
8.5 years ago

Hi I want to obtain all the exons for mouse (mm10) with chromosome start_pos and end_pos followed by gene name. In UCSC table browser I have selected group: genes and gene prediction track; track: refseq genes; table refseq. In the output format tab i just selected custom track. I was getting the output as follows after i select the Exons plus and hit the get custom track in file button

chr1    134199221    134203590    NM_001039510_exon_0_0_chr1_134199222_r    0    -
chr1    134234014    134234412    NM_001039510_exon_1_0_chr1_134234015_r    0    -
chr1    134235227    134235431    NM_001039510_exon_2_0_chr1_134235228_r    0    -

However, i wanted it as

   chr1        134199221    134203590    Adora1_exon_0_0_chr1_134199222_r    0    -
   chr1       134234014    134234412    Adora1_exon_1_0_chr1_134234015_r    0    -
    chr1    134235227    134235431   Adora1_exon_2_0_chr1_134235228_r    0    -
ucsc genes • 1.6k views
ADD COMMENT
0
Entering edit mode

you need to select the name2 in the output (while selecting custom columns)!!

ADD REPLY
0
Entering edit mode

Unless you know differently, I don't think it's possible to get exon coordinates in BED format with gene names instead of transcript IDs

ADD REPLY
0
Entering edit mode

I don't get custom column selection option here in this case

ADD REPLY
0
Entering edit mode
8.5 years ago
always_learning ★ 1.1k

Can you please let us know from Where you got this dat ?

chr1       134199221    134203590    Adora1_exon_0_0_chr1_134199222_r    0    -
chr1       134234014    134234412    Adora1_exon_1_0_chr1_134234015_r    0    -
chr1       134235227    134235431    Adora1_exon_2_0_chr1_134235228_r    0    -

I understand this is different name of same gene !!

ADD COMMENT
0
Entering edit mode
8.5 years ago
Mitch Bekritsky ★ 1.3k

When I've wanted to substitute gene names for transcript IDs in a BED file, I used a dictionary I got from selecting "selected fields from primary and related tables" under output format in the table browser, then selecting only the "name" and "name2" fields on the page linked from there, which will provide a list of transcript IDs and their corresponding gene names. You may have to write a quick script to do the substitution in perl, R, etc.

ADD COMMENT
0
Entering edit mode

Yeah i was doing the same thing but just wanted to know whether there is a method without code within the ucsc browser!

ADD REPLY
0
Entering edit mode

Not sure that there is one...that's two of us who'd like to know now!

ADD REPLY

Login before adding your answer.

Traffic: 3042 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6