* samtools-0.1.5-5 (r395)

[samtools.git] / bam_pileup.c

diff --git a/bam_pileup.c b/bam_pileup.c
index 3ffd52851e3be64bed628081b2b57d8ca6378b8b..d2d2aa6a8f1890b22344d330d0ca7e612c93b244 100644 (file)
--- a/bam_pileup.c
+++ b/bam_pileup.c
@@ -212,3 +212,19 @@ int bam_plbuf_push(const bam1_t *b, bam_plbuf_t *buf)
        }
        return 0;
 }
+
+int bam_pileup_file(bamFile fp, int mask, bam_pileup_f func, void *func_data)
+{
+       bam_plbuf_t *buf;
+       int ret;
+       bam1_t *b;
+       b = bam_init1();
+       buf = bam_plbuf_init(func, func_data);
+       bam_plbuf_set_mask(buf, mask);
+       while ((ret = bam_read1(fp, b)) >= 0)
+               bam_plbuf_push(b, buf);
+       bam_plbuf_push(0, buf);
+       bam_plbuf_destroy(buf);
+       bam_destroy1(b);
+       return 0;
+}