Skip to content

Commit 168ab55

Browse files
committed
Conditionalize GenotypeGVCFs and ExcessHet/InbreedingCoef
1 parent 519e0f6 commit 168ab55

File tree

2 files changed

+13
-4
lines changed

2 files changed

+13
-4
lines changed

SequenceAnalysis/resources/web/SequenceAnalysis/panel/VariantProcessingPanel.js

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -268,6 +268,12 @@ Ext4.define('SequenceAnalysis.panel.VariantProcessingPanel', {
268268
label: 'Intervals to Exclude',
269269
description: 'A BED or similar file with intervals to skip',
270270
defaultValue: null
271+
},{
272+
fieldXtype: 'ldk-expdatafield',
273+
name: 'skipExcessHetAndInbreeding',
274+
label: 'Skip Excess Het And Inbreeding',
275+
description: 'If checked, the ExcessHet and InbreedingCoeff annotations will be skipped, which can be important when using force-output-intervals',
276+
defaultValue: null
271277
},{
272278
fieldXtype: 'checkbox',
273279
name: 'doCopyInputs',

SequenceAnalysis/src/org/labkey/sequenceanalysis/analysis/GenotypeGVCFHandler.java

Lines changed: 7 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -636,11 +636,14 @@ private void processOneInput(JobContext ctx, PipelineJob job, ReferenceGenome ge
636636
}
637637

638638
// NOTE: added to side-step https://github.com/broadinstitute/gatk/issues/7938
639-
toolParams.add("-AX");
640-
toolParams.add("InbreedingCoeff");
639+
if (ctx.getParams().optBoolean("variantCalling.GenotypeGVCFs.skipExcessHetAndInbreeding"))
640+
{
641+
toolParams.add("-AX");
642+
toolParams.add("InbreedingCoeff");
641643

642-
toolParams.add("-AX");
643-
toolParams.add("ExcessHet");
644+
toolParams.add("-AX");
645+
toolParams.add("ExcessHet");
646+
}
644647

645648
if (!ctx.getParams().isNull("variantCalling.GenotypeGVCFs.maxGenotypeCount"))
646649
{

0 commit comments

Comments
 (0)