Function shifts footprints (GRanges) using specified offsets for every of the specified lengths. Reads that do not conform to the specified lengths are filtered out and rejected. Reads are resized to single base in 5' end fashion, treated as p site. This function takes account for junctions in cigars of the reads. Length of the footprint is saved in size' parameter of GRanges output. Footprints are also sorted according to their genomic position, ready to be saved as a ofst, bed or wig file.
a data.frame / data.table with minimum 2 columns,
fraction (selected_lengths) and selected_shifts (relative position).
logical, default TRUE. If False will keep original order of reads, and not sort output reads in increasing genomic location per chromosome and strand.
The two columns in the shift data.frame/data.table argument are:
- fraction Numeric vector of lengths of footprints you select for shifting.
- offsets_start Numeric vector of shifts for corresponding selected_lengths. eg. c(-10, -10) with selected_lengths of c(31, 32) means length of 31 will be shifted left by 10. Footprints of length 32 will be shifted right by 10.
NOTE: It will remove softclips from valid width, the CIGAR 3S30M is qwidth 33, but will remove 3S so final read width is 30 in ORFik.
GRanges object of shifted footprints, sorted and
resized to 1bp of p-site,
with metacolumn "size" indicating footprint size before shifting and
resizing, sorted in increasing order.
1 2 3 4 5 6 7 8 9 10 11 12 13 14
## Basic run # Transcriptome annotation -> gtf_file <- system.file("extdata", "annotations.gtf", package = "ORFik") # Ribo seq data -> riboSeq_file <- system.file("extdata", "ribo-seq.bam", package = "ORFik") ## Not run: footprints <- readBam(riboSeq_file) # detect the shifts automagically shifts <- detectRibosomeShifts(footprints, gtf_file) # shift the RiboSeq footprints shiftedReads <- shiftFootprints(footprints, shifts) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.