diff --git a/cluster.yaml b/cluster.yaml index ebcfc934aa58a85618537322b9212d0a6fce65da..43a7f9684bb46ee56f6476e68c411650e79ee6ac 100644 --- a/cluster.yaml +++ b/cluster.yaml @@ -18,7 +18,7 @@ genomestrip: mem: 8 preprocess: - mem: 8 + mem: 16 time: "72:00:00" cnvpipeline: diff --git a/snakecnv/bin/gstrip_preprocess.py b/snakecnv/bin/gstrip_preprocess.py index 5e121351906e8ddec32649a5c618f582a0be3bcd..c65492e9a797c1f1dd3dd5763e83031edcbafb76 100755 --- a/snakecnv/bin/gstrip_preprocess.py +++ b/snakecnv/bin/gstrip_preprocess.py @@ -40,7 +40,7 @@ def launch(bamlist, reference, gendermap, outdir, outprefix): old_libpath = ":" + os.environ["LD_LIBRARY_PATH"] if "LD_LIBRARY_PATH" in os.environ else "" os.environ["LD_LIBRARY_PATH"] = sv_dir + "/bwa" + old_libpath - mx = "-Xmx4g" + mx = "-Xmx16g" classpath = "{sv_dir}/lib/SVToolkit.jar:{sv_dir}/lib/gatk/GenomeAnalysisTK.jar:{sv_dir}/lib/gatk/Queue.jar". \ format(sv_dir=sv_dir) @@ -58,7 +58,7 @@ def launch(bamlist, reference, gendermap, outdir, outprefix): .format(classpath=classpath, mx=mx, sv_dir=sv_dir), shell=True).returncode if exit_code != 0: return fail("An error has occurred. See above.") - + command = ( """LC_ALL=C java -cp {classpath} {mx} \ org.broadinstitute.gatk.queue.QCommandLine \ diff --git a/snakecnv/tools/pindel.snk b/snakecnv/tools/pindel.snk index d290f7407c7379d1ca23931f7bdbec78a0badb6e..acb74cda579dc12cec0ca09b080f8dfbf8078e48 100644 --- a/snakecnv/tools/pindel.snk +++ b/snakecnv/tools/pindel.snk @@ -41,10 +41,11 @@ rule pindel: shell: "pindel -f {input.genome} -i {input.statfile}" " --min_perfect_match_around_BP 20 " + " --max_range_index 3 " " --exclude {input.excluded}" " --number_of_threads {threads}" " --report_interchromosomal_events false" - " --window_size 1" # 1Mb + " --window_size 1" # 1Mb " --minimum_support_for_event 3 -c {wildcards.chrom}:{wildcards.chrbatch}" " -o {wildcards.batch}/pindel/{wildcards.chrom}/{wildcards.chrbatch}/pindel_{wildcards.chrom}" " 1>{log.stdout} 2>{log.stderr}"