Usage
simulate_experiment_countmat(fasta = NULL, gtf = NULL, seqpath = NULL, readmat, outdir = ".", fraglen = 250, fragsd = 25, readlen = 100, error_rate = 0.005, paired = TRUE, seed = NULL, ...)
Arguments
fasta
path to FASTA file containing transcripts from which to simulate
reads. See details.
gtf
path to GTF file containing transcript structures from which reads
should be simulated. See details.
seqpath
path to folder containing one FASTA file (.fa
extension) for each chromosome in gtf
. See details.
readmat
matrix with rows representing transcripts and columns
representing samples. Entry i,j specifies how many reads to simulate from
transcript i for sample j.
outdir
character, path to folder where simulated reads should be
written, without a slash at the end of the folder name. By default, reads
written to the working directory.
fraglen
Mean RNA fragment length. Sequences will be read off the
end(s) of these fragments.
fragsd
Standard deviation of fragment lengths.
error_rate
Sequencing error rate. Must be between 0 and 1. A uniform
error model is assumed.
paired
If TRUE
, paired-end reads are simulated; else single-end
reads are simulated.
seed
Optional seed to set before simulating reads, for
reproducibility.
...
Further arguments to pass to seq_gtf
, if gtf
is not
NULL
.