alpine (version 0.99.1)

getReadLength: Get read length

Description

Gets the length of the first read in a BAM file

Usage

getReadLength(bam.files)

Arguments

bam.files

a character vector pointing to BAM files

Value

a numeric vector, one number per BAM file, the length of the first read in the file

Examples

# these next lines just write out a BAM file from R
# typically you would already have a BAM file
library(alpineData)
library(GenomicAlignments)
library(rtracklayer)
gap <- ERR188088()
dir <- system.file(package="alpineData", "extdata")
bam.file <- c("ERR188088" = file.path(dir,"ERR188088.bam"))
export(gap, con=bam.file)

getReadLength(bam.file)