Search
Question: Gviz: Possible 5' and 3' marker errors
0
gravatar for stianlagstad
11 months ago by
stianlagstad90
stianlagstad90 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 11 months ago • written 11 months ago by stianlagstad90
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: 330 users visited in the last hour