get_reads {polyester} | R Documentation

## 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

[Package polyester version 1.23.1 Index]