# NOT RUN {
##-------------------------------------##
## A) Example query data
##-------------------------------------##
## 1 2 3 4 5 6 7 ##
qry <- data.frame(id = 1:7, seqid = "1",
lstart = c(10100L, 11800L, 12220L, 12220L, 12220L, 32000L, 40000L),
lend = c(10100L, 12000L, 12225L, 12227L, 12227L, 32100L, 40100L),
rstart = c(10200L, 12200L, 12057L, 12613L, 12650L, 32200L, 40200L),
rend = c(10300L, 12250L, 12179L, 12620L, 12700L, 32300L, 40300L))
## 1 2 3 4 5 6 7 ##
##-------------------------------------##
## B) Example reference data
##-------------------------------------##
# B.1) Ensembl genome:
ensfile <- system.file("extdata", "hs.ensembl.62.small.RData", package="refGenome")
ens <- loadGenome(ensfile)
gp <- getGenePositions(ens)
# B.2) Ensembl junctions:
junc <- getSpliceTable(ens)
##-------------------------------------##
## C) Do overlap
##-------------------------------------##
res <- overlapJuncs(qry, junc)
# }
Run the code above in your browser using DataLab