//bootstrap - to set confidenceScore
int numToSelect = queryKmers.size() / 8;
-// tax = bootstrapResults(queryKmers, index, numToSelect);
+ tax = bootstrapResults(queryKmers, index, numToSelect);
return tax;
}
//find taxonomy with highest probability that this sequence is from it
- cout << genusNodes.size() << endl;
+// cout << genusNodes.size() << endl;
for (int k = 0; k < genusNodes.size(); k++) {
prob += wordGenusProb[queryKmer[i]][k];
}
- cout << phyloTree->get(genusNodes[k]).name << '\t' << prob << endl;
+// cout << phyloTree->get(genusNodes[k]).name << '\t' << prob << endl;
//is this the taxonomy with the greatest probability?
if (prob > maxProbability) {