cout << "DONE." << endl << endl; cout.flush();
}
+/**************************************************************************************************/
+
+KmerDB::~KmerDB(){
+
+ //for (int i = 0; i < templateSequences.size(); i++) { delete templateSequences[i]; }
+ // templateSequences.clear();
+}
/**************************************************************************************************/
-Sequence* KmerDB::findClosestSequence(Sequence* candidateSeq){
+Sequence KmerDB::findClosestSequence(Sequence* candidateSeq){
Kmer kmer(kmerSize);
for(int i=0;i<numSeqs;i++){ // for all of the template sequences...
- string seq = templateSequences[i]->getUnaligned(); // ...take the unaligned sequence...
+ string seq = templateSequences[i].getUnaligned(); // ...take the unaligned sequence...
int numKmers = seq.length() - kmerSize + 1;
vector<int> seenBefore(maxKmer+1,0);