shmulateSeq - Simulate mutations in a single sequence
Generates random mutations in a sequence iteratively using a targeting model. Targeting probabilities at each position are updated after each iteration.
shmulateSeq(sequence, mutations, targetingModel = HH_S5F)
- sequence string in which mutations are to be introduced.
- number of mutations to be introduced into
- 5-mer TargetingModel object to be used for computing probabilities of mutations at each position. Defaults to HH_S5F.
A string defining the mutated sequence.
# Define example input sequence sequence <- "NGATCTGACGACACGGCCGTGTATTACTGTGCGAGAGATAGTTTA" # Simulate using the default human 5-mer targeting model shmulateSeq(sequence, mutations=6)