fixed a minor problem in the example coming with kprobaln.c

diff --git a/bcftools/vcfutils.pl b/bcftools/vcfutils.pl
index d0b797143cf56fe0a67383ec4f610096eaecf2c0..e2c31078d3204b0c2110aaa9376c999e86d373ff 100755 (executable)
--- a/bcftools/vcfutils.pl
+++ b/bcftools/vcfutils.pl
@@ -155,7 +155,7 @@ Note: This command discards indels. Output: QUAL #non-indel #SNPs #transitions #
        next if (length($t[3]) != 1 || uc($t[3]) eq 'N');
        $t[3] = uc($t[3]); $t[4] = uc($t[4]);
        my @s = split(',', $t[4]);
-       $t[5] = 3 if ($t[5] < 0);
+       $t[5] = 3 if ($t[5] eq '.' || $t[5] < 0);
        next if (length($s[0]) != 1);
        my $hit;
        if ($is_vcf) {

diff --git a/kprobaln.c b/kprobaln.c
index 1a3d7b63b94cbe00a4c41f52a59c510287a2e4ce..c3d05dd1b0fe0d8c99f7a94515edee7d739233e2 100644 (file)
--- a/kprobaln.c
+++ b/kprobaln.c
@@ -285,7 +285,7 @@ int main(int argc, char *argv[])
        iqual = malloc(l_query);
        memset(iqual, q, l_query);
        kpa_par_def.bw = b;
-       kpa_glocal(ref, l_ref, query, l_query, iqual, &kpa_par_def, 0, 0, 0);
+       kpa_glocal(ref, l_ref, query, l_query, iqual, &kpa_par_def, 0, 0);
        free(iqual);
        return 0;
 }