Question: Rice genome ideogram construction
gravatar for sabrina.elias
2.9 years ago by
sabrina.elias0 wrote:

Hi All

I was having difficulty in constructing an ideogram for rice genome using gviz, since in UCSC rice genome is not included. I also wanted to export annotation information for specific regions. As gviz provide a nice view for this presentation i was trying to use it, but could not do it properly. 

Florian helped me by modifying my codes as below.

Posting it here thinking it might be helpful for someone else also. 

I was having difficulty with gieStain option also. If anyone can help with that.







idTable <- data.frame(chrom=tolower(seqnames(BSgenome.Osativa.MSU.MSU7)), chromStart=0, chromEnd=seqlengths(BSgenome.Osativa.MSU.MSU7),name=seqnames(BSgenome.Osativa.MSU.MSU7), gieStain="gneg", stringsAsFactors=FALSE)

idTrack <- IdeogramTrack(bands=idTable, genome="MSU7")


mart = useEnsembl(biomart="plants_mart", host="", dataset = "osativa_eg_gene")

bm <- useMart(host = "", biomart="plants_mart", dataset = "osativa_eg_gene")

fm <- Gviz:::.getBMFeatureMap()

fm["cdsl"] <- c("cds_start", "cds_length")

fm["symbol"] <- c("external_gene_id", "external_gene_name")

biomTrack <- BiomartGeneRegionTrack(chromosome = 1, start = 20000000, end = 21000000, biomart=bm, featureMap=fm)

plotTracks(list(idTrack, biomTrack), transcriptAnnotation="symbol")


ideogram gviz plants_mart • 479 views
ADD COMMENTlink written 2.9 years ago by sabrina.elias0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 276 users visited in the last hour