Commit 43b105e9 authored by Mathieu Giraud's avatar Mathieu Giraud

core/germline.h, core/segment.cpp: remove MAX12 e-value penalty

Discards 7b1032e4.
As noted by @mikael-s, the MAX12 index has a higher index load
than the regular germlines. There is no need to further penalize it.
parent 351d768f
......@@ -10,7 +10,6 @@
#include "tools.h"
#define PSEUDO_GERMLINE_MAX12 "xxx"
#define PSEUDO_GERMLINE_MAX12_EVALUE_PENALTY 1000
using namespace std;
......
......@@ -270,9 +270,6 @@ KmerSegmenter::KmerSegmenter(Sequence seq, Germline *germline, double threshold,
// Anyway, it allows to gather such reads and their reverse complement into a unique window...
// ... except when the read is quite different outside the window
strand = nb_strand[0] > nb_strand[1] ? -1 : 1 ;
// The pseudo-germline should never take precedence over regular germlines of similar e-value
multiplier *= PSEUDO_GERMLINE_MAX12_EVALUE_PENALTY ;
}
else
......
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