[Genome] trying to extract DNA sequences for introns only

Kayla Smith kayla at soe.ucsc.edu
Thu May 15 15:09:13 PDT 2008


Hello Sheila,

There are two things going on here which will help you understand your 
results:

1.  The Known Genes track will have items in it that overlap each other. 
  That is, a region that is an exon for one item, may be an intron for 
another item.

2.  When you select "Encode Regions" in the Table Browser, your results 
will include the introns from any Known Gene that overlapped an encode 
region --even if that particular intron does not itself overlap an 
encode region.  In order to get around this problem, you can use the 
Table Browser to do a "basewise" intersection with the encodeRegions table.

Here is a link to a session that I made with 2 custom tracks.  The 
second one, "KGintronsencode_2" is probably what you are looking for:
http://genome.ucsc.edu/cgi-bin/hgTracks?hgS_doOtherUser=submit&hgS_otherUserName=Kayla&hgS_otherUserSessionName=kayla_intron_encode_for_Sheila

If you want to get sequence, go to the Table Browser, select "group: 
Custom Tracks", track: "KGintronsencode_2", output format = "sequence"

I hope this information is helpful to you.  Please don't hesitate to 
contact us again if you require further assistance.

Kayla Smith
UCSC Genome Bioinformatics Group




Sheila Reynolds wrote:
> I'm trying to extract the DNA sequences for introns only from the Genome
> Browser.  Through the Tables option I tried to ask for the introns only from
> the knownGene track for the ENCODE regions, but if I look at one of the
> regions that is returned (for example chr6:41874646-41971039), it has some
> exons in it (if I am interpreting the genome browser depiction correctly).
> 
> I'd appreciate any tips, I'm rather new at this.
> 
> thanks,
> 
> Sheila Reynolds
> Dept of Electrical Engineering
> University of Washington
> _______________________________________________
> Genome maillist  -  Genome at soe.ucsc.edu
> http://www.soe.ucsc.edu/mailman/listinfo/genome



More information about the Genome mailing list