mtandon09 / CCBR_GATK4_Exome_Seq_Pipeline

An easy-to-use, flexible variant calling pipeline for use on the Biowulf cluster at NIH
https://mtandon09.github.io/CCBR_GATK4_Exome_Seq_Pipeline/
MIT License
4 stars 3 forks source link

Patch: BCF tools in SOB Detection rule #29

Closed skchronicles closed 2 years ago

skchronicles commented 2 years ago

The following BCF tool sub commands (annotate and filter) work fine when given a empty VCF file (i.e. an input file only containing a VCF header). With that being said, there is no need to reset the $bcf_annotate_option an empty string. This only causes issues later on when BCF tools is parsing command line arguments. It interprets an empty string as a file with its input name being an empty string: image

As so, BCF tools will work with bcf_annotate_option="-e 'INFO/pArtifact < 0.05' " when the VCF file is empty so there is no need to reset the variable.