Commit 37fc6576 authored by Mathieu Giraud's avatar Mathieu Giraud

core/germline.{h,cpp}: Germline.mark_as_ambiguous(other)

Marks as AFFECT_AMBIGUOUS kmers from another germline
parent 8746e51b
......@@ -98,6 +98,20 @@ void Germline::update_index()
cout << " --- index " << index << " updated " << affect_5 << "/" << affect_4 << "/" << affect_3 << endl;
}
void Germline::mark_as_ambiguous(Germline *other)
{
index->insert(other->rep_5, AFFECT_AMBIGUOUS_SYMBOL);
if (other->affect_4.size())
index->insert(other->rep_4, AFFECT_AMBIGUOUS_SYMBOL);
index->insert(other->rep_3, AFFECT_AMBIGUOUS_SYMBOL);
cout << " --- index " << index << " - " << code << " marked " << other->code << endl;
}
Germline::~Germline()
{
if (index)
......
......@@ -47,6 +47,8 @@ class Germline {
void new_index(string seed);
void use_index(IKmerStore<KmerAffect> *index);
void mark_as_ambiguous(Germline *other);
list <string> f_reps_5 ;
list <string> f_reps_4 ;
list <string> f_reps_3 ;
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment