From 8328f1ce0952a9717a86ce8de3172edef96f501e Mon Sep 17 00:00:00 2001 From: naumenko-sa Date: Wed, 24 Oct 2018 14:28:44 -0400 Subject: [PATCH] modified: cre.vt.normalize.sh --- cre.vt.normalize.sh | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/cre.vt.normalize.sh b/cre.vt.normalize.sh index f2dd2d2..c9b3509 100755 --- a/cre.vt.normalize.sh +++ b/cre.vt.normalize.sh @@ -5,7 +5,9 @@ #PBS -d . #PBS -l vmem=20g,mem=20g +# do not fail if a variant in inconsistent with the reference + bname=`basename $vcf .vcf.gz` -vt normalize -r ~/work/tools/bcbio/genomes/Hsapiens/GRCh37/seq/GRCh37.fa $vcf | vt uniq - -o $bname.normalized.vcf.gz +vt normalize -r ~/work/tools/bcbio/genomes/Hsapiens/GRCh37/seq/GRCh37.fa $vcf -n | vt uniq - -o $bname.normalized.vcf.gz tabix $bname.normalized.vcf.gz