Skip to content

Commit b4ec046

Browse files
committed
Account for multiple lofreq VCF naming styles
1 parent 277d7ed commit b4ec046

File tree

1 file changed

+3
-1
lines changed

1 file changed

+3
-1
lines changed

SequenceAnalysis/src/org/labkey/sequenceanalysis/run/analysis/MergeLoFreqVcfHandler.java

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -548,7 +548,9 @@ private Map<String, Integer> getContigToOffset(SAMSequenceDictionary dict)
548548

549549
private int getReadDepth(File vcf, Map<String, Integer> contigToOffset, String contig, int position1, JobContext ctx) throws PipelineJobException
550550
{
551-
File gatkDepth = new File(vcf.getParentFile(), vcf.getName().replaceAll(".all.vcf.gz", ".coverage"));
551+
//NOTE: Earlier iterations used multiple VCF naming conventions
552+
String basename = vcf.getName().split("lofreq")[0];
553+
File gatkDepth = new File(vcf.getParentFile(), basename + "lofreq.coverage");
552554
if (!gatkDepth.exists())
553555
{
554556
throw new PipelineJobException("File not found: " + gatkDepth.getPath());

0 commit comments

Comments
 (0)