Skip to content

Commit ee6ea2c

Browse files
committed
Make CellRangerVDJ merge both filtered and all contig files
1 parent 350e9a2 commit ee6ea2c

File tree

1 file changed

+6
-5
lines changed

1 file changed

+6
-5
lines changed

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

Lines changed: 6 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -389,7 +389,10 @@ public AlignmentStep.AlignmentOutput performAlignment(Readset rs, List<File> inp
389389
File csvAB = processOutputsForType(id, rs, referenceGenome, outdir, output, "vdj_t");
390390
File csvGD = processOutputsForType(id, rs, referenceGenome, outdir, output, "vdj_t_gd");
391391

392-
File combinedCSV = processAndMergeCSVs(csvAB, csvGD, getPipelineCtx().getLogger());
392+
processAndMergeCSVs(csvAB, csvGD, getPipelineCtx().getLogger());
393+
394+
// Repeat for filtered contig annotations:
395+
processAndMergeCSVs(new File(csvAB.getParentFile(), "filtered_contig_annotations.csv"), new File(csvGD.getParentFile(), "filtered_contig_annotations.csv"), getPipelineCtx().getLogger());
393396

394397
//NOTE: this folder has many unnecessary files and symlinks that get corrupted when we rename the main outputs
395398
File directory = new File(outdir.getParentFile(), "SC_MULTI_CS");
@@ -788,9 +791,9 @@ protected File getExe()
788791
return SequencePipelineService.get().getExeForPackage("CELLRANGERPATH", "cellranger");
789792
}
790793

791-
private static File processAndMergeCSVs(File abCSV, File gdCSV, Logger log) throws PipelineJobException
794+
private static void processAndMergeCSVs(File abCSV, File gdCSV, Logger log) throws PipelineJobException
792795
{
793-
File output = new File(abCSV.getParentFile(), "all_contig_annotations_combined.csv");
796+
File output = new File(abCSV.getParentFile(), FileUtil.getBaseName(abCSV) + "_combined.csv");
794797

795798
try (PrintWriter writer = PrintWriters.getPrintWriter(output))
796799
{
@@ -801,8 +804,6 @@ private static File processAndMergeCSVs(File abCSV, File gdCSV, Logger log) thro
801804
{
802805
throw new PipelineJobException(e);
803806
}
804-
805-
return output;
806807
}
807808

808809
private static void processCSV(PrintWriter writer, boolean printHeader, File inputCsv, Logger log, List<String> acceptableChains, String chainType) throws IOException

0 commit comments

Comments
 (0)