Entering edit mode
stianlagstad
▴
90
@stianlagstad-9723
Last seen 4.8 years ago
The code below:
transcriptID <- "ENST00000395443" mart <- biomaRt::useMart( biomart = "ENSEMBL_MART_ENSEMBL", dataset = "hsapiens_gene_ensembl", host = "dec2013.archive.ensembl.org") # Create transcript track biomartTrack <- Gviz::BiomartGeneRegionTrack( filters = list( ensembl_transcript_id = transcriptID), biomart = mart) # Create axis track axisTrack <- Gviz::GenomeAxisTrack() displayPars(axisTrack) <- list( add53 = TRUE, labelPos = "below", cex = 1.5 ) # works as expected Gviz::plotTracks( c(biomartTrack, axisTrack), chromosome = "chr3")
Produces this plot:
And if I add reverseStrand=TRUE, I get this:
# add53 in same position? Gviz::plotTracks( c(biomartTrack, axisTrack), chromosome = "chr3", reverseStrand = TRUE)
As you can see, everything is reversed except for the 5' and 3' indicators. As the arrows on the axis change, shouldn't the add53 parameter react to reverseStrand also? Or should I react to that by using add35 instead?
The code in my original post above still produce the exact same plots (i.e. the 5' and 3' indicators are wrong with reverseStrand=TRUE). Using Gviz version 1.16. Could you check that the fix worked?