X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=chimerarealigner.cpp;fp=chimerarealigner.cpp;h=bfd283db0d1b66db7df4606d980cc9d54dad58a7;hb=9489965363593bb2a3e94f801b4079a32ddf8732;hp=48d9bc16a14b3afa8ad7588141cc9ffaaf2e03c2;hpb=62aba27bfa0c2036493b012deff73246d9ac130f;p=mothur.git diff --git a/chimerarealigner.cpp b/chimerarealigner.cpp index 48d9bc1..bfd283d 100644 --- a/chimerarealigner.cpp +++ b/chimerarealigner.cpp @@ -49,15 +49,15 @@ void ChimeraReAligner::reAlign(Sequence* query, vector parents) { parentParts.push_back(parent); - if (q.length() > longest) { longest = q.length(); } - if (p.length() > longest) { longest = p.length(); } + if (queryFrag->getUnaligned().length() > longest) { longest = queryFrag->getUnaligned().length(); } + if (parent->getUnaligned().length() > longest) { longest = parent->getUnaligned().length(); } } //align each peice to correct parent from results for (int i = 0; i < queryParts.size(); i++) { if ((queryParts[i]->getUnaligned() == "") || (parentParts[i]->getUnaligned() == "")) {;} else { - alignment = new NeedlemanOverlap(-2.0, match, misMatch, longest+1); //default gapopen, match, mismatch, longestbase + Alignment* alignment = new NeedlemanOverlap(-2.0, 1.0, -1.0, longest+1); //default gapopen, match, mismatch, longestbase Nast nast(alignment, queryParts[i], parentParts[i]); delete alignment;