if(interactive()){
mart = useMart("ensembl", dataset="hsapiens_gene_ensembl")
genomeAxis = new("GenomeAxis", add53=TRUE)
plusStrand = new("GeneRegion", chromosome = "17", start = 30450000, end = 30550000, strand = "+", biomart = mart)
gdPlot(list(genomeAxis, plusStrand), minBase = 30450000, maxBase = 30550000)
}
Run the code above in your browser using DataLab