BuildConsensus

Builds a consensus sequence for each set of input sequences

usage: BuildConsensus [--version] [-h] -s SEQ_FILES [SEQ_FILES ...]
                      [-o OUT_FILES [OUT_FILES ...]] [--outdir OUT_DIR]
                      [--outname OUT_NAME] [--log LOG_FILE] [--failed]
                      [--fasta] [--delim DELIMITER DELIMITER DELIMITER]
                      [--nproc NPROC] [-n MIN_COUNT] [--bf BARCODE_FIELD]
                      [-q MIN_QUAL] [--freq MIN_FREQ] [--maxgap MAX_GAP]
                      [--pf PRIMER_FIELD] [--prcons PRIMER_FREQ]
                      [--cf COPY_FIELDS [COPY_FIELDS ...]]
                      [--act {min,max,sum,set,majority} [{min,max,sum,set,majority} ...]]
                      [--dep] [--maxdiv MAX_DIVERSITY | --maxerror MAX_ERROR]
--version

show program’s version number and exit

-h, --help

show this help message and exit

-s <seq_files>

A list of FASTA/FASTQ files containing sequences to process.

-o <out_files>

Explicit output file name(s). Note, this argument cannot be used with the –failed, –outdir, or –outname arguments. If unspecified, then the output filename will be based on the input filename(s).

--outdir <out_dir>

Specify to changes the output directory to the location specified. The input file directory is used if this is not specified.

--outname <out_name>

Changes the prefix of the successfully processed output file to the string specified. May not be specified with multiple input files.

--log <log_file>

Specify to write verbose logging to a file. May not be specified with multiple input files.

--failed

If specified create files containing records that fail processing.

--fasta

Specify to force output as FASTA rather than FASTQ.

--delim <delimiter>

A list of the three delimiters that separate annotation blocks, field names and values, and values within a field, respectively.

--nproc <nproc>

The number of simultaneous computational processes to execute (CPU cores to utilized).

-n <min_count>

The minimum number of sequences needed to define a valid consensus.

--bf <barcode_field>

Position of description barcode field to group sequences by.

-q <min_qual>

Consensus quality score cut-off under which an ambiguous character is assigned; does not apply when quality scores are unavailable.

--freq <min_freq>

Fraction of character occurrences under which an ambiguous character is assigned.

--maxgap <max_gap>

If specified, this defines a cut-off for the frequency of allowed gap values for each position. Positions exceeding the threshold are deleted from the consensus. If not defined, positions are always retained.

--pf <primer_field>

Specifies the field name of the primer annotations

--prcons <primer_freq>

Specify to define a minimum primer frequency required to assign a consensus primer, and filter out sequences with minority primers from the consensus building step.

--cf <copy_fields>

Specifies a set of additional annotation fields to copy into the consensus sequence annotations.

--act {min,max,sum,set,majority}

List of actions to take for each copy field which defines how each annotation will be combined into a single value. The actions “min”, “max”, “sum” perform the corresponding mathematical operation on numeric annotations. The action “set” combines annotations into a comma delimited list of unique values and adds an annotation named <FIELD>_COUNT specifying the count of each item in the set. The action “majority” assigns the most frequent annotation to the consensus annotation and adds an annotation named <FIELD>_FREQ specifying the frequency of the majority value.

--dep

Specify to calculate consensus quality with a non-independence assumption

--maxdiv <max_diversity>

Specify to calculate the nucleotide diversity of each read group (average pairwise error rate) and remove groups exceeding the given diversity threshold. Diversity is calculate for all positions within the read group, ignoring any character filtering imposed by the -q, –freq and –maxgap arguments. Mutually exclusive with –maxerror.

--maxerror <max_error>

Specify to calculate the error rate of each read group (rate of mismatches from consensus) and remove groups exceeding the given error threshold. The error rate is calculated against the final consensus sequence, which may include masked positions due to the -q and –freq arguments and may have deleted positions due to the –maxgap argument. Mutually exclusive with –maxdiv.

output files:
consensus-pass
consensus reads.
consensus-fail
raw reads failing consensus filtering criteria.
output annotation fields:
PRIMER
a comma delimited list of unique primer annotations found within the barcode read group.
PRCOUNT
a comma delimited list of the corresponding counts of unique primer annotations.
PRCONS
the majority primer within the barcode read group.
PRFREQ
the frequency of the majority primer.
CONSCOUNT
the count of reads within the barcode read group which contributed to the consensus sequence. This is the total size of the read group, minus sequence excluded due to user defined filtering criteria.