Skip to content
Snippets Groups Projects
Commit 0b5ebc9d authored by Kedi Cao's avatar Kedi Cao
Browse files

re-arrange rules

parent 0edb88a1
No related merge requests found
import os
rule sort_vcf:
input:
RESULTS_DIR + "/place-finish.done",
ins=RESULTS_DIR + "/insertions_unsorted.vcf"
output:
temp(RESULTS_DIR + "/insertions.vcf")
conda:
os.path.join(WORKFLOW_PATH,"snakemodules/envs/bcftools.yml")
resources:
mem_mb = resources["standard_2G"]["mem"],
runtime = resources["standard_2G"]["time"]
threads:
threads["single"]
log:
out="logs/genotype/sort_vcf.out",
err="logs/genotype/sort_vcf.err"
benchmark:
"benchmarks/genotype/sort_vcf.txt"
shell:
"bcftools sort {input.ins} "
" --output-file {output}"
" --output-type v"
" > {log.out} 2> {log.err}"
rule popins2_genotype:
input:
......
......@@ -156,4 +156,29 @@ rule popins2_place_finish:
" --prefix {WORK_DIR} "
" --insertions {rules.popins2_place_refalign.output.vcf}"
" --reference {REFERENCE} "
" > {log.out} 2> {log.err}"
rule sort_vcf:
input:
RESULTS_DIR + "/place-finish.done",
ins=RESULTS_DIR + "/insertions_unsorted.vcf"
output:
temp(RESULTS_DIR + "/insertions.vcf")
conda:
os.path.join(WORKFLOW_PATH,"snakemodules/envs/bcftools.yml")
resources:
mem_mb = resources["standard_2G"]["mem"],
runtime = resources["standard_2G"]["time"]
threads:
threads["single"]
log:
out="logs/genotype/sort_vcf.out",
err="logs/genotype/sort_vcf.err"
benchmark:
"benchmarks/genotype/sort_vcf.txt"
shell:
"bcftools sort {input.ins} "
" --output-file {output}"
" --output-type v"
" > {log.out} 2> {log.err}"
\ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment