Random sequences are generated to create a new ref_genome object. Note that this function will never generate empty sequences.

create_genome(n_seqs, len_mean, len_sd = 0, pi_tcag = rep(0.25, 4),
n_threads = 1)

## Arguments

n_seqs Number of sequences. Mean for the gamma distribution of sequence sizes. Standard deviation for the gamma distribution of sequence sizes. If set to <= 0, all sequences will be the same length. Defaults to 0. Vector of length 4 containing the nucleotide equilibrium frequencies for "T", "C", "A", and "G", respectively. Defaults to rep(0.25, 4). Number of threads to use for parallel processing. This argument is ignored if OpenMP is not enabled. Defaults to 1.

## Value

A ref_genome object.

## Examples


genome <- create_genome(10, 100e3, 100, pi_tcag = c(0.1, 0.2, 0.3, 0.4))