cwl/gvcf2fasta/src/fixvcf-get_bed_varonlyvcf.sh

   1 # Copyright (C) The Lightning Authors. All rights reserved.
   2 #
   3 # SPDX-License-Identifier: AGPL-3.0
   4
   5 #!/bin/bash
   6
   7 set -eo pipefail
   8
   9 sampleid="$1"
  10 vcf="$2"
  11 gqcutoff="$3"
  12 genomebed="$4"
  13
  14 bcftools view --trim-alt-alleles $vcf | egrep -v "\*|<NON_REF>" | tee \
  15   >( /gvcf_regions/gvcf_regions.py --min_GQ $gqcutoff - > "$sampleid".bed ) \
  16   >( rtg vcffilter -i - -o - --remove-overlapping | awk '{if ($5 != ".") print $0}' | bgzip -c > "$sampleid"_varonly.vcf.gz ) \
  17   > /dev/null
  18
  19 bedtools subtract -a $genomebed -b "$sampleid".bed > "$sampleid"_nocall.bed
  20 rm "$sampleid".bed
  21 tabix "$sampleid"_varonly.vcf.gz