getReadLength: determine the most common read length for input bam_file....

Description Usage Arguments Value Examples

View source: R/functions_fetch_bam.R

Description

determine the most common read length for input bam_file. uses 50 randomly selected regions from query_gr. If fewer than 20 reads are present, loads all of query_gr.

Usage

1
getReadLength(bam_file, query_gr)

Arguments

bam_file

indexed bam file

query_gr

GRanges to read from bam file

Value

numeric of most common read length.

Examples

1
2
3
qgr = CTCF_in_10a_overlaps_gr[1:5]
bam_file = system.file("extdata/test.bam", package = "seqsetvis", mustWork = TRUE)
getReadLength(bam_file, qgr)

seqsetvis documentation built on Nov. 8, 2020, 5:57 p.m.