1,555
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search830 - MaqIndex::GetWord() (view source)
Revision as of 13:39, 28 October 2013
, 13:39, 28 October 2013Created page with "<source lang="cpp"> long long MaqIndex::GetWord(String & sequence, unsigned int position) { long long word = 0; if (position + mask.Length() > sequence.Length()) return -..."
<source lang="cpp">
long long MaqIndex::GetWord(String & sequence, unsigned int position)
{
long long word = 0;
if (position + mask.Length() > sequence.Length())
return -1;
for (unsigned int i = 0; i < mask.Length(); i++)
if (mask[i] == '1')
{
char base = sequence[position + i];
if (base < 1 || base > 4) return -1;
word = word * 4 + base - 1;
}
return word;
}
long long MaqIndex::GetReferenceWord(unsigned int position)
{
return GetWord(reference, position);
}
</source>
long long MaqIndex::GetWord(String & sequence, unsigned int position)
{
long long word = 0;
if (position + mask.Length() > sequence.Length())
return -1;
for (unsigned int i = 0; i < mask.Length(); i++)
if (mask[i] == '1')
{
char base = sequence[position + i];
if (base < 1 || base > 4) return -1;
word = word * 4 + base - 1;
}
return word;
}
long long MaqIndex::GetReferenceWord(unsigned int position)
{
return GetWord(reference, position);
}
</source>
