vector<string> queryAlign(numRefSeqs);
vector<string> refAlign(numRefSeqs);
- vector<vector<int> > leftDiffs(numRefSeqs, 0);
- vector<vector<int> > rightDiffs(numRefSeqs, 0);
- vector<vector<int> > leftMaps(numRefSeqs, 0);
- vector<vector<int> > rightMaps(numRefSeqs, 0);
+ vector<vector<int> > leftDiffs(numRefSeqs);
+ vector<vector<int> > rightDiffs(numRefSeqs);
+ vector<vector<int> > leftMaps(numRefSeqs);
+ vector<vector<int> > rightMaps(numRefSeqs);
int bestRefIndex = -1;
int bestRefDiffs = numeric_limits<int>::max();
vector<int> bestRight(seqLength, -1);
for(int l=0;l<seqLength;l++){
-
+
for(int i=0;i<numRefSeqs;i++){
if(rightDiffs[i][l] < singleRight[l]){
singleRight[l] = rightDiffs[i][l];