New -p switch for increased sensititivy of indel-calling. Here -m and -F are used...

[samtools.git] / bam2bcf.c

diff --git a/bam2bcf.c b/bam2bcf.c
index 6ac5dce182a2353ada163dd65bdb5ca530d59127..502d832455fb09bd2b9a2e1a1b21d077cafd869c 100644 (file)
--- a/bam2bcf.c
+++ b/bam2bcf.c
@@ -26,6 +26,7 @@ bcf_callaux_t *bcf_call_init(double theta, int min_baseQ)
        bca->e = errmod_init(1. - theta);
        bca->min_frac = 0.002;
        bca->min_support = 1;
+    bca->per_sample_flt = 0;
        return bca;
 }
 
@@ -54,7 +55,6 @@ int bcf_call_glfgen(int _n, const bam_pileup1_t *pl, int ref_base, bcf_callaux_t
                bca->bases = (uint16_t*)realloc(bca->bases, 2 * bca->max_bases);
        }
        // fill the bases array
-       memset(r, 0, sizeof(bcf_callret1_t));
        for (i = n = r->n_supp = 0; i < _n; ++i) {
                const bam_pileup1_t *p = pl + i;
                int q, b, mapQ, baseQ, is_diff, min_dist, seqQ;