[FIX] sam/bam with new extension modes

h-2 · h-2 · commit 770c963afc6c · 2016-08-19T12:51:59.000+02:00
diff --git a/src/lambda.hpp b/src/lambda.hpp
@@ -2038,6 +2038,9 @@ iterateMatchesFullSimd(TLocalHolder & lH)
         appendValue(depSetV, source(bm.alignRow1));
 
         _setFrames(bm, *it, lH);
+
+        bm._n_qId = it->qryId / qNumFrames(TGlobalHolder::blastProgram);
+        bm._n_sId = it->subjId / sNumFrames(TGlobalHolder::blastProgram);
     }
 
     // fill up last batch
@@ -2219,6 +2222,9 @@ iterateMatchesFullSerial(TLocalHolder & lH)
         }
 
         _setFrames(bm, m, lH);
+
+        bm._n_qId = it->qryId / qNumFrames(TGlobalHolder::blastProgram);
+        bm._n_sId = it->subjId / sNumFrames(TGlobalHolder::blastProgram);
     }
 
     _writeRecord(record, lH);

Original file line number	Diff line number	Diff line change
`@@ -2038,6 +2038,9 @@ iterateMatchesFullSimd(TLocalHolder & lH)`
`2038`	`2038`	`appendValue(depSetV, source(bm.alignRow1));`
`2039`	`2039`
`2040`	`2040`	`_setFrames(bm, *it, lH);`
	`2041`	`+`
	`2042`	`+ bm._n_qId = it->qryId / qNumFrames(TGlobalHolder::blastProgram);`
	`2043`	`+ bm._n_sId = it->subjId / sNumFrames(TGlobalHolder::blastProgram);`
`2041`	`2044`	`}`
`2042`	`2045`
`2043`	`2046`	`// fill up last batch`
`@@ -2219,6 +2222,9 @@ iterateMatchesFullSerial(TLocalHolder & lH)`
`2219`	`2222`	`}`
`2220`	`2223`
`2221`	`2224`	`_setFrames(bm, m, lH);`
	`2225`	`+`
	`2226`	`+ bm._n_qId = it->qryId / qNumFrames(TGlobalHolder::blastProgram);`
	`2227`	`+ bm._n_sId = it->subjId / sNumFrames(TGlobalHolder::blastProgram);`
`2222`	`2228`	`}`
`2223`	`2229`
`2224`	`2230`	`_writeRecord(record, lH);`