Changeset 3055 for trunk/yat/omic/BamRead.cc

Timestamp:

Jun 16, 2013, 2:21:30 AM (10 years ago)

Author:

Peter

Message:

closes #746. New function BamRead::name(1)

File:

• trunk/yat/omic/BamRead.cc (modified) (1 diff)

trunk/yat/omic/BamRead.cc

@@ -181 +181 @@
   const char* BamRead::name(void) const
   { return bam1_qname(bam_); }
+
+
+  void BamRead::name(const std::string& n)
+  {
+    int offset = n.size() + 1 - core().l_qname;
+    assert(bam_);
+    reserve(bam_->data_len + offset);
+    // move remaining data
+    if (offset)
+      memmove(bam_->data + core().l_qname + offset,
+              bam_->data + core().l_qname,
+              bam_->data_len - core().l_qname);
+    core().l_qname += offset;
+    // copy new name
+    memcpy(bam1_qname(bam_), n.c_str(), n.size()+1);
+    bam_->data_len += offset;
+    assert(bam_->data_len <= bam_->m_data);
+  }