Skip to content

Commit

Permalink
let's try removing a call to Picard SortSam that seems redundant
Browse files Browse the repository at this point in the history
  • Loading branch information
tomkinsc committed Jan 25, 2018
1 parent 947b5c1 commit 74214cb
Showing 1 changed file with 2 additions and 9 deletions.
11 changes: 2 additions & 9 deletions tools/bwa.py
Original file line number Diff line number Diff line change
Expand Up @@ -8,6 +8,7 @@
import os
import os.path
import subprocess
import shutil

import tools
import tools.samtools
Expand Down Expand Up @@ -194,15 +195,7 @@ def align_mem_one_rg(self, inBam, refDb, outBam, rgid=None, options=None,
# os.unlink(headerFile)
# os.system("samtools view -h {} > /Users/tomkinsc/Desktop/test_reheader.bam".format(reheadered_bam))

# sort
sorter = tools.picard.SortSamTool()
sorter.execute(
tmp_bam_aligned,
outBam,
sort_order='coordinate',
picardOptions=['CREATE_INDEX=true', 'VALIDATION_STRINGENCY=SILENT'],
JVMmemory=JVMmemory
)
shutil.move(tmp_bam_aligned, outBam)
#os.unlink(reheadered_bam)

def mem(self, inReads, refDb, outAlign, options=None, min_score_to_filter=None,
Expand Down

0 comments on commit 74214cb

Please sign in to comment.