Difference between revisions of "830 - BestAlignment::EvaluateAlignment() (main loop)"
From Genome Analysis Wiki
Jump to navigationJump to search (Created page with "<source lang="cpp"> void BestAlignment::EvaluateAlignment(MaqIndex & index, FastQ & read, unsigned int offset, bool fwd) { unsigned int start, stop; index.LookupWord(read....") |
(No difference)
|
Latest revision as of 13:46, 28 October 2013
void BestAlignment::EvaluateAlignment(MaqIndex & index, FastQ & read, unsigned int offset, bool fwd)
{
unsigned int start, stop;
index.LookupWord(read.sequence, offset, start, stop);
for (unsigned int i = start; i <= stop; i++)
EvaluateAlignment(index.reference, read, index.sortedPositions[i] - offset, fwd);
}