Skip to content

Commit

Permalink
Consistent naming of reference genome files
Browse files Browse the repository at this point in the history
* Dropped "_dengue" from names, since it's redundant in the context of the project
* Added "_genome" to the end of the reference genome files, to parallel the "_E"
  at the end of the reference E sequence files.
  • Loading branch information
j23414 committed May 2, 2024
1 parent 2f48ddb commit ab5dff4
Show file tree
Hide file tree
Showing 8 changed files with 6 additions and 6 deletions.
File renamed without changes.
2 changes: 1 addition & 1 deletion phylogenetic/rules/annotate_phylogeny.smk
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +45,7 @@ rule translate:
tree = "results/tree_{serotype}_{gene}.nwk",
node_data = "results/nt-muts_{serotype}_{gene}.json",
# The genbank references for the E gene are dynamically generated files located within the results folder.
reference = lambda wildcard: "config/reference_dengue_{serotype}.gb" if wildcard.gene in ["genome"] else "results/config/reference_dengue_{serotype}_{gene}.gb"
reference = lambda wildcard: "config/reference_{serotype}_{gene}.gb" if wildcard.gene in ["genome"] else "results/config/reference_{serotype}_{gene}.gb"
output:
node_data = "results/aa-muts_{serotype}_{gene}.json"
shell:
Expand Down
2 changes: 1 addition & 1 deletion phylogenetic/rules/prepare_sequences.smk
Original file line number Diff line number Diff line change
Expand Up @@ -82,7 +82,7 @@ rule align:
"""
input:
sequences = "results/filtered_{serotype}.fasta",
reference = "config/reference_dengue_{serotype}.gb"
reference = "config/reference_{serotype}_genome.gb"
output:
alignment = "results/aligned_{serotype}_genome.fasta"
shell:
Expand Down
8 changes: 4 additions & 4 deletions phylogenetic/rules/prepare_sequences_E.smk
Original file line number Diff line number Diff line change
Expand Up @@ -22,10 +22,10 @@ rule generate_E_reference_files:
Generating reference files for the E gene
"""
input:
reference = "config/reference_dengue_{serotype}.gb",
reference = "config/reference_{serotype}_genome.gb",
output:
fasta = "results/config/reference_dengue_{serotype}_E.fasta",
genbank = "results/config/reference_dengue_{serotype}_E.gb",
fasta = "results/config/reference_{serotype}_E.fasta",
genbank = "results/config/reference_{serotype}_E.gb",
params:
gene = "E",
shell:
Expand All @@ -43,7 +43,7 @@ rule nextclade3_cut_E:
"""
input:
sequences = "data/sequences_{serotype}.fasta",
reference = "results/config/reference_dengue_{serotype}_E.fasta"
reference = "results/config/reference_{serotype}_E.fasta"
output:
sequences = "results/sequences_{serotype}_E.fasta"
params:
Expand Down

0 comments on commit ab5dff4

Please sign in to comment.