X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=sffinfocommand.cpp;h=cb5f87410e64583afa710ba034448e15cffb8fd1;hb=a3d9877d9c62058dbe167ceca159f4b76744f1cf;hp=8c0e473a87e86453c8721b0f5c5b5bde947b703b;hpb=b4f80c1d2be78a8743a408a2b6d462b07f9f71ff;p=mothur.git diff --git a/sffinfocommand.cpp b/sffinfocommand.cpp index 8c0e473..cb5f874 100644 --- a/sffinfocommand.cpp +++ b/sffinfocommand.cpp @@ -584,16 +584,23 @@ int SffInfoCommand::printFastaSeqData(ofstream& out, seqRead& read, Header& head string seq = read.bases; - if (trim) { - seq = seq.substr((header.clipQualLeft-1), (header.clipQualRight-header.clipQualLeft+1)); + if(header.clipQualRight != 0){ + seq = seq.substr((header.clipQualLeft-1), (header.clipQualRight-header.clipQualLeft+1)); + } + else{ + seq = seq.substr(header.clipQualLeft-1); + } }else{ //if you wanted the sfftxt then you already converted the bases to the right case if (!sfftxt) { //make the bases you want to clip lowercase and the bases you want to keep upper case for (int i = 0; i < (header.clipQualLeft-1); i++) { seq[i] = tolower(seq[i]); } for (int i = (header.clipQualLeft-1); i < (header.clipQualRight-1); i++) { seq[i] = toupper(seq[i]); } - for (int i = (header.clipQualRight-1); i < seq.length(); i++) { seq[i] = tolower(seq[i]); } + + if(header.clipQualRight != 0){ + for (int i = (header.clipQualRight-1); i < seq.length(); i++) { seq[i] = tolower(seq[i]); } + } } }