getATACfragments: getATACfragments

View source: R/utils_atac.R

getATACfragmentsR Documentation

getATACfragments

Description

Given a GAlignmentPairs or GAlignments object, return a GRanges object containing the fragment (or insertion site for single-end data)

Usage

getATACfragments(ga, offsetTn = TRUE)

Arguments

ga

GAlignmentPairs (or GAlignments for single-end sequencing) object

offsetTn

apply offsets for Tn5 dimer cut site (+4 bp on genomic + strand; -4 bp on genomic - strand)

Value

GRanges object containing derived insertions. For paired-end data (recommended), the width of the resulting ranges corresponds to the insert size for single-end data, the width is set to 1bp

Author(s)

Fabian Mueller


GreenleafLab/ChrAccR documentation built on March 22, 2023, 11:42 p.m.