[BioC] Consensus sequence

James W. MacDonald jmacdon at uw.edu
Tue Sep 9 04:32:25 CEST 2014

I have a use case that I am not sure the best way to proceed.

I have a non-model organism for which we would like to know the sequence
for a set of genes. It's a bird species, and I can for example align to the
zebra finch genome, then run bam_tally over a region to get  variants where
my species differs from zebra finch, and then infer the sequence based on
the reference and the variants.

But that seems harder than it should be. Is there some function that I am
missing that I can feed a GRanges and get back the consensus sequence for
that region, based on say a simple vote of the reads that overlap it?



	[[alternative HTML version deleted]]

More information about the Bioconductor mailing list