Commit da43f726 authored by Mikaël Salson's avatar Mikaël Salson

germline: display the real information on germlines

The information on germlines were recovered through the index which does
not make sense with an Aho-Corasick automaton. In such a case the index
is the same but the germlines have different seeds and therefore index
loads should differ.

Since getIndexLoad depends on the KmerAffect, it also differs on V and J
parent 0de26304
...@@ -214,10 +214,14 @@ ostream &operator<<(ostream &out, const Germline &germline) ...@@ -214,10 +214,14 @@ ostream &operator<<(ostream &out, const Germline &germline)
<< setw(3) << germline.delta_min << setw(3) << germline.delta_min
<< " "; << " ";
size_t seed_w = seed_weight(germline.seed);
if (germline.index) { if (germline.index) {
out << " 0x" << hex << setw(2) << setfill('0') << germline.index->id << dec << setfill(' ') << " " ; out << " 0x" << hex << setw(2) << setfill('0') << germline.index->id << dec << setfill(' ') << " " ;
out << fixed << setprecision(3) << setw(8) << 100 * germline.index->getIndexLoad() << "%"; out << fixed << setprecision(3) << setw(8)
out << " l" << germline.index->getS() << " k" << germline.index->getK() << " " << germline.index->getSeed() ; // TODO: there should be a << for index << 100 * germline.index->getIndexLoad(KmerAffect(germline.affect_5, 1, seed_w)) << "%" << " "
<< 100 * germline.index->getIndexLoad(KmerAffect(germline.affect_3, 1, seed_w)) << "%";
out << " l" << germline.seed.length() << " k" << seed_w << " " << germline.seed ;
} }
out << endl; out << endl;
......
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