Skip to content

Instantly share code, notes, and snippets.

@mvelinder
Created October 27, 2020 17:17
Show Gist options
  • Save mvelinder/11e3ff21e226a9c00f6e7757e220aac0 to your computer and use it in GitHub Desktop.
Save mvelinder/11e3ff21e226a9c00f6e7757e220aac0 to your computer and use it in GitHub Desktop.
bcftools csq and vcfanno to avoid vep
set -oexu pipefail
VCF=$1
bcftools csq \
-f /scratch/ucgd/lustre/work/u0691312/reference/ucgd_reference/GRCh38/human_g1k_v38_decoy_phix.fasta \
-g /scratch/ucgd/lustre/work/u0691312/reference/Homo_sapiens.GRCh38.100.gff3.gz \
$VCF \
| vcfanno -p 12 \
/scratch/ucgd/lustre/work/u0691312/reference/vcfanno_revel.toml \
/dev/stdin
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment