Struct Simulate 

pub struct Simulate {

    pub reference: ReferenceOptions,
    pub vcf: VcfOptions,
    pub bed: BedOptions,
    pub output: OutputPrefixOptions,
    pub seed: SeedOptions,
    pub coverage: f64,
    pub read_length: usize,
    pub fragment_mean: usize,
    pub fragment_stddev: usize,
    pub min_fragment_length: usize,
    pub single_end: bool,
    pub adapter_r1: String,
    pub adapter_r2: String,
    pub min_error_rate: f64,
    pub max_error_rate: f64,
    pub golden_bam: bool,
    pub golden_vcf: bool,
    pub simple_names: bool,
    pub compression: u8,
    pub threads: usize,
}

Simulate sequencing reads from a reference genome.

Generates paired-end or single-end FASTQ files with optional ground-truth BAM and VCF outputs for benchmarking alignment and variant calling pipelines. Variants are applied from an input VCF to construct haplotype sequences, and reads are sampled with a position-dependent Illumina error model.

Fields

reference: ReferenceOptions

vcf: VcfOptions

bed: BedOptions

output: OutputPrefixOptions

seed: SeedOptions

coverage: f64

Mean sequencing coverage depth.

read_length: usize

Length of each read in bases.

fragment_mean: usize

Mean insert size (outer distance between read pair ends).

fragment_stddev: usize

Standard deviation of the insert size distribution.

min_fragment_length: usize

Minimum fragment length in bases; sampled lengths below this are clamped up. Fragments shorter than the read length are padded with adapter sequence. Must be at least 1.

single_end: bool

Generate single-end reads instead of paired-end.

adapter_r1: String

Adapter sequence appended to read 1 when the fragment is shorter than the read length.

adapter_r2: String

Adapter sequence appended to read 2 when the fragment is shorter than the read length.

min_error_rate: f64

Minimum per-base error rate, applied at the start of reads.

max_error_rate: f64

Maximum per-base error rate, applied at the end of reads.

golden_bam: bool

Write a ground-truth BAM file with correct alignments.

golden_vcf: bool

Write a ground-truth VCF annotated with simulated coverage.

simple_names: bool

Use simple read names (holodeck:N) instead of encoding truth coordinates in the read name.

compression: u8

BGZF compression level (0-12). Lower values are faster with larger output files; higher values produce smaller files at the cost of speed. Level 0 is no compression; 1 is fastest; 12 is maximum compression.

threads: usize

Number of threads for parallel BGZF output compression.

Trait Implementations

impl Args for Simulate

fn group_id() -> Option<Id>

Report the ArgGroup::id for this set of arguments

fn augment_args<'b>(__clap_app: Command) -> Command

Append to Command so it can instantiate Self via FromArgMatches::from_arg_matches_mut

fn augment_args_for_update<'b>(__clap_app: Command) -> Command

Append to Command so it can instantiate self via FromArgMatches::update_from_arg_matches_mut

impl Command for Simulate

fn execute(&self) -> Result<()>

Execute the command.

impl CommandFactory for Simulate

fn command<'b>() -> Command

Build a Command that can instantiate Self.

fn command_for_update<'b>() -> Command

Build a Command that can update self.

impl Debug for Simulate

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl FromArgMatches for Simulate

fn from_arg_matches(__clap_arg_matches: &ArgMatches) -> Result<Self, Error>

Instantiate Self from ArgMatches, parsing the arguments as needed.

fn from_arg_matches_mut( __clap_arg_matches: &mut ArgMatches, ) -> Result<Self, Error>

Instantiate Self from ArgMatches, parsing the arguments as needed.

fn update_from_arg_matches( &mut self, __clap_arg_matches: &ArgMatches, ) -> Result<(), Error>

Assign values from ArgMatches to self.

fn update_from_arg_matches_mut( &mut self, __clap_arg_matches: &mut ArgMatches, ) -> Result<(), Error>

Assign values from ArgMatches to self.

impl Parser for Simulate

fn parse() -> Self

Parse from std::env::args_os(), exit on error.

fn try_parse() -> Result<Self, Error>

Parse from std::env::args_os(), return Err on error.

fn parse_from<I, T>(itr: I) -> Self

where I: IntoIterator<Item = T>, T: Into<OsString> + Clone,

Parse from iterator, exit on error.

fn try_parse_from<I, T>(itr: I) -> Result<Self, Error>

where I: IntoIterator<Item = T>, T: Into<OsString> + Clone,

Parse from iterator, return Err on error.

fn update_from<I, T>(&mut self, itr: I)

where I: IntoIterator<Item = T>, T: Into<OsString> + Clone,

Update from iterator, exit on error.

fn try_update_from<I, T>(&mut self, itr: I) -> Result<(), Error>

where I: IntoIterator<Item = T>, T: Into<OsString> + Clone,

Update from iterator, return Err on error.