Quickly test if a BAM file has paired end reads
Iterate through a BAM file until a paired-end read is encountered or the end of file is reached; report the occurrence of paired-end reads to the user.
testPairedEndBam(file, index=file, ...)
file |
character(1) BAM file name, or a |
index |
(optional) character(1) name of the index file of the 'BAM' file being processed; this is given without the '.bai' extension. |
... |
Additional arguments, currently unused. |
A logical vector of length 1 containing TRUE is returned if BAM file contained paired end reads, FALSE otherwise.
Martin Morgan mailto:mtmorgan@fhcrc.org, Sonali Arora mailto:sarora@fhcrc.org
fl <- system.file("extdata", "ex1.bam", package="Rsamtools") testPairedEndBam(fl)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.