Help text does not explain that a reference sequence is needed to evaluate substitutions.
Are we arbitrarily using first sequence in input FASTA?
Ideally the script should generate a consensus sequence as the reference.
Yes, currently hypermut.py and the LANL's hypermut designate the first sequence as the reference sequence. I will update the help text and add an option for generating the consensus sequence.
Help text does not explain that a reference sequence is needed to evaluate substitutions. Are we arbitrarily using first sequence in input FASTA? Ideally the script should generate a consensus sequence as the reference.