Skip to content

Commit

Permalink
Update k-score sensitivity fitting for no-cntxt profile searches
Browse files Browse the repository at this point in the history
  • Loading branch information
milot-mirdita committed Sep 27, 2021
1 parent 3e92abf commit f3f5b13
Showing 1 changed file with 6 additions and 6 deletions.
12 changes: 6 additions & 6 deletions src/prefiltering/Prefiltering.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -1001,14 +1001,14 @@ int Prefiltering::getKmerThreshold(const float sensitivity, const bool isProfile
}
} else if(isProfile && hasContextPseudoCnts == false) {
if (kmerSize == 5) {
float base = 117.35;
kmerThrBest = base - (sensitivity * 5.85);
float base = 108.8;
kmerThrBest = base - (sensitivity * 4.7);
} else if (kmerSize == 6) {
float base = 135.75;
kmerThrBest = base - (sensitivity * 6.75);
float base = 134.35;
kmerThrBest = base - (sensitivity * 6.15);
} else if (kmerSize == 7) {
float base = 152.5;
kmerThrBest = base - (sensitivity * 8);
float base = 149.15;
kmerThrBest = base - (sensitivity * 6.85);
} else {
Debug(Debug::ERROR) << "The k-mer size " << kmerSize << " is not valid.\n";
EXIT(EXIT_FAILURE);
Expand Down

0 comments on commit f3f5b13

Please sign in to comment.