Command:
cons -plurality 0.0 -setcase 9.9 cons-bug.fasta
Output:
ATGAGggGGGGggggg-aGTGAGCGACGTCGCGTGAAGCGCGGCTATCTTCGCTTCGCCA
Both "G" and "A" get score 3 at column 17, but cons produced "-" for this column.
"-" is not described as possible output in the cons documentation. We propose to print one of input letters ("G" or "A") instead.
EMBOSS 6.3.1 (Debian)
input file for cons
hanks for spoting this one.
This situation is addressed by prettyplot using the collisions option which allows one of the alternatives to be used in the consensus. We used the same functions for cons but some additional parameters will be needed to ensure a consensus is filled in.
Certainly something has to be chosen because the cons output should be a valid sequence and '-' is treated as a gap character.