NAME

App::SimulateReads::Command::Simulate::Genome - simulate subcommand class. Simulate genome sequencing

VERSION

version 0.11

SYNOPSIS

simulate_reads simulate genome [options] <fasta-file>

Arguments:
 a fasta-file 

Options:
 -h, --help               brief help message
 -M, --man                full documentation
 -v, --verbose            print log messages
 -p, --prefix             prefix output [default:"out"]	
 -o, --output-dir         output directory [default:"."]
 -j, --jobs               number of jobs [default:"1"; Integer]
 -z, --gzip               compress output file
 -c, --coverage           fastq-file coverage [default:"8", Number]
 -t, --sequencing-type    single-end or paired-end reads
                          [default:"paired-end"]
 -q, --quality-profile    illumina sequencing system profiles
                          [default:"hiseq"]
 -e, --sequencing-error   sequencing error rate
                          [default:"0.005"; Number]
 -r, --read-size          the read size [default:"101"; Integer]
 -m, --fragment-mean      the fragment mean size for paired-end reads
                          [default:"300"; Integer]
 -d, --fragment-stdd      the fragment standard deviation size for
                          paired-end reads [default:"50"; Integer]

DESCRIPTION

Simulate genome sequencing.

OPTIONS

--help

Print a brief help message and exits.

--man

Prints the manual page and exits.

--verbose

Prints log information to standard error

--prefix

Concatenates the prefix to the output-file name.

--output-dir

Creates output-file inside output-dir. If output-dir does not exist, it is created recursively

--jobs

Sets the number of child jobs to be created

--gzip

Compress the output-file with gzip algorithm. It is possible to pass --no-gzip if one wants uncompressed output-file

--read-size

Sets the read size. For now the unique valid value is 101

--coverage

Calculates the number of reads based on the sequence coverage: number_of_reads = (sequence_size * coverage) / read_size. This is the default option for genome sequencing simulation

--sequencing-type

Sets the sequencing type to single-end or paired-end

--fragment-mean

If the sequencing-type is set to paired-end, it sets the fragment mean

--fragment-stdd

If the sequencing-type is set to paired-end, it sets the fragment standard deviation

--sequencing-error

Sets the sequencing error rate. Valid values are between zero and one

--quality-profile

Sets the illumina sequencing system profile for quality. For now, the unique valid values are hiseq and poisson

AUTHOR

Thiago L. A. Miller <tmiller@mochsl.org.br>

COPYRIGHT AND LICENSE

This software is Copyright (c) 2018 by Teaching and Research Institute from Sírio-Libanês Hospital.

This is free software, licensed under:

The GNU General Public License, Version 3, June 2007