X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=seqerrorcommand.cpp;fp=seqerrorcommand.cpp;h=ecd6f0ceb723290337f5856f724710eff7a751ee;hb=5ca84bbe4393adecc30dc05d64e79091ef35c40f;hp=641a3abce04bde1b8002e11af636fcd5e58bc6e2;hpb=2a29fceeeb8754c3fd97ba830d2fbed5c4349ee8;p=mothur.git diff --git a/seqerrorcommand.cpp b/seqerrorcommand.cpp index 641a3ab..ecd6f0c 100644 --- a/seqerrorcommand.cpp +++ b/seqerrorcommand.cpp @@ -680,10 +680,10 @@ int SeqErrorCommand::driver(string filename, string qFileName, string rFileName, int numSeqs = 0; map::iterator it; - qScoreErrorMap['m'].assign(41, 0); - qScoreErrorMap['s'].assign(41, 0); - qScoreErrorMap['i'].assign(41, 0); - qScoreErrorMap['a'].assign(41, 0); + qScoreErrorMap['m'].assign(101, 0); + qScoreErrorMap['s'].assign(101, 0); + qScoreErrorMap['i'].assign(101, 0); + qScoreErrorMap['a'].assign(101, 0); errorForward['m'].assign(maxLength,0); errorForward['s'].assign(maxLength,0); @@ -719,8 +719,8 @@ int SeqErrorCommand::driver(string filename, string qFileName, string rFileName, qualForwardMap.resize(maxLength); qualReverseMap.resize(maxLength); for(int i=0;i > qScoreErrorMap){ outputNames.push_back(errorQualityFileName); outputTypes["errorquality"].push_back(errorQualityFileName); errorQualityFile << "qscore\tmatches\tsubstitutions\tinsertions\tambiguous" << endl; - for(int i=0;i<41;i++){ + for(int i=0;i<101;i++){ errorQualityFile << i << '\t' << qScoreErrorMap['m'][i] << '\t' << qScoreErrorMap['s'][i] << '\t' << qScoreErrorMap['i'][i] << '\t'<< qScoreErrorMap['a'][i] << endl; } errorQualityFile.close();