We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent df76adb commit ce4f079Copy full SHA for ce4f079
1 file changed
src/bash_workflows/check_aligner.sh
@@ -0,0 +1,26 @@
1
+#!/bin/bash
2
+
3
+set -euo pipefail
4
5
+reads_path=$1
6
+reads_index=$2
7
+sample_id=$3
8
+reference_fasta=${4:-""}
9
+reference_index=${5:-""}
10
+reference_dict=${6:-""}
11
12
+if [[ -n "$reference_fasta" ]]; then
13
+ ref_opt="-R $reference_fasta"
14
+else
15
+ ref_opt=""
16
+fi
17
18
+gatk PrintReadsHeader \
19
+ -I "${reads_path}" \
20
+ --read-index "${reads_index}" \
21
+ -O "${sample_id}".header.sam \
22
+ $ref_opt
23
24
+awk '$0~"@PG" && $0~"ID: DRAGEN SW build" && $0~"VN: 05.021.604.3.7.8"' "${sample_id}".header.sam \
25
+ | wc -l \
26
+ > is_dragen_3_7_8.txt
0 commit comments