random_seq: Create random sequences

In FastqCleaner: A Shiny Application for Quality Control, Filtering and Trimming of FASTQ Files

Description

Create a DNAStringSet object with random sequences

Usage

random_seq(slength, swidth, nuc = c("DNA", "RNA"), prob = c(0.25, 0.25,
  0.25, 0.25))

Arguments

slength	Number of sequences
swidth	Width of the sequences
nuc	Create sequences of DNA (nucleotides = c('A', 'C', 'G', 'T')) or RNA (nucleotides = c('A, 'C', 'G', 'U'))?. Default: 'DNA'
prob	A vector of four probability values used to set the frequency of the nucleotides 'A', 'C', 'G', 'T', for DNA, or 'A', 'C', 'G', 'U', for RNA. For example = c(0.25, 0.25, 0.5, 0). Default is = c(0.25, 0.25, 0.25, 0.25) (equiprobability for the 4 bases). If the sum of the probabilities is > 1, the values will be nomalized to the range [0, 1].

Value

DNAStringSet object

Author(s)

Leandro Roser [e-mail]

Examples

# For reproducible examples, make a call to set.seed before 
# running each random function

set.seed(10)
s1 <- random_seq(slength = 10, swidth = 20)
s1

set.seed(10)
s2 <- random_seq(slength = 10, swidth = 20, 
prob = c(0.6, 0.1, 0.3, 0))
s2

FastqCleaner documentation built on Nov. 8, 2020, 5:05 p.m.