readAlignmentatCpGIsland: Reads a sequence alignment from a BAM file at the given...

Description Usage Arguments Value Author(s)

Description

Reads paired-end sequences from a BAM file given the coordinates of a CpG Isalnd and returns the corresponding sequence alignment and base qualities

Usage

1
readAlignmentatCpGIsland(bam, bai, chromosome, start, end, strand)

Arguments

bam

character string specifying path to BAM file

bai

character string specifying path to index file. Should not contain extension

chromosome

character [X, Y, MT] or numeric [1-22] string of query chomosome

start

numeric string defining start of genomic coordinate to query

end

numeric string defining end of genomic coordinate to query

strand

character string specifying which reads to return. Forward (+) or reverse (-)

Value

Invisibly a list containing a GAlignmentPairs metadata of the sequences returned and two DNAStringSet objects of sequences and base qualities

Author(s)


David Brown
brownd7@mskcc.org


ndbrown6/EMseq documentation built on Jan. 20, 2021, 2:36 a.m.