Skip to content

Commit e2b2360

Browse files
committed
Ensure BAM index present
1 parent 9b02414 commit e2b2360

File tree

1 file changed

+4
-0
lines changed

1 file changed

+4
-0
lines changed

SequenceAnalysis/src/org/labkey/sequenceanalysis/run/preprocessing/FilterReadsStep.java

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -14,11 +14,13 @@
1414
import org.labkey.api.sequenceanalysis.pipeline.PreprocessingStep;
1515
import org.labkey.api.sequenceanalysis.pipeline.ReferenceGenome;
1616
import org.labkey.api.sequenceanalysis.pipeline.SequenceAnalysisJobSupport;
17+
import org.labkey.api.sequenceanalysis.pipeline.SequencePipelineService;
1718
import org.labkey.api.sequenceanalysis.pipeline.ToolParameterDescriptor;
1819
import org.labkey.api.util.FileUtil;
1920
import org.labkey.api.util.Pair;
2021
import org.labkey.sequenceanalysis.run.alignment.BWAMemWrapper;
2122
import org.labkey.sequenceanalysis.run.analysis.UnmappedReadExportHandler;
23+
import org.labkey.sequenceanalysis.util.SequenceUtil;
2224

2325
import java.io.File;
2426
import java.util.ArrayList;
@@ -75,6 +77,8 @@ public Output processInputFile(File inputFile, @Nullable File inputFile2, File o
7577
AlignmentOutputImpl alignmentOutput = new AlignmentOutputImpl();
7678
String basename = SequenceAnalysisService.get().getUnzippedBaseName(inputFile.getName()) + ".filterAlign";
7779
wrapper.performMemAlignment(getPipelineCtx().getJob(), alignmentOutput, inputFile, inputFile2, outputDir, genome, basename, bwaArgs);
80+
81+
SequencePipelineService.get().ensureBamIndex(inputFile, getPipelineCtx().getLogger(), true);
7882

7983
output.addIntermediateFile(alignmentOutput.getBAM());
8084
output.addIntermediateFile(new File(alignmentOutput.getBAM().getPath() + ".bai"));

0 commit comments

Comments
 (0)