get_reads {polyester} | R Documentation |
simulate the sequencing process by returning the sequence of one or both ends of provided fragments
get_reads(tFrags, readlen, paired = TRUE)
tFrags |
DNAStringSet representing fragments |
readlen |
Read length. |
paired |
If |
DNAStringSet representing simulated RNA-seq reads
simulate_experiment
,
simulate_experiment_countmat
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