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