X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=chimera.cpp;fp=chimera.cpp;h=85b0b38c8daab23e44c5c676ed1f38009afb7a62;hp=9ad287b68d8457e6dffd6e101bbca007c64745eb;hb=9c23307c583d4e8595f75278c13e708788f2f058;hpb=81276c241b984898f8d30ad123c00592ee6db7b8 diff --git a/chimera.cpp b/chimera.cpp index 9ad287b..85b0b38 100644 --- a/chimera.cpp +++ b/chimera.cpp @@ -22,14 +22,14 @@ string Chimera::createFilter(vector seqs, float t) { vector t; t.resize(seqs[0]->getAligned().length(), 0); vector g; g.resize(seqs[0]->getAligned().length(), 0); vector c; c.resize(seqs[0]->getAligned().length(), 0); - + filterString = (string(seqs[0]->getAligned().length(), '1')); //for each sequence for (int i = 0; i < seqs.size(); i++) { string seqAligned = seqs[i]->getAligned(); - + for (int j = 0; j < seqAligned.length(); j++) { //if this spot is a gap if ((seqAligned[j] == '-') || (seqAligned[j] == '.')) { gaps[j]++; } @@ -47,10 +47,8 @@ string Chimera::createFilter(vector seqs, float t) { if(gaps[i] == seqs.size()) { filterString[i] = '0'; numColRemoved++; } else if (((a[i] < threshold) && (t[i] < threshold) && (g[i] < threshold) && (c[i] < threshold))) { filterString[i] = '0'; numColRemoved++; } - //cout << "a = " << a[i] << " t = " << t[i] << " g = " << g[i] << " c = " << c[i] << endl; + cout << "a = " << a[i] << " t = " << t[i] << " g = " << g[i] << " c = " << c[i] << endl; } - -//cout << "filter = " << filterString << endl; mothurOut("Filter removed " + toString(numColRemoved) + " columns."); mothurOutEndLine(); return filterString; @@ -95,7 +93,7 @@ vector Chimera::readSeqs(string file) { openInputFile(file, in); vector container; int count = 0; - int length = 0; + length = 0; unaligned = false; //read in seqs and store in vector