Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

fix: processing vembrane config #330

Merged
merged 9 commits into from
Oct 17, 2024
Merged
Show file tree
Hide file tree
Changes from 1 commit
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion workflow/rules/candidate_calling.smk
Original file line number Diff line number Diff line change
Expand Up @@ -20,7 +20,7 @@ rule freebayes:
),
threads: max(workflow.cores - 1, 1) # use all available cores -1 (because of the pipe) for calling
wrapper:
"v2.7.0/bio/freebayes"
"v4.4.0/bio/freebayes"
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is this necessary for anything (like a bug fix in a newer version of freebayes or the wrapper)? Because otherwise I would leave this as is, so as not to trigger rerunning of the candidate calling in workflows that update to a newer release after merging this PR...



rule delly:
Expand Down
10 changes: 6 additions & 4 deletions workflow/rules/common.smk
Original file line number Diff line number Diff line change
Expand Up @@ -1169,10 +1169,12 @@ def get_vembrane_config(wildcards, input):

def append_items(items, renames, field_func, header_func=None):
for item in items:
if item in renames:
header_name = renames[item]["name"]
if "expr" in renames[item]:
parts_field = renames[item]["expr"]
if type(item) == dict or item in renames:
entry = item if type(item) == dict else renames[item]
header_name = entry["name"]
if "expr" in entry:
parts_field = entry["expr"]
# TODO How can this case be covered in the config file?
FelixMoelder marked this conversation as resolved.
Show resolved Hide resolved
else:
parts_field = field_func(item)
else:
Expand Down
2 changes: 1 addition & 1 deletion workflow/rules/mapping.smk
Original file line number Diff line number Diff line change
Expand Up @@ -83,7 +83,7 @@ rule mark_duplicates:
#https://broadinstitute.github.io/picard/faq.html
mem_mb=3000,
wrapper:
"v2.5.0/bio/picard/markduplicates"
"v4.2.0/bio/picard/markduplicates"
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

As above: if such an update is not really necessary for anything in this pull request, I would defer updating to whenever this is really necessary. Otherwise this will trigger reruns in previously workflows where we update to a newer workflow version...



rule calc_consensus_reads:
Expand Down
2 changes: 1 addition & 1 deletion workflow/rules/plugins.smk
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@ rule download_revel:
conda:
"../envs/curl.yaml"
shell:
"curl https://rothsj06.dmz.hpc.mssm.edu/revel-v1.3_all_chromosomes.zip -o {output} &> {log}"
"curl https://zenodo.org/records/7072866/files/revel-v1.3_all_chromosomes.zip -o {output} &> {log}"
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Nice, let's hope this stays more robust in the long-term and that newer versions will also end up on Zenodo!



rule process_revel_scores:
Expand Down
Loading