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)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.