vector<int> topMatches;
ofstream queryFile;
- openOutputFile((queryFileName+seq->getName()), queryFile);
+ m->openOutputFile((queryFileName+seq->getName()), queryFile);
queryFile << '>' << seq->getName() << endl;
queryFile << seq->getUnaligned() << endl;
queryFile.close();
system(blastCommand.c_str());
ifstream m8FileHandle;
- openInputFile(blastFileName+seq->getName(), m8FileHandle, "no error");
+ m->openInputFile(blastFileName+seq->getName(), m8FileHandle, "no error");
string dummy;
int templateAccession;
- gobble(m8FileHandle);
+ m->gobble(m8FileHandle);
while(!m8FileHandle.eof()){
m8FileHandle >> dummy >> templateAccession >> searchScore;
//get rest of junk in line
while (!m8FileHandle.eof()) { char c = m8FileHandle.get(); if (c == 10 || c == 13){ break; } }
- gobble(m8FileHandle);
+ m->gobble(m8FileHandle);
topMatches.push_back(templateAccession);
}
m8FileHandle.close();
vector<int> topMatches;
ofstream queryFile;
- openOutputFile((queryFileName+seq->getName()), queryFile);
+ m->openOutputFile((queryFileName+seq->getName()), queryFile);
queryFile << '>' << seq->getName() << endl;
queryFile << seq->getUnaligned() << endl;
queryFile.close();
system(blastCommand.c_str());
ifstream m8FileHandle;
- openInputFile(blastFileName+seq->getName(), m8FileHandle, "no error");
+ m->openInputFile(blastFileName+seq->getName(), m8FileHandle, "no error");
string dummy;
int templateAccession;
- gobble(m8FileHandle);
+ m->gobble(m8FileHandle);
while(!m8FileHandle.eof()){
m8FileHandle >> dummy >> templateAccession >> searchScore;
//get rest of junk in line
while (!m8FileHandle.eof()) { char c = m8FileHandle.get(); if (c == 10 || c == 13){ break; } }
- gobble(m8FileHandle);
+ m->gobble(m8FileHandle);
topMatches.push_back(templateAccession);
//cout << templateAccession << endl;
}
try {
ofstream unalignedFastaFile;
- openOutputFileAppend(dbFileName, unalignedFastaFile);
+ m->openOutputFileAppend(dbFileName, unalignedFastaFile);
// generating a fasta file with unaligned template
unalignedFastaFile << '>' << count << endl; // sequences, which will be input to formatdb