Entering edit mode
                    stianlagstad
        
    
        ▴
    
    90
        @stianlagstad-9723
        Last seen 5.7 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?