getATACfragments

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

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

Fabian Mueller