Bugfix to DeepVariant docker command

BimberLab · Jan 22, 2024 · c522d20 · c522d20
1 parent 207b83e
commit c522d20
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/SequenceAnalysis/src/org/labkey/sequenceanalysis/run/analysis/DeepVariantAnalysis.java b/SequenceAnalysis/src/org/labkey/sequenceanalysis/run/analysis/DeepVariantAnalysis.java
@@ -204,6 +204,9 @@ private File ensureLocalCopy(File input, File workingDirectory, PipelineOutputTr
         public void execute(File inputBam, File refFasta, File outputGvcf, PipelineOutputTracker tracker, String binVersion, List<String> extraArgs) throws PipelineJobException
         {
             File workDir = outputGvcf.getParentFile();
+            File outputVcf = new File(outputGvcf.getPath().replaceAll(".g.vcf", ".vcf"));
+            tracker.addIntermediateFile(outputVcf);
+            tracker.addIntermediateFile(new File(outputVcf.getPath() + ".tbi"));
 
             File inputBamLocal = ensureLocalCopy(inputBam, workDir, tracker);
             ensureLocalCopy(SequenceUtil.getExpectedIndex(inputBam), workDir, tracker);
@@ -221,6 +224,7 @@ public void execute(File inputBam, File refFasta, File outputGvcf, PipelineOutpu
             bashArgs.add("--ref=/work/" + refFastaLocal.getName());
             bashArgs.add("--reads=/work/" + inputBamLocal.getName());
             bashArgs.add("--output_gvcf=/work/" + outputGvcf.getName());
+            bashArgs.add("--output_vcf=/work/" + outputVcf.getName());
             Integer maxThreads = SequencePipelineService.get().getMaxThreads(getLogger());
             if (maxThreads != null)
             {