get_reads: get sequencing reads from fragments

Description Usage Arguments Value See Also Examples

Description

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

Usage

1
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

1
2
3
4
5
6
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

alyssafrazee/polyester-release documentation built on May 12, 2019, 2:32 a.m.