shiftReadsByFrame: Shift reads by reading frame

Description Usage Arguments Value Examples

View source: R/shiftReadsByFrame.R

Description

Shift reads P site position by reading frame. After shifting, all reading frame will be set as 0

Usage

1
shiftReadsByFrame(reads, txdb)

Arguments

reads

Output of getPsiteCoordinates

txdb

A TxDb object.

Value

Reads with reading frame information

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
library(Rsamtools)
bamfilename <- system.file("extdata", "RPF.WT.1.bam",
                           package="ribosomeProfilingQC")
yieldSize <- 10000000
bamfile <- BamFile(bamfilename, yieldSize = yieldSize)
pc <- getPsiteCoordinates(bamfile, bestpsite=11)
pc.sub <- pc[pc$qwidth %in% c(29, 30)]
library(GenomicFeatures)
library(BSgenome.Drerio.UCSC.danRer10)
txdb <- makeTxDbFromGFF(system.file("extdata",
          "Danio_rerio.GRCz10.91.chr1.gtf.gz",
          package="ribosomeProfilingQC"),
          organism = "Danio rerio",
          chrominfo = seqinfo(Drerio)["chr1"],
          taxonomyId = 7955)
pc.sub <- shiftReadsByFrame(pc.sub, txdb)

ribosomeProfilingQC documentation built on March 13, 2021, 2:01 a.m.