narrowAlignments: Narrow a set of aligned reads to a target region

Description Usage Arguments Value Author(s) Examples

Description

Aligned reads are narrowed to the target region. In the case of reads with deletions spanning the boundaries of the target, reads are narrowed to the start of the deletion,

Usage

1
2
3
4
5
narrowAlignments(alns, target, ...)

## S4 method for signature 'GAlignments,GRanges'
narrowAlignments(alns, target, ...,
  reverse.complement, verbose = FALSE)

Arguments

alns

A GAlignments object including a metadata column "seq" containing the sequence

target

A GRanges object

...

additional arguments

reverse.complement

Should the aligned reads be reverse complemented?

verbose

(Default: FALSE)

Value

The narrowed alignments (GAlignments),

Author(s)

Helen Lindsay

Examples

1
2
3
4
5
6
bam_fname <- system.file("extdata", "gol_F1_clutch_2_embryo_4_s.bam",
                         package = "CrispRVariants")
bam <- GenomicAlignments::readGAlignments(bam_fname, use.names = TRUE)
target <- GenomicRanges::GRanges("18", IRanges::IRanges(4647377, 4647399),
          strand = "+")
narrowAlignments(bam, target, reverse.complement = FALSE)

markrobinsonuzh/CrispRVariants documentation built on May 21, 2019, 12:23 p.m.