Scroll to navigation

BWA-MEM2_MEM(1) User Commands BWA-MEM2_MEM(1)

NAME

bwa-mem2_mem - Sequence alignment using Burrows-Wheeler Transform

SYNOPSIS

bwa-mem2 mem [options] <idxbase> <in1.fq> [in2.fq]

DESCRIPTION

Looking to launch executable "/build/reproducible-path/bwa-mem2-2.2.1/./bwa-mem2.avx2", simd = .avx2 Launching executable "/build/reproducible-path/bwa-mem2-2.2.1/./bwa-mem2.avx2" ----------------------------- Executing in AVX2 mode!! ----------------------------- * SA compression enabled with xfactor: 8

OPTIONS

Algorithm options:
Output SAM file name
number of threads [1]
minimum seed length [19]
band width for banded alignment [100]
off-diagonal X-dropoff [100]
look for internal seeds inside a seed longer than {-k} * FLOAT [1.5]
seed occurrence for the 3rd round seeding [20]
skip seeds with more than INT occurrences [500]
drop chains shorter than FLOAT fraction of the longest overlapping chain [0.50]
discard a chain if seeded bases shorter than INT [0]
perform at most INT rounds of mate rescues for each read [50]
skip mate rescue
output file name missing
skip pairing; mate rescue performed unless -S also in use

Scoring options:

score for a sequence match, which scales options -TdBOELU unless overridden [1]
penalty for a mismatch [4]
gap open penalties for deletions and insertions [6,6]
gap extension penalty; a gap of size k cost '{-O} + {-E}*k' [1,1]
penalty for 5'- and 3'-end clipping [5,5]
penalty for an unpaired read pair [17]

Input/output options:

smart pairing (ignoring in2.fq)
read group header line such as '@RG\tID:foo\tSM:bar' [null]
insert STR to header if it starts with @; or insert lines in FILE [null]
treat ALT contigs as part of the primary assembly (i.e. ignore <idxbase>.alt file)
-5
for split alignment, take the alignment with the smallest coordinate as primary
don't modify mapQ of supplementary alignments
process INT input bases in each batch regardless of nThreads (for reproducibility) []
verbose level: 1=error, 2=warning, 3=message, 4+=debugging [3]
minimum score to output [30]
if there are <INT hits with score >80% of the max score, output all in XA [5,200]
output all alignments for SE or unpaired PE
append FASTA/FASTQ comment to SAM output
output the reference FASTA header in the XR tag
use soft clipping for supplementary alignments
mark shorter split hits as secondary

-I FLOAT[,FLOAT[,INT[,INT]]]

specify the mean, standard deviation (10% of the mean if absent), max (4 sigma from the mean if absent) and min of the insert size distribution. FR orientation only. [inferred]

Note: Please read the man page for detailed description of the command line and options.

April 2025 bwa-mem2 mem 2.2.1