Not a member of Pastebin yet?
                        Sign Up,
                        it unlocks many cool features!                    
                - # replicate highlight box reverseStrand error
 - 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()
 - # works as expected
 - Gviz::plotTracks(
 - c(biomartTrack, axisTrack),
 - chromosome = "chr3")
 - # works as expected
 - Gviz::plotTracks(
 - c(biomartTrack, axisTrack),
 - chromosome = "chr3",
 - reverseStrand = TRUE)
 - # Create highlight track
 - highlightTrack <- Gviz::HighlightTrack(
 - trackList = biomartTrack,
 - start = biomartTrack@start,
 - width = 2000,
 - chromosome = "chr3")
 - # works as expected
 - Gviz::plotTracks(c(highlightTrack, axisTrack))
 - # this only reverses the transcript, not the highlighttrack
 - Gviz::plotTracks(
 - c(highlightTrack, axisTrack),
 - reverseStrand = TRUE)
 
Advertisement
 
                    Add Comment                
                
                        Please, Sign In to add comment