Skip to content

Commit

Permalink
modifying TMP_DIR
Browse files Browse the repository at this point in the history
  • Loading branch information
LiaOb21 committed Jan 3, 2024
1 parent cb88e67 commit 16536eb
Showing 1 changed file with 1 addition and 2 deletions.
3 changes: 1 addition & 2 deletions workflow/rules/picard.smk
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,6 @@ PICARD = r"""{CONDA_PREFIX}/share/picard-*/picard.jar"""
rule picard:
input:
tmp_bam="results/arima_mapping_pipeline/TMP_DIR/{sample}.bam",
TMP_DIR = "results/arima_mapping_pipeline/TMP_DIR"
output:
bam_paired="results/arima_mapping_pipeline/PAIR_DIR/{sample}.bam",
mark_dup="results/arima_mapping_pipeline/REP_DIR/{sample}_rep1.bam",
Expand All @@ -23,7 +22,7 @@ rule picard:
mkdir -p results/arima_mapping_pipeline/PAIR_DIR/
mkdir -p results/arima_mapping_pipeline/REP_DIR/
java -Xmx4G -Djava.io.tmpdir=temp/ -jar {PICARD} AddOrReplaceReadGroups INPUT={input.tmp_bam} OUTPUT={output.bam_paired} ID={wildcards.sample} LB={wildcards.sample} SM={wildcards.sample} PL=ILLUMINA PU=none >> {log} 2>&1
java -Xmx30G -XX:-UseGCOverheadLimit -Djava.io.tmpdir=temp/ -jar {PICARD} MarkDuplicates INPUT={output.bam_paired} OUTPUT={output.mark_dup} METRICS_FILE={output.metrics} TMP_DIR={input.TMP_DIR} ASSUME_SORTED=TRUE VALIDATION_STRINGENCY=LENIENT REMOVE_DUPLICATES=TRUE >> {log} 2>&1
java -Xmx30G -XX:-UseGCOverheadLimit -Djava.io.tmpdir=temp/ -jar {PICARD} MarkDuplicates INPUT={output.bam_paired} OUTPUT={output.mark_dup} METRICS_FILE={output.metrics} TMP_DIR=results/arima_mapping_pipeline/TMP_DIR ASSUME_SORTED=TRUE VALIDATION_STRINGENCY=LENIENT REMOVE_DUPLICATES=TRUE >> {log} 2>&1
samtools index {output.mark_dup} >> {log} 2>&1
perl scripts/get_stats.pl {output.mark_dup} > {output.stats} >> {log} 2>&1
"""

0 comments on commit 16536eb

Please sign in to comment.