psitecal: P-site position of reads

In kerenli/RiboDiPA_updating: Differential pattern analysis for Ribo-seq data

P-site position of reads

Description

This function outputs the P-site position, provided the CIGAR string of the alignment and the start position of read. It is implemented as a C++ function using the Rcpp package.

Usage

psiteCal(cigar, start, psitemap)

Arguments

cigar	A vector of CIGAR strings.
start	A vector of read start positions.
psitemap	A vector of relative P-site positions, which describe the offset from the 5 prime most nucleotide of the read to the P-site.

Value

A vector of P-site positions for the reads.

Examples

ex.cigar <- c("21M74731N7M", "2S11M57302N12M3S", "28M", "27M1S")
ex.start <- c(177640, 249163, 249286, 249290)
ex.psitemap <- c(9, 18, 9, 9)
psiteCal(ex.cigar, ex.start, ex.psitemap)

kerenli/RiboDiPA_updating documentation built on June 25, 2022, 9:50 p.m.