Skip to content

Commit

Permalink
Default to --infer-ambiguous
Browse files Browse the repository at this point in the history
This keeps default augur behavior unchanged despite the TreeTime upgrade. By adding the --keep-ambiguous flag, new behavior can be requested. This seemed safer in terms of not causing unexpected behavior for users.
  • Loading branch information
trvrb committed Feb 12, 2020
1 parent 909d847 commit aab713c
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions augur/ancestral.py
Original file line number Diff line number Diff line change
Expand Up @@ -106,9 +106,9 @@ def register_arguments(parser):
parser.add_argument('--vcf-reference', type=str, help='fasta file of the sequence the VCF was mapped to')
parser.add_argument('--output-vcf', type=str, help='name of output VCF file which will include ancestral seqs')
ambiguous = parser.add_mutually_exclusive_group()
ambiguous.add_argument('--keep-ambiguous', action="store_false", dest='infer-ambiguous',
ambiguous.add_argument('--keep-ambiguous', action="store_false", default=False, dest='infer-ambiguous',
help='do not infer nucleotides at ambiguous (N) sites on tip sequences (leave as N).')
ambiguous.add_argument('--infer-ambiguous', action="store_true", default=False,
ambiguous.add_argument('--infer-ambiguous', action="store_true",
help='infer nucleotides at ambiguous (N,W,R,..) sites on tip sequences and replace with most likely state.')
parser.add_argument('--keep-overhangs', action="store_true", default=False,
help='do not infer nucleotides for gaps (-) on either side of the alignment')
Expand Down

0 comments on commit aab713c

Please sign in to comment.