rangeToFastq: Extract read information from alignments for given genomic range as fastq.

Description

The function rangeToFastq writes all (or selected) aligns from a bamRange into a compressed fastq file.

Usage

1

Arguments

object

bamReader. Must be opened.

filename

Name of output 'fastq' file.

which

(Optional) Logical vector. When given, the routine checks for each of the given values the value of which. When TRUE the corresponding align will be written to the output file, otherwise will be skipped.

append

(Optional) When TRUE, the routine will append to an existing file. Otherwise existing files will be overwritten.

Details

range2fastq ist deprecated and will be replaced by rangeToFastq soon.

Author(s)

Wolfgang Kaisers

Examples

1
2
3
4
5
6
7
bam <- system.file("extdata", "accepted_hits.bam", package="rbamtools")
reader <- bamReader(bam, idx=TRUE)
coords <- as.integer(c(0, 0, 249250621))
range <- bamRange(reader, coords)
## Not run: rangeToFastq(range, "range.fq.gz")
index <- sample(1:size(range), 100)
## Not run: rangeToFastq(range, "range_subset.fq.gz", which=index)

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.