[Seqan-dev] obtain consensus sequence

[Erik Aronesty <erik@q32.com>]

i have a list of sequences that i ran a globalmsa on

i can output the alignment and see that it worked ok

i would like obtain the "consensus sequence" for it... 

i was going to walk through each alignment and record the "majority rules" for each position 

is this the best way to do it... it seems there a tool that's built (consensus.h)... but not sure how to use it:

I do this:

.... appendValue(rows(sub),fq.seq+p) ...

Then this:

 globalMsaAlignment(sub);

Then I'd like to get the consensus sequence out of "sub" such that at each position, the majority rules sequence.