Skip to content

Commit 7f721ea

Browse files
add -w flag to grep for varID searches in 04b (#139)
1 parent 2166cd2 commit 7f721ea

File tree

2 files changed

+3
-3
lines changed

2 files changed

+3
-3
lines changed

wdl/CombineReassess.wdl

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -64,7 +64,7 @@ task Vcf2Bed {
6464
svtk vcf2bed ~{vcf} ~{vcf}.bed
6565
awk '{if($6!="")print $0}' ~{vcf}.bed >~{vcf}.nonempty.bed
6666
cut -f 4 ~{regeno_file} >regeno_variants.txt
67-
fgrep -f regeno_variants.txt ~{vcf}.nonempty.bed> nonempty.txt
67+
fgrep -wf regeno_variants.txt ~{vcf}.nonempty.bed> nonempty.txt
6868
>>>
6969
output {
7070
File nonempty="nonempty.txt"

wdl/Module04b.wdl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -741,9 +741,9 @@ task ConcatRegenotypedVcfs {
741741
command <<<
742742
set -euo pipefail
743743
zcat ~{regeno_vcf} |fgrep "#" > head.txt
744-
zcat ~{regeno_vcf} |fgrep -f ~{regeno_variants} >body.txt
744+
zcat ~{regeno_vcf} |fgrep -wf ~{regeno_variants} >body.txt
745745
cat head.txt body.txt|bgzip -c > regeno.vcf.gz
746-
zcat ~{depth_vcf} |fgrep -f ~{regeno_variants} -v |bgzip -c > no_variant.vcf.gz
746+
zcat ~{depth_vcf} |fgrep -wf ~{regeno_variants} -v |bgzip -c > no_variant.vcf.gz
747747
vcf-concat regeno.vcf.gz no_variant.vcf.gz \
748748
| vcf-sort -c \
749749
| bgzip -c > ~{batch}.depth.regeno_final.vcf.gz

0 commit comments

Comments
 (0)