###################################################
### chunk number 1: 
###################################################
#line 39 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
library(org.Bt.eg.db)
genes <- c('280705', '280706', '100327208')
symbols <- org.Bt.egSYMBOL[genes]
toTable(symbols)


###################################################
### chunk number 2: 
###################################################
#line 56 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
symbols <- c('BCSE','TF','TG')
entrez <- org.Bt.egSYMBOL2EG[symbols]
entrez <- toTable(entrez)[,'gene_id']
entrez
refseq <- org.Bt.egREFSEQ[entrez]
toTable(refseq)


###################################################
### chunk number 3: 
###################################################
#line 72 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
library(AnnotationFuncs)
translate(c(280705, 280706, 100327208, 123), org.Bt.egSYMBOL)
translate(c('BCSE','TF','TG'), from=org.Bt.egSYMBOL2EG, to=org.Bt.egREFSEQ)


###################################################
### chunk number 4: 
###################################################
#line 102 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
library(AnnotationFuncs)
library(org.Bt.eg.db)
genes <- c(280705, 280706, 100327208)
translate(genes, org.Bt.egGENENAME)
translate(genes, org.Bt.egCHR)
translate(genes, org.Bt.egENSEMBL)
translate(genes, org.Bt.egREFSEQ, remove.missing=F)


###################################################
### chunk number 5: 
###################################################
#line 116 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
symbols <- c("SERPINA1","KERA","CD5")
translate(symbols, org.Bt.egSYMBOL2EG)
translate(symbols, revmap(org.Bt.egSYMBOL))


###################################################
### chunk number 6: 
###################################################
#line 125 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
translate(symbols, from=org.Bt.egSYMBOL2EG, to=org.Bt.egGENENAME)
# As a curiosity, if you specify another organism, you are warned:
library(org.Hs.eg.db)
translate(symbols, from=org.Bt.egSYMBOL2EG, to=org.Hs.egGENENAME)
warnings()


###################################################
### chunk number 7: 
###################################################
#line 137 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
symbols <- c("SERPINA1","KERA","CD5")
translate(symbols, org.Bt.egSYMBOL2EG, org.Bt.egREFSEQ)
translate(symbols, org.Bt.egSYMBOL2EG, org.Bt.egREFSEQ, reduce='first')
translate(symbols, org.Bt.egSYMBOL2EG, org.Bt.egREFSEQ, reduce='last')


###################################################
### chunk number 8: 
###################################################
#line 146 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
translate(symbols, org.Bt.egSYMBOL2EG, org.Bt.egREFSEQ, return.list=F)


###################################################
### chunk number 9: 
###################################################
#line 160 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
symbols <- c("SERPINA1","KERA","CD5")
refseq <- translate(symbols, from=org.Bt.egSYMBOL2EG, to=org.Bt.egREFSEQ)
mRNA <- pickRefSeq(refseq, priorities=c('NM','XM'))	# Equals pickRefSeq.mRNA
mRNA
proteins <- pickRefSeq(refseq, priorities=c('NP','XP')) # Equals pickRefSeq.Protein
proteins


###################################################
### chunk number 10:  eval=FALSE
###################################################
## #line 178 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
## symbols <- c("SERPINA1","KERA","CD5")
## GOs <- translate(symbols, from=org.Bt.egSYMBOL2EG, to=org.Bt.egGO)
## # Pick biological process:
## pickGO(GOs, category='BP')
## # Pick only those biological processes for Entrez Gene ID 280730
## # which have been inferred from sequence similarity or electronic annotation:
## pickGO(translate(280730, org.Bt.egGO), category='BP', evidence=c('ISS','IEA'))


###################################################
### chunk number 11: 
###################################################
#line 198 "vignettes/AnnotationFuncs/inst/doc/Userguide.Rnw"
toLatex(sessionInfo())