Bug fix in BamReader::Jump()

diff --git a/BamReader.cpp b/BamReader.cpp
index be6b72b51b5da6ae8325d85eec6335635e1b4993..b1b1b5de2ab79a76416bd5267f6dc5cc1ecccd02 100644 (file)
--- a/BamReader.cpp
+++ b/BamReader.cpp
@@ -540,8 +540,10 @@ bool BamReader::BamReaderPrivate::Jump(int refID, int position) {
         \r
         // if this alignment corresponds to desired position\r
         // return success of seeking back to 'current offset'\r
-        if ( (bAlignment.RefID == refID && bAlignment.Position + bAlignment.Length > position) || (bAlignment.RefID > refID) )\r
+        if ( (bAlignment.RefID == refID && bAlignment.Position + bAlignment.Length > position) || (bAlignment.RefID > refID) ) {\r
+            if ( o != offsets.begin() ) --o;\r
             return mBGZF.Seek(*o);\r
+        }\r
     }\r
     \r
     return result;\r