Skip to content

Commit

Permalink
fix: correctly parse amino acids that have a name with a dot in it
Browse files Browse the repository at this point in the history
  • Loading branch information
florianzwagemaker committed Nov 28, 2022
1 parent cd7ee96 commit 3e14dd5
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion ViroConstrictor/workflow/scripts/group_aminoacids.py
Original file line number Diff line number Diff line change
Expand Up @@ -15,7 +15,7 @@
for f in input_files:
for s in list(SeqIO.parse(f, "fasta")):
sample = s.id.split(".")[0]
aa_feat = s.id.split(".")[1]
aa_feat = ".".join(s.id.split(".")[1:])
# select the rows that contain the sample name from snakemake.params.space
data = space.loc[space["sample"] == sample]
data = data.explode("AA_FEAT_NAMES").reset_index(drop=True)[
Expand Down

0 comments on commit 3e14dd5

Please sign in to comment.