Skip to content

Commit 18013c8

Browse files
committed
Bugfix to MergeBamAlignment
1 parent e90c143 commit 18013c8

File tree

1 file changed

+1
-3
lines changed

1 file changed

+1
-3
lines changed

SequenceAnalysis/src/org/labkey/sequenceanalysis/run/util/MergeBamAlignmentWrapper.java

Lines changed: 1 addition & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -94,7 +94,7 @@ else if (header.getReadGroups().size() > 1)
9494
FastqToSamWrapper fq = new FastqToSamWrapper(getLogger());
9595
fq.setOutputDir(alignedBam.getParentFile());
9696
fq.setStringency(ValidationStringency.SILENT);
97-
File unmappedReadsBam = fq.execute(pair.getKey(), pair.getValue(), SAMFileHeader.SortOrder.queryname, rg);
97+
File unmappedReadsBam = fq.execute(pair.getKey(), pair.getValue(), SAMFileHeader.SortOrder.unsorted, rg);
9898
if (!unmappedReadsBam.exists())
9999
{
100100
throw new PipelineJobException("BAM file not created, expected: " + unmappedReadsBam.getPath());
@@ -121,8 +121,6 @@ else if (header.getReadGroups().size() > 1)
121121
params.add("--UNMAPPED_BAM");
122122
params.add(finalUnmappedReadsBam.getPath());
123123

124-
//TODO: bisulfiteSequence
125-
126124
params.add("--REFERENCE_SEQUENCE");
127125
params.add(refFasta.getPath());
128126

0 commit comments

Comments
 (0)