1,555
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search830 - BestAlignment::EvaluateAlignment() (main loop) (view source)
Revision as of 13:46, 28 October 2013
, 13:46, 28 October 2013Created page with "<source lang="cpp"> void BestAlignment::EvaluateAlignment(MaqIndex & index, FastQ & read, unsigned int offset, bool fwd) { unsigned int start, stop; index.LookupWord(read...."
<source lang="cpp">
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);
}
</source>
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);
}
</source>
