Changeset 2988

Timestamp:

Feb 23, 2013, 6:24:31 AM (11 years ago)

Author:

Peter

Message:

bam1_seq_seti was not defined in samtools 0.1.18. This workarounds fixes that.

File:

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

trunk/yat/omic/BamRead.cc

@@ -26 +26 @@
 #include <cassert>
 #include <sstream>
+
+// bam1_seq_seti is not defined in 0.1.18 (or older), so
+// backport. Code is taken bam.h in version 0.1.19-dev.
+#ifndef bam1_seq_seti
+#define bam1_seq_seti(s, i, c) ( (s)[(i)>>1] = \
+  ((s)[(i)>>1] & 0xf<<(((i)&1)<<2)) | (c)<<((~(i)&1)<<2) )
+#endif
 
 namespace theplu {