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
10X specific output fields:
CELL, C_CALL, CONSCOUNT, UMICOUNT, V_CALL_10X, D_CALL_10X, J_CALL_10X, JUNCTION_10X, JUNCTION_10X_AA

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 ...]
                      [--10x CELLRANGER_FILE [CELLRANGER_FILE ...]]
                      [--asis-id] [--asis-calls] [--partial] [--extended]
--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.

--10x <cellranger_file>

Table file containing 10X annotations (with .csv or .tsv extension).

--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.

--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.

--extended

Specify to include additional aligner specific fields in the output. Adds the <VDJ>_SCORE, <VDJ>_IDENTITY, <VDJ>_EVALUE, and <VDJ>_CIGAR; FWR1_IMGT, FWR2_IMGT, FWR3_IMGT, and FWR4_IMGT; CDR1_IMGT, CDR2_IMGT, and CDR3_IMGT; CDR3_IGBLAST_NT and CDR3_IGBLAST_AA (requires IgBLAST version 1.5 or greater).

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 ...]
                   [--10x CELLRANGER_FILE [CELLRANGER_FILE ...]] [--asis-id]
                   [--partial] [--extended]
--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.

--10x <cellranger_file>

Table file containing 10X annotations (with .csv or .tsv extension).

--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.

--extended

Specify to include additional aligner specific fields in the output. Adds the path score of the iHMMune-Align hidden Markov model as VDJ_SCORE; FWR1_IMGT, FWR2_IMGT, FWR3_IMGT, and FWR4_IMGT; CDR1_IMGT, CDR2_IMGT, and CDR3_IMGT.

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 ...]]
                   [--10x CELLRANGER_FILE [CELLRANGER_FILE ...]] [--asis-id]
                   [--partial] [--extended]
--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.

--10x <cellranger_file>

Table file containing 10X annotations (with .csv or .tsv extension).

--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.

--extended

Specify to include additional aligner specific fields in the output. Adds <VDJ>_SCORE and <VDJ>_IDENTITY>; FWR1_IMGT, FWR2_IMGT, FWR3_IMGT, and FWR4_IMGT; CDR1_IMGT, CDR2_IMGT, and CDR3_IMGT; and N1_LENGTH, N2_LENGTH, P3V_LENGTH, P5D_LENGTH, P3D_LENGTH, P5J_LENGTH, and D_FRAME.