Changeset 3409 for trunk/test/bam_header.cc

Timestamp:

Apr 17, 2015, 3:07:04 AM (8 years ago)

Author:

Peter

Message:

refs #817. Fix bug that struct was not updated in text(1) function (in libbam mode); add test to catch bug.

File:

• trunk/test/bam_header.cc (modified) (2 diffs)

trunk/test/bam_header.cc

@@ -1 +1 @@
 // $Id$
 //
-// Copyright (C) 2013, 2014 Peter Johansson
+// Copyright (C) 2013, 2014, 2015 Peter Johansson
 //
 // This program is free software; you can redistribute it and/or modify
@@ -100 +100 @@
   }
 
+  int32_t n_chr = hdr.n_targets();
+
+  // check that we can add another chromosomevia a @SG field
+  std::stringstream ss(hdr.text());
+  std::string str3;
+  std::string line;
+  int32_t n_targ = 0;
+  while (getline(ss, line)) {
+    str3 += line + "\n";
+    if (line.substr(0,3) == "@SQ") {
+      ++n_targ;
+      if (n_targ == n_chr)
+        str3 += "@SQ\tSN:Zbababababa\tLN:123\n";
+    }
+  }
+  hdr.text(str3);
+
+  int32_t n = hdr.n_targets();
+  if (n != n_chr+1) {
+    suite.add(false);
+    suite.err() << "error: incorrect number of targets: " << n
+                << " expected " << (n_chr+1) << "\n";
+    suite.err() << "header text:\n" << hdr.text() << "===\n";
+  }
+
 #endif
 }