polyester (version 1.8.3)

get_reads: get sequencing reads from fragments

Description

simulate the sequencing process by returning the sequence of one or both ends of provided fragments

Usage

get_reads(tFrags, readlen, paired = TRUE)

Arguments

tFrags

DNAStringSet representing fragments

readlen

Read length.

paired

If FALSE, return only the first readlen bases of each element of tFrags in the result; if TRUE, also return last readlen bases.

Value

DNAStringSet representing simulated RNA-seq reads

See Also

simulate_experiment, simulate_experiment_countmat

Examples

library(Biostrings)
  data(srPhiX174)
  set.seed(174)
  srPhiX174_reads = get_reads(srPhiX174, readlen=15, paired=FALSE)
  srPhiX174_reads
  # set of single-end, 15bp reads, treating srPhiX174 as the fragments