Skip to content

Commit 0925a6f

Browse files
committed
Adding new docker image with bwa-mem2 wrapper
1 parent e4aadb3 commit 0925a6f

File tree

2 files changed

+19
-4
lines changed

2 files changed

+19
-4
lines changed

config/containers.json

+2-2
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,7 @@
11
{
22
"images": {
33
"genome-seek": "docker://skchronicles/genome-seek:v0.1.0",
4-
"genome-seek_trim_map": "docker://skchronicles/genome-seek_trim_map:v0.1.0",
4+
"genome-seek_trim_map": "docker://skchronicles/genome-seek_trim_map:v0.1.1",
55
"genome-seek_somatic": "docker://skchronicles/genome-seek_somatic:v0.1.0",
66
"genome-seek_qc": "docker://skchronicles/genome-seek_qc:v0.1.0",
77
"genome-seek_sv": "docker://skchronicles/genome-seek_sv:v0.2.0",
@@ -19,4 +19,4 @@
1919
"sequenza": "docker://sequenza/sequenza:3.0.0",
2020
"vcf2maf": "docker://skchronicles/ncbr_vcf2maf:v0.1.0"
2121
}
22-
}
22+
}

workflow/rules/trim_map.smk

+17-2
Original file line numberDiff line numberDiff line change
@@ -76,8 +76,23 @@ rule bwa_mem2:
7676
if [ ! -d "{params.tmpdir}" ]; then mkdir -p "{params.tmpdir}"; fi
7777
tmp=$(mktemp -d -p "{params.tmpdir}")
7878
trap 'rm -rf "${{tmp}}"' EXIT
79-
80-
bwa-mem2 mem \\
79+
80+
# This wrapper to bwa-mem2 selects a
81+
# compatible pre-compiled binary for
82+
# Intel and AMD (Zen4) processors. The
83+
# AVX-512 instruction set was added to
84+
# processors with the release of Zen4.
85+
# This instruction set is NOT compatible
86+
# with the Intel-compiled AVX-512bw
87+
# binary in the docker image. The wrapper
88+
# detects whether it is running on an
89+
# Intel or AMD processor and selects the
90+
# appropriate binary. If AMD, it will use
91+
# the AVX2 binary. If Intel, it uses the
92+
# default entry point for bwa-mem2 which
93+
# will select the best available based
94+
# on available CPU information.
95+
bwa-mem2-wrapper mem \\
8196
-t {threads} \\
8297
-K 100000000 \\
8398
-M \\

0 commit comments

Comments
 (0)