Skip to content

Commit

Permalink
Merge pull request #19 from Akazhiel/master
Browse files Browse the repository at this point in the history
Fixed bug in gatk
  • Loading branch information
jfnavarro authored Feb 8, 2021
2 parents 676aa8f + 30af1a0 commit a9a89a6
Showing 1 changed file with 4 additions and 8 deletions.
12 changes: 4 additions & 8 deletions dna_pipeline.py
Original file line number Diff line number Diff line change
Expand Up @@ -131,10 +131,6 @@ def main(R1_NORMAL,

cmd2 = '{} --java-options "-Xmx32g" MarkDuplicatesSpark -I sample2_header.bam -O sample2_dedup.bam'.format(GATK)

# Wait for the processes to finish in parallel
p1.wait()
p2.wait()

else:
cmd1 = '{} --java-options "-Xmx32g" MarkDuplicates -I sample1_header.bam -O sample1_dedup.bam ' \
'--CREATE_INDEX true -M sample1_dup_metrics.txt'.format(GATK)
Expand All @@ -143,10 +139,10 @@ def main(R1_NORMAL,
'--CREATE_INDEX true -M sample2_dup_metrics.txt'.format(GATK)

# Wait for the processes to finish in parallel
p1 = exec_command(cmd1, detach=True)
p2 = exec_command(cmd2, detach=True)
p1.wait()
p2.wait()
p1 = exec_command(cmd1, detach=True)
p2 = exec_command(cmd2, detach=True)
p1.wait()
p2.wait()

intervals_cmd = '--intervals {}'.format(INTERVALS) if INTERVALS else ''

Expand Down

0 comments on commit a9a89a6

Please sign in to comment.