Skip to content

Commit 3ef3d58

Browse files
committed
Update cellranger arguments
1 parent d8955d0 commit 3ef3d58

File tree

2 files changed

+6
-4
lines changed

2 files changed

+6
-4
lines changed

singlecell/src/org/labkey/singlecell/run/CellRangerFeatureBarcodeHandler.java

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -252,6 +252,8 @@ private void processType(List<Pair<File, File>> inputFastqs, AlignmentOutputImpl
252252
File indexDir = makeDummyIndex(ctx);
253253
args.add("--transcriptome=" + indexDir.getPath());
254254

255+
args.add("--create-bam=false");
256+
255257
wrapper.setWorkingDir(ctx.getOutputDir());
256258

257259
//Note: we can safely assume only this server is working on these files, so if the _lock file exists, it was from a previous failed job.
@@ -284,12 +286,10 @@ private void processType(List<Pair<File, File>> inputFastqs, AlignmentOutputImpl
284286
File outsdir = new File(crDir, "outs");
285287

286288
File bam = new File(outsdir, "possorted_genome_bam.bam");
287-
if (!bam.exists())
289+
if (bam.exists())
288290
{
289-
throw new PipelineJobException("Unable to find file: " + bam.getPath());
291+
throw new PipelineJobException("Did not expect to find BAM: " + bam.getPath());
290292
}
291-
bam.delete();
292-
SequenceAnalysisService.get().getExpectedBamOrCramIndex(bam).delete();
293293

294294
wrapper.deleteSymlinks(wrapper.getLocalFastqDir(ctx.getOutputDir()));
295295

singlecell/src/org/labkey/singlecell/run/CellRangerGexCountStep.java

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -352,6 +352,8 @@ public AlignmentOutput performAlignment(Readset rs, List<File> inputFastqs1, @Nu
352352
File indexDir = AlignerIndexUtil.getIndexDir(referenceGenome, getIndexCachedDirName(getPipelineCtx().getJob()));
353353
args.add("--transcriptome=" + indexDir.getPath());
354354

355+
args.add("--create-bam=true");
356+
355357
getWrapper().setWorkingDir(outputDirectory);
356358

357359
//Note: we can safely assume only this server is working on these files, so if the _lock file exists, it was from a previous failed job.

0 commit comments

Comments
 (0)