Advertisement
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
Advertisement