Skip to content

Commit

Permalink
fixup! Use only the 21L sequence from references
Browse files Browse the repository at this point in the history
  • Loading branch information
victorlin committed Sep 30, 2024
1 parent 97a3309 commit 229eebb
Showing 1 changed file with 5 additions and 1 deletion.
6 changes: 5 additions & 1 deletion nextstrain_profiles/nextstrain-gisaid-21L/prefilter.smk
Original file line number Diff line number Diff line change
Expand Up @@ -91,8 +91,10 @@ rule gisaid_21L_aligned:
r"""
exec 2> {log:q}
echo "21L" > reference_ids.txt
< {input.references:q} \
seqkit grep --by-name -f <(echo "21L") \
seqkit grep --by-name -f reference_ids.txt \
| zstd \
> {output.aligned}
Expand All @@ -101,4 +103,6 @@ rule gisaid_21L_aligned:
| seqkit grep --by-name -f {input.strains:q} \
| zstd -T$(({threads} - 2)) \
>> {output.aligned:q}
rm reference_ids.txt
"""

0 comments on commit 229eebb

Please sign in to comment.