Add test for ancestral --genes file

nextstrain · Dec 20, 2023 · aae0640 · aae0640
1 parent b6f63ff
commit aae0640
Show file tree

Hide file tree

Showing 2 changed files with 24 additions and 0 deletions.
diff --git a/tests/functional/ancestral/cram/infer-amino-acid-sequences-genes-file.t b/tests/functional/ancestral/cram/infer-amino-acid-sequences-genes-file.t
@@ -0,0 +1,22 @@
+Setup
+
+  $ source "$TESTDIR"/_setup.sh
+
+Infer ancestral nucleotide and amino acid sequences, using a genes file.
+
+  $ ${AUGUR} ancestral \
+  >  --tree $TESTDIR/../data/tree.nwk \
+  >  --alignment $TESTDIR/../data/aligned.fasta \
+  >  --annotation $TESTDIR/../data/zika_outgroup.gb \
+  >  --genes $TESTDIR/../data/genes.txt \
+  >  --translations $TESTDIR/../data/aa_sequences_%GENE.fasta \
+  >  --output-node-data ancestral_mutations.json \
+  >  --output-sequences ancestral_sequences.fasta \
+  >  --output-translations ancestral_aa_sequences_%GENE.fasta > /dev/null
+
+Check that the reference length was correctly exported as the nuc annotation
+
+  $ grep -E "\"(ENV|PRO|nuc)\": {" ancestral_mutations.json
+      "ENV": {
+      "PRO": {
+      "nuc": {
diff --git a/tests/functional/ancestral/data/genes.txt b/tests/functional/ancestral/data/genes.txt
@@ -0,0 +1,2 @@
+ENV
+PRO