dadaist2-rundada
Note
This is a new wrapper, introduced in 1.2.0, and experimental
Wrapper for DADA2 without any modification of the input reads. The input can be supplied either:
- As a single directory containing the reads (
-i DIRECTORY
), or - As two separate directory, one for the forward reads (
-f FOR_DIR
) and one for the reverse reads (-r REV_DIR
).
The latter is used as a compatibility layer and will be used by dadaist2 itself to invoke the wrapper.
Synopsis
usage: dadaist2-rundada [-h] [-i INPUT_DIR] [-f FOR_DIR] [-r REV_DIR] -o OUTPUT_DIR [--tmp TMP] [--fortag FORTAG] [--revertag REVERTAG] [--sample-separator SAMPLE_SEPARATOR]
[--sample-extension SAMPLE_EXTENSION] [-q TRUNC_QUAL] [-j] [-p] [--trunc-len-1 TRUNC_LEN_1] [--trunc-len-2 TRUNC_LEN_2] [--trim-left-1 TRIM_LEFT_1] [--trim-left-2 TRIM_LEFT_2]
[--maxee-1 MAXEE_1] [--maxee-2 MAXEE_2] [--chimera {none,pooled,consensus}] [--min-parent-fold MIN_PARENT_FOLD] [--n-learn N_LEARN] [-t THREADS] [--keep-temp] [--save-rds]
[--save-plots] [--log LOG] [--copy] [--skip-checks] [--verbose]
Run DADA2
optional arguments:
-h, --help show this help message and exit
Main:
-i INPUT_DIR, --input-dir INPUT_DIR
Input directory with both R1 and R2
-f FOR_DIR, --for-dir FOR_DIR
Input directory with R1 reads
-r REV_DIR, --rev-dir REV_DIR
Input directory with R2 reads
-o OUTPUT_DIR, --output-dir OUTPUT_DIR
Output directory
--tmp TMP Temporary directory
Input filtering:
--fortag FORTAG String defining a file as forward [default: _R1]
--revertag REVERTAG String defining a file as reverse [default: _R2]
--sample-separator SAMPLE_SEPARATOR
String acting as samplename separator [default: _]
--sample-extension SAMPLE_EXTENSION
String acting as samplename extension [default: .fastq.gz]
DADA2 parameters:
-q TRUNC_QUAL, --trunc-qual TRUNC_QUAL
Truncate at the first occurrence of a base with Q lower [default: 8]
-j, --join Join without merging
-p, --pool Pool samples
--trunc-len-1 TRUNC_LEN_1
Position at which to truncate forward reads [default: 0]
--trunc-len-2 TRUNC_LEN_2
Position at which to truncate reverse reads [default: 0]
--trim-left-1 TRIM_LEFT_1
Number of nucleotide to trim from the beginning of forward reads [default: 0]
--trim-left-2 TRIM_LEFT_2
Number of nucleotide to trim from the beginning of reverse reads [default: 0]
--maxee-1 MAXEE_1 Maximum expected errors in forward reads [default: 1.0]
--maxee-2 MAXEE_2 Maximum expected errors in reverse reads [default: 1.00]
--chimera {none,pooled,consensus}
Chimera handling can be none, pooled or consensus [default: pooled]
--min-parent-fold MIN_PARENT_FOLD
Minimum abundance of parents of a chimeric sequence (>1.0) [default: 1.0]
--n-learn N_LEARN Number of reads to learn the model, 0 for all [default: 0]
Other parameters:
-t THREADS, --threads THREADS
Number of threads
--keep-temp Keep temporary files
--save-rds Save RDS file with DADA2 output
--save-plots Save Quality plots of the input reads (PDF)
--log LOG Log file
--copy Copy input files instead of symbolic linking
--skip-checks Do not check installation of dependencies
--verbose Verbose mode
Output files
The output directory will contain:
- dada2.stats (table with the statistics of reads loss)
- dada2.tsv (main feature table)
- dada2.rds (R object with the table. if
--save-rds
is specified) - quality_R1.pdf, quality_R2.pdf (quality plots, if
--save-plots
is specified) - dada2.execution.log, dada2.execution.txt (wrapper log files)