get_reads: get sequencing reads from fragments

Description Usage Arguments Value See Also Examples

View source: R/get_reads.R

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

polyester documentation built on Nov. 8, 2020, 8:09 p.m.