broadinstitute / gatk-sv

A structural variation pipeline for short-read sequencing
BSD 3-Clause "New" or "Revised" License
160 stars 71 forks source link

Make SplitVariants task in TasksGenotypeBatch.wdl compatible with downstream analysis #650

Closed kirtanav98 closed 4 months ago

kirtanav98 commented 4 months ago

The SplitVariants task used to have some lines to switch columns 5 and 6 of the bed file output, which is read in downstream tasks of TrainRDGenotyping.GenotypePESR. This causes the TrainRDGenotyping.GenotypePESR to error out reporting.

Error: WARNING: Incorrect CNV type specified 1: stop("WARNING: Incorrect CNV type specified") The python script splitvariants.py was modified to switch the columns to the appropriate order to be compatible with downstream analysis requirements. The script passed testing with both womtool and cromshell on the 1kg reference panel with the bca flag set to both True and False.