[Fix] Get correct number of header entries for subject sequences

sarahet · sarahet · commit d48cc93e1b3a · 2020-01-30T05:55:44.000+01:00
diff --git a/src/search_output.hpp b/src/search_output.hpp
@@ -330,7 +330,7 @@ myWriteHeader(TGH & globalHolder, TLambdaOptions const & options)
         if (sIsTranslated(TGH::blastProgram))
         {
             //TODO can we get around a copy?
-            subjSeqLengths = globalHolder.untransSubjSeqLengths;
+            subjSeqLengths = prefix(globalHolder.untransSubjSeqLengths, length(globalHolder.untransSubjSeqLengths) - 1);
         } else
         {
             // compute lengths ultra-fast

Original file line number	Diff line number	Diff line change
`@@ -330,7 +330,7 @@ myWriteHeader(TGH & globalHolder, TLambdaOptions const & options)`
`330`	`330`	`if (sIsTranslated(TGH::blastProgram))`
`331`	`331`	`{`
`332`	`332`	`//TODO can we get around a copy?`
`333`		`- subjSeqLengths = globalHolder.untransSubjSeqLengths;`
	`333`	`+ subjSeqLengths = prefix(globalHolder.untransSubjSeqLengths, length(globalHolder.untransSubjSeqLengths) - 1);`
`334`	`334`	`} else`
`335`	`335`	`{`
`336`	`336`	`// compute lengths ultra-fast`