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

Add E gene trees #18

Closed
wants to merge 14 commits into from
Prev Previous commit
Next Next commit
Rename rule to describe what it does instead of the internal tool it …
…uses.

The purpose of this rule is to align sequences to an E gene reference sequence
and extract the E gene region, if any. The rule is renamed to reflect this.
j23414 committed May 2, 2024
commit 6a5cd62fc2bdc3b0334326fc41b7b2ba4ccb60aa
4 changes: 2 additions & 2 deletions phylogenetic/rules/prepare_sequences_E.smk
Original file line number Diff line number Diff line change
@@ -14,7 +14,7 @@ This part of the workflow usually includes the following steps:
See Augur's usage docs for these commands for more details.
"""

ruleorder: nextclade3_cut_E > decompress
ruleorder: align_and_extract_E > decompress
ruleorder: filter_E > align

rule generate_E_reference_files:
@@ -37,7 +37,7 @@ rule generate_E_reference_files:
--gene {params.gene}
"""

rule nextclade3_cut_E:
rule align_and_extract_E:
"""
Cutting sequences to the length of the E gene reference sequence
"""