From: pschloss <pschloss>
Date: Wed, 23 Mar 2011 12:07:55 +0000 (+0000)
Subject: small mod to chimerarefseq.cpp
X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=commitdiff_plain;h=a11ad41867e19c38377e04b1afc1ca82900c4f82

small mod to chimerarefseq.cpp
---

diff --git a/refchimeratest.cpp b/refchimeratest.cpp
index 4d4a658..aab17e6 100644
--- a/refchimeratest.cpp
+++ b/refchimeratest.cpp
@@ -60,12 +60,8 @@ int RefChimeraTest::analyzeQuery(string queryName, string querySeq){
 	int nMera = 0;
 	string chimeraRefSeq = "";
 	
-	if(bestSequenceMismatch - minMismatchToChimera <= 3){
-		nMera = 1;
-		chimeraRefSeq = referenceSeqs[bestMatch];
-	}
-	else {
-
+	if(bestSequenceMismatch - minMismatchToChimera > 3 || (minMismatchToChimera == 0 && bestSequenceMismatch != 0)){
+	
 		nMera = 2;
 		chimeraRefSeq = stitchBimera(leftParentBi, rightParentBi, breakPointBi);
 		
@@ -81,6 +77,12 @@ int RefChimeraTest::analyzeQuery(string queryName, string querySeq){
 //		}
 		
 	}
+	else{
+		nMera = 1;
+		chimeraRefSeq = referenceSeqs[bestMatch];
+	}
+	
+	
 	double distToChimera = calcDistToChimera(querySeq, chimeraRefSeq);
 	
 //	double loonIndex = calcLoonIndex(querySeq, referenceSeqs[leftParentBi], referenceSeqs[rightParentBi], breakPointBi, binMatrix);