Search
Question: Gviz: Possible 5' and 3' marker errors
0
gravatar for stianlagstad
4 weeks ago by
stianlagstad70
stianlagstad70 wrote:

This code:
```
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
)

Gviz::plotTracks(
  c(biomartTrack, axisTrack),
  chromosome = "chr3",
  add53 = TRUE,
  add35 = TRUE)
```

Produces this:
Both strands are tagged with 5' on the left and 3' on the right, which is wrong. Is this a bug?

Thank you.

 

ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by stianlagstad70
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 114 users visited in the last hour