getting snp and restriction Enzyme site from ucsc genome browswer
1
0
Entering edit mode
prabin.dm • 0
@prabindm-9986
Last seen 20 months ago
United States

I am using rtracklayer package to access rs ID and restriction enzyme sites in a genomic interval so that I can make genotyping decisions.

library(rtracklayer)
session <- browserSession()
genome(session) <- "mm9"
trackNames(session) ## list the track names
# for Restr Enzymes use "cutters"

query <- ucscTableQuery(session,"cutters" ,
                        GRangesForUCSCGenome("mm9", "chr7",
                                             IRanges(149627713, 149630315 )))

This yeilds a mesage saying Error in normArgTrack(track, trackids) : Unknown track: cutters

I read that "Note that not all tables are available as tracks.". So I am guessing "cutters" are available as a table but not as a tack. So how to I know which track to query so that I can get cutter table.

rtracklayer • 961 views
ADD COMMENT
0
Entering edit mode
@michael-lawrence-3846
Last seen 2.9 years ago
United States

The table name is "cutters", so you would need to pass that as the table= argument, leaving track=NULL so that it searches all tracks. But in this case UCSC does not provide that table anyway. You can get it from REBASE.

ADD COMMENT

Login before adding your answer.

Traffic: 799 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