Refactored wiggle code and added rsem-bam2readdepth program

[rsem.git] / bam2readdepth.cpp

diff --git a/bam2readdepth.cpp b/bam2readdepth.cpp
new file mode 100644 (file)
index 0000000..c7f0adb
--- /dev/null
+++ b/bam2readdepth.cpp
@@ -0,0 +1,13 @@
+#include <iostream>
+#include "wiggle.h"
+
+int main(int argc, char* argv[]) {
+       if (argc != 2) {
+               printf("Usage: rsem-bam2readdepth sorted_bam_input\n");
+        std::exit(1);
+       }
+    ReadDepthWriter depth_writer(std::cout);
+    build_wiggles(argv[1], depth_writer);
+
+       return 0;
+}