Skip to content

Commit

Permalink
fixing cleavage site export to include wildcards
Browse files Browse the repository at this point in the history
  • Loading branch information
lhmoncla committed Apr 22, 2024
1 parent c49732c commit 94b2641
Showing 1 changed file with 4 additions and 8 deletions.
12 changes: 4 additions & 8 deletions Snakefile
Original file line number Diff line number Diff line change
Expand Up @@ -8,11 +8,6 @@ path_to_fauna = '../fauna'
rule all:
input:
auspice_json = expand("auspice/flu_avian_{subtype}_{segment}_{time}.json", subtype=SUBTYPES, segment=SEGMENTS,time=TIME)
# auspice_json = expand("auspice/flu_avian_h5nx_{segment}_{time}.json", segment=SEGMENTS,time=TIME),
# auspice_json = expand("auspice/flu_avian_h5n1_{segment}_{time}.json", segment=SEGMENTS,time=TIME),
# auspice_json = expand("auspice/flu_avian_h7n9_{segment}_all-time.json", segment=SEGMENTS),
# auspice_json = expand("auspice/flu_avian_h9n2_{segment}_all-time.json", segment=SEGMENTS)


#auspice_json = expand("auspice/avian-flu_{subtype}_{segment}.json", subtype=SUBTYPES, segment=SEGMENTS)
#sequences = expand("results/sequences_{subtype}_{segment}.fasta", subtype=SUBTYPES, segment=SEGMENTS),
Expand All @@ -31,6 +26,7 @@ rule files:

files = rules.files.params


def download_by(w):
db = {'h5nx': 'subtype:h5n1,h5n2,h5n3,h5n4,h5n5,h5n6,h5n7,h5n8,h5n9', 'h5n1': 'subtype:h5n1', 'h7n9': 'subtype:h7n9', 'h9n2': 'subtype:h9n2'}
return(db[w.subtype])
Expand Down Expand Up @@ -305,10 +301,10 @@ rule traits:
rule cleavage_site:
message: "determining sequences that harbor furin cleavage sites"
input:
alignment = "results/aligned_{subtype}_ha.fasta"
alignment = "results/aligned_{subtype}_ha_{time}.fasta"
output:
cleavage_site_annotations = "results/cleavage-site_{subtype}_ha.json",
cleavage_site_sequences = "results/cleavage-site-sequences_{subtype}_ha.json"
cleavage_site_annotations = "results/cleavage-site_{subtype}_ha_{time}.json",
cleavage_site_sequences = "results/cleavage-site-sequences_{subtype}_ha_{time}.json"
shell:
"""
python scripts/annotate-ha-cleavage-site.py \
Expand Down

0 comments on commit 94b2641

Please sign in to comment.