MakeDb

Create tab-delimited database file to store sequence alignment information

usage: MakeDb [--version] [-h]  ...
--version

show program’s version number and exit

-h, --help

show this help message and exit

output files:
db-pass
database of alignment records with functionality information, V and J calls, and a junction region.
db-fail
database with records that fail due to no functionality information (did not pass IMGT), no V call, no J call, or no junction region.
universal output fields:
SEQUENCE_ID, SEQUENCE_INPUT, SEQUENCE_VDJ, SEQUENCE_IMGT, FUNCTIONAL, IN_FRAME, STOP, MUTATED_INVARIANT, INDELS, V_CALL, D_CALL, J_CALL, V_SEQ_START, V_SEQ_LENGTH, D_SEQ_START, D_SEQ_LENGTH, D_GERM_START, D_GERM_LENGTH, J_SEQ_START, J_SEQ_LENGTH, J_GERM_START, J_GERM_LENGTH, NP1_LENGTH, NP2_LENGTH, JUNCTION_LENGTH, JUNCTION, GERMLINE_IMGT, FWR1_IMGT, FWR2_IMGT, FWR3_IMGT, FWR4_IMGT, CDR1_IMGT, CDR2_IMGT, CDR3_IMGT
imgt specific output fields:
V_GERM_START_IMGT, V_GERM_LENGTH_IMGT, N1_LENGTH, N2_LENGTH, P3V_LENGTH, P5D_LENGTH, P3D_LENGTH, P5J_LENGTH, D_FRAME, V_SCORE, V_IDENTITY, J_SCORE, J_IDENTITY,
igblast specific output fields:
V_GERM_START_VDJ, V_GERM_LENGTH_VDJ, V_EVALUE, V_SCORE, V_IDENTITY, V_CIGAR, D_EVALUE, D_SCORE, D_IDENTITY, D_CIGAR, J_EVALUE, J_SCORE, J_IDENTITY, J_CIGAR, CDR3_IGBLAST, CDR3_IGBLAST_AA
ihmm specific output fields:
V_GERM_START_VDJ, V_GERM_LENGTH_VDJ, VDJ_SCORE

MakeDb igblast

Process IgBLAST output.

usage: MakeDb igblast [--version] [-h] [-o OUT_FILES [OUT_FILES ...]]
                      [--outdir OUT_DIR] [--outname OUT_NAME] [--log LOG_FILE]
                      [--failed] [--format {changeo,airr}] -i ALIGNER_FILES
                      [ALIGNER_FILES ...] -r REPO [REPO ...] -s SEQ_FILES
                      [SEQ_FILES ...] [--partial] [--scores] [--regions]
                      [--cdr3] [--asis-id] [--asis-calls]
--version

show program’s version number and exit

-h, --help

show this help message and exit

-o <out_files>

Explicit output file name. 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.

--format {changeo,airr}

Specify input and output format.

-i <aligner_files>

IgBLAST output files in format 7 with query sequence (IgBLAST argument ‘-outfmt “7 std qseq sseq btop”’).

-r <repo>

List of folders and/or fasta files containing the same germline set used in the IgBLAST alignment. These reference sequences must contain IMGT-numbering spacers (gaps) in the V segment.

-s <seq_files>

List of input FASTA files (with .fasta, .fna or .fa extension), containing sequences.

--partial

If specified, include incomplete V(D)J alignments in the pass file instead of the fail file. An incomplete alignment is defined as a record for which a valid IMGT-gapped sequence cannot be built or that is missing a V gene assignment, J gene assignment, junction region, or productivity call.

--scores

Specify if alignment score metrics should be included in the output. Adds the <VDJ>_SCORE, <VDJ>_IDENTITY, <VDJ>_EVALUE, <VDJ>_CIGAR columns.

--regions

Specify if IMGT FWR and CDRs should be included in the output. Adds the FWR1_IMGT, FWR2_IMGT, FWR3_IMGT, FWR4_IMGT, CDR1_IMGT, CDR2_IMGT, and CDR3_IMGT columns.

--cdr3

Specify if the CDR3 sequences generated by IgBLAST should be included in the output. Adds the columns CDR3_IGBLAST_NT and CDR3_IGBLAST_AA. Requires IgBLAST version 1.5 or greater.

--asis-id

Specify to prevent input sequence headers from being parsed to add new columns to database. Parsing of sequence headers requires headers to be in the pRESTO annotation format, so this should be specified when sequence headers are incompatible with the pRESTO annotation scheme. Note, unrecognized header formats will default to this behavior.

--asis-calls

Specify to prevent gene calls from being parsed into standard allele names in both the IgBLAST output and reference database. Note, this requires the sequence identifiers in the reference sequence set and the IgBLAST database to be exact string matches.

