Sha256: 2b275ac35d005733b072326698dbf8026371da01d352f310c10b410b4782b723
Contents?: true
Size: 879 Bytes
Versions: 183
Compression:
Stored size: 879 Bytes
Contents
module hamming;
import std.string;
import std.algorithm.comparison : mismatch;
import std.exception;
int distance (string lhs, string rhs)
{
if (lhs.length != rhs.length)
{
throw new Exception(format("DNA sequence with different lengths: %s - %s", lhs.length, rhs.length));
}
int count = 0;
auto m = mismatch(lhs, rhs);
while (!m[0].empty)
{
++count;
m = mismatch(m[0][1 .. $], m[1][1 .. $]);
}
return count;
}
unittest
{
assert(distance("A", "A") == 0);
assert(distance("A", "G") == 1);
assert(distance("AG", "CT") == 2);
assert(distance("AT", "CT") == 1);
assert(distance("GGACG", "GGTCG") == 1);
assertThrown(distance("AAAG", "AAA"));
assertThrown(distance("AAA", "AAAG"));
assert(distance("GATACA", "GCATAA") == 4);
assert(distance("GGACGGATTCTG", "AGGACGGATTCT") == 9);
}
void main ()
{
}
Version data entries
183 entries across 183 versions & 1 rubygems