Skip to content

Commit

Permalink
Fix snakemake typo
Browse files Browse the repository at this point in the history
  • Loading branch information
ahmig committed Jan 4, 2024
1 parent 740fcdc commit 9b6ba7f
Showing 1 changed file with 13 additions and 10 deletions.
23 changes: 13 additions & 10 deletions workflow/rules/fetch.smk
Original file line number Diff line number Diff line change
Expand Up @@ -2,30 +2,33 @@ rule fetch_alignment_reference:
threads: 1
conda: "../envs/fetch.yaml"
params:
ref = config["ALIGNMENT_REFERENCE"]
id = config["ALIGNMENT_REFERENCE"],
db = "nucleotide",
format = "fasta"
output:
fasta = OUTDIR/"reference.fasta"
log:
LOGDIR / "fetch_alignment_reference" / "log.txt"
shell:
"esearch -db nucleotide -query {params.ref} | efetch -format fasta > {output.fasta} 2> {log}"
wrapper:
"v3.3.3/bio/entrez/efetch"


rule fetch_reference_gb:
threads: 1
conda: "../envs/fetch.yaml"
params:
ref = config["ALIGNMENT_REFERENCE"]
id = config["ALIGNMENT_REFERENCE"],
db = "gb",
format = "fasta"
output:
fasta = OUTDIR/"reference.gb"
log:
LOGDIR / "fetch_reference_gb" / "log.txt"
shell:
"esearch -db nucleotide -query {params.ref} | efetch -format gb > {output.fasta} 2> {log}"

wrapper:
"v3.3.3/bio/entrez/efetch"


rule fetch_mapping_references:
rule fetch_mapping_references: # TODO: use wrapper
threads: 1
conda: "../envs/fetch.yaml"
input:
Expand All @@ -45,13 +48,13 @@ rule fetch_mapping_references:
rule fetch_alignment_annotation:
threads: 1
params:
ref = config["ALIGNMENT_REFERENCE"]
id = config["ALIGNMENT_REFERENCE"]
output:
temp(OUTDIR/"reference.gff3")
log:
LOGDIR / "fetch_alignment_annotation" / "log.txt"
shell:
"curl 'https://www.ncbi.nlm.nih.gov/sviewer/viewer.cgi?db=nuccore&report=gff3&id={params.ref}' -o {output} -s 2>{log}"
"curl 'https://www.ncbi.nlm.nih.gov/sviewer/viewer.cgi?db=nuccore&report=gff3&id={params.id}' -o {output} -s 2>{log}"


rule fetch_problematic_vcf:
Expand Down

0 comments on commit 9b6ba7f

Please sign in to comment.