MakeDb ihmm

Process iHMMune-Align output.

usage: MakeDb ihmm [--version] [-h] [-o OUT_FILES [OUT_FILES ...]]
                   [--outdir OUT_DIR] [--outname OUT_NAME] [--log LOG_FILE]
                   [--failed] [--format {changeo,airr}] -i ALIGNER_FILES
                   [ALIGNER_FILES ...] -r REPO [REPO ...] -s SEQ_FILES
                   [SEQ_FILES ...] [--asis-id] [--partial] [--scores]
                   [--regions]
--version

show program’s version number and exit

-h, --help

show this help message and exit

-o <out_files>

Explicit output file name. 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.

--format {changeo,airr}

Specify input and output format.

-i <aligner_files>

iHMMune-Align output file.

-r <repo>

List of folders and/or FASTA files containing the set of germline sequences used by iHMMune-Align. These reference sequences must contain IMGT-numbering spacers (gaps) in the V segment.

-s <seq_files>

List of input FASTA files (with .fasta, .fna or .fa extension) containing sequences.

--asis-id

Specify to prevent input sequence headers from being parsed to add new columns to database. Parsing of sequence headers requires headers to be in the pRESTO annotation format, so this should be specified when sequence headers are incompatible with the pRESTO annotation scheme. Note, unrecognized header formats will default to this behavior.

--partial

If specified, include incomplete V(D)J alignments in the pass file instead of the fail file. An incomplete alignment is defined as a record for which a valid IMGT-gapped sequence cannot be built or that is missing a V gene assignment, J gene assignment, junction region, or productivity call.

--scores

Specify if alignment score metrics should be included in the output. Adds the path score of the iHMMune-Align hidden Markov model to VDJ_SCORE.

--regions

Specify if IMGT FWRs and CDRs should be included in the output. Adds the FWR1_IMGT, FWR2_IMGT, FWR3_IMGT, FWR4_IMGT, CDR1_IMGT, CDR2_IMGT, and CDR3_IMGT columns.

MakeDb imgt

Process IMGT/HighV-Quest output
(does not work with V-QUEST).
usage: MakeDb imgt [--version] [-h] [-o OUT_FILES [OUT_FILES ...]]
                   [--outdir OUT_DIR] [--outname OUT_NAME] [--log LOG_FILE]
                   [--failed] [--format {changeo,airr}]
                   [-i ALIGNER_FILES [ALIGNER_FILES ...]]
                   [-s [SEQ_FILES [SEQ_FILES ...]]] [-r REPO [REPO ...]]
                   [--asis-id] [--partial] [--scores] [--regions] [--junction]
--version

show program’s version number and exit

-h, --help

show this help message and exit

-o <out_files>

Explicit output file name. 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.

--format {changeo,airr}

Specify input and output format.

-i <aligner_files>

Either zipped IMGT output files (.zip or .txz) or a folder containing unzipped IMGT output files (which must include 1_Summary, 2_IMGT-gapped, 3_Nt-sequences, and 6_Junction).

-s <seq_files>

List of FASTA files (with .fasta, .fna or .fa extension) that were submitted to IMGT/HighV-QUEST. If unspecified, sequence identifiers truncated by IMGT/HighV-QUEST will not be corrected.

-r <repo>

List of folders and/or fasta files containing the germline sequence set used by IMGT/HighV-QUEST. These reference sequences must contain IMGT-numbering spacers (gaps) in the V segment. If unspecified, the germline sequence reconstruction will not be included in the output.

--asis-id

Specify to prevent input sequence headers from being parsed to add new columns to database. Parsing of sequence headers requires headers to be in the pRESTO annotation format, so this should be specified when sequence headers are incompatible with the pRESTO annotation scheme. Note, unrecognized header formats will default to this behavior.

--partial

If specified, include incomplete V(D)J alignments in the pass file instead of the fail file. An incomplete alignment is defined as a record that is missing a V gene assignment, J gene assignment, junction region, or productivity call.

--scores

Specify if alignment score metrics should be included in the output. Adds the <VDJ>_SCORE and <VDJ>_IDENTITY> columns.

--regions

Specify if IMGT FWRs and CDRs should be included in the output. Adds the FWR1_IMGT, FWR2_IMGT, FWR3_IMGT, FWR4_IMGT, CDR1_IMGT, CDR2_IMGT, and CDR3_IMGT columns.

--junction

Specify if detailed junction fields should be included in the output. Adds the columns N1_LENGTH, N2_LENGTH, P3V_LENGTH, P5D_LENGTH, P3D_LENGTH, P5J_LENGTH, D_FRAME.