Skip to content
Snippets Groups Projects
Commit 89eba335 authored by Floreal Cabanettes's avatar Floreal Cabanettes
Browse files

Move windowed files in subdir for pindel (try fixing problem...)

parent a2da3720
No related branches found
No related tags found
1 merge request!12Split pindel detection in 20mb windows
......@@ -71,7 +71,7 @@ class SnakemakeUtils:
def get_pindel_chr_batches(self, wildcards):
inputs = []
for batch in self.chr_batches[wildcards.chrom]:
inputs.append("{batch}/pindel/pindel_{chrom}_{chrbatch}_{svtype}.gz".format(
inputs.append("{batch}/pindel/{chrbatch}/pindel_{chrom}_{svtype}.gz".format(
batch=wildcards.batch, chrom=wildcards.chrom, chrbatch="-".join(map(str, batch)),
svtype=wildcards.svtype
))
......
......@@ -30,9 +30,9 @@ rule pindel:
genome = REFERENCE,
fai = REFERENCE + ".fai"
output:
"{batch}/pindel/pindel_{chrom}_{chrbatch}_D.gz",
"{batch}/pindel/pindel_{chrom}_{chrbatch}_INV.gz",
"{batch}/pindel/pindel_{chrom}_{chrbatch}_TD.gz",
"{batch}/pindel/{chrbatch}/pindel_{chrom}_D.gz",
"{batch}/pindel/{chrbatch}/pindel_{chrom}_INV.gz",
"{batch}/pindel/{chrbatch}/pindel_{chrom}_TD.gz",
threads:
get_threads("pindel", 4)
log:
......@@ -50,7 +50,7 @@ rule pindel:
"""
for svtype in D TD INV RP SI LI BP CloseEndMapped INT_final;
do
gzip -f {wildcards.batch}/pindel/pindel_{wildcards.chrom}_{wildcards.chrbatch}_${{svtype}}
gzip -f {wildcards.batch}/pindel/{wildcards.chrbatch}/pindel_{wildcards.chrom}_${{svtype}}
done
"""
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment