From Genome Analysis Wiki
Jump to navigationJump to search
345 bytes added
, 13:34, 28 October 2013
<source lang="cpp">
int MaqIndex::CompareToIndex(String & sequence, unsigned int offset, unsigned int index)
{
index = sortedPositions[index];
for (unsigned int i = 0; i < mask.Length(); i++)
if (mask[i] == '1' && sequence[offset + i] != reference[index + i])
return sequence[offset + i] - reference[index + i];
return 0;
}
</source>