getReadsInRange: Retrieves Reads in Specified Range

Description Usage Arguments Value Examples

Description

'getReadsInRange' retrieves reads in range based on arguments refSeq, refStart, refEnd and idx.

Usage

1
getReadsInRange(cmpH5, refSeq, refStart = 1, refEnd = NA, idx = 1:nrow(cmpH5))

Arguments

cmpH5

An object of class PacBioCmpH5.

refSeq

Reference sequence from which to retrieve reads.

refStart

Reference start point from which to retrieve reads.

refEnd

Reference end point from which to retrieve reads.

idx

The indices of the alignments to return.

Value

'getReadsInRange' returns a vector of integers.

Examples

1
2
3
4
  cmpH5 <- PacBioCmpH5(system.file("h5_files",  "aligned_reads.cmp.h5", package = "pbh5"))
  readsInRange <- getReadsInRange(cmpH5, 1, 100, 200)
  class(readsInRange)
  head(readsInRange)

PacificBiosciences/R-pbh5 documentation built on May 7, 2019, 11:54 p.m.