[Seqan-dev] obtain consensus sequence
: Erik Aronesty <
: Wed, 12 Oct 2011 16:03:08 -0400
: SeqAn Development <
: [Seqan-dev] obtain consensus sequence
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 I'd like to get the consensus sequence out of "sub" such that at each position, the majority rules sequence.
Re: [Seqan-dev] obtain consensus sequence
Manuel Holtgrewe <email@example.com>
seqan-dev - October 2011 - Archives indexes sorted by:
[ thread ]
[ subject ]
[ author ]
[ date ]
Complete archive of the seqan-dev mailing list
More info on this list...