I trying to make a nice plot such as in here:
I mean, chromosome at the top and gene (I need only one gene) below, exons presetned as blocks and introns as line with arrows.
My plot looks like this:
and my code:
library(org.Hs.eg.db) eid <- select(org.Hs.eg.db, "GNG7", "ENTREZID", "SYMBOL")[["ENTREZID"]] txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene txid <- select(txdb, eid, "TXNAME", "GENEID")[["TXNAME"]] exons <- exonsBy(txdb, by="tx", use.names=TRUE) gng7exons <- exons[names(exons) %in% txid] gng7exons[] anno <- AnnotationTrack(gng7exons) plotTracks(list(gtrack,itrack,anno), from = 2645916, to = 2646703, cex = 0.8)
It is even possible to add SNP's, cpg islands, etc to this plot? Summarizing, I'm looking for something like in UCSC and I have no idea how to get it...