Merge pull request #156 from sarahet/fix_bam_header

h-2 · web-flow · commit 52fd0bb00f85 · 2020-01-30T14:22:57.000Z
[Fix] Get correct number of header entries for subject sequences
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`