I am using Gviz to plot several tracks. I create a GenomeAxisTrack()
and then several AnnotationTrack()
s that I feed to plotTracks()
. This produces a stacked output of tracks with annotated arrows as I want but the displays are linear. In some cases I have the complete genomes and I know they are circular (like some bacterial and viral chromosomes), how can I display the plots as circular?
I tried feeding to AnnotationTrack()
; a GRanges()
object with a Seqinfo()
object specifying the isCircular=TRUE
flag, but this didn't change the resulting plot.
Here is minimal code to create a linear plot, I would like one where the first and last arrow join up in a circle:
library(Gviz)
library(stringr)
track <- AnnotationTrack(start=c(1,5,7), end=c(2,6,10), strand=c('+','+','-'), stacking="dense", fill=c('red','blue', 'black'))
png(file=paste("test","png", sep="."))
plotTracks(track)
dev.off()
Note, I am not interested in alternative software recommendations, I would like to solve this with gviz if it is possible. Also asked this here.