X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=trimseqscommand.cpp;h=bb4506385774260141692e0ee5991ac5e10a9ee1;hb=905cc2b0bd18c5ce611b048d785e93859865a5ea;hp=ef7dfdd0f45d649cc6296d91e106e75b9a6fd167;hpb=dabb32850eae34e26b066a7da55716ad685786a3;p=mothur.git diff --git a/trimseqscommand.cpp b/trimseqscommand.cpp index ef7dfdd..bb45063 100644 --- a/trimseqscommand.cpp +++ b/trimseqscommand.cpp @@ -778,10 +778,14 @@ bool TrimSeqsCommand::compareDNASeq(string oligo, string seq){ bool TrimSeqsCommand::stripQualThreshold(Sequence& seq, ifstream& qFile){ try { string rawSequence = seq.getUnaligned(); - int seqLength = rawSequence.length(); - string name; + int seqLength; // = rawSequence.length(); + string name, temp, temp2; - qFile >> name; + qFile >> name >> temp; + + splitAtEquals(temp2, temp); //separates length=242, temp=length, temp2=242 + convert(temp, seqLength); //converts string to int + if (name.length() != 0) { if(name.substr(1) != seq.getName()) { m->mothurOut("sequence name mismatch btwn fasta and qual file"); m->mothurOutEndLine(); } } while (!qFile.eof()) { char c = qFile.get(); if (c == 10 || c == 13){ break; } }