Fix -max-coverage=1.
authorTom Clegg <tom@tomclegg.ca>
Fri, 18 Sep 2020 13:47:27 +0000 (09:47 -0400)
committerTom Clegg <tom@tomclegg.ca>
Fri, 18 Sep 2020 13:47:27 +0000 (09:47 -0400)
Arvados-DCO-1.1-Signed-off-by: Tom Clegg <tom@tomclegg.ca>

filter.go

index 4560a281a5d43fe0380d784d63bfa1e62f59546d..89aa778f8bf5e0b9c03b9e8bfe45fc0b7fc709c3 100644 (file)
--- a/filter.go
+++ b/filter.go
@@ -36,7 +36,7 @@ func (cmd *filterer) RunCommand(prog string, args []string, stdin io.Reader, std
        inputFilename := flags.String("i", "-", "input `file`")
        outputFilename := flags.String("o", "-", "output `file`")
        maxvariants := flags.Int("max-variants", -1, "drop tiles with more than `N` variants")
-       mincoverage := flags.Float64("min-coverage", 1, "drop tiles with coverage less than `P` across all haplotypes (0 < P ≤ 1)")
+       mincoverage := flags.Float64("min-coverage", 0, "drop tiles with coverage less than `P` across all haplotypes (0 < P ≤ 1)")
        maxtag := flags.Int("max-tag", -1, "drop tiles with tag ID > `N`")
        err = flags.Parse(args)
        if err == flag.ErrHelp {
@@ -138,7 +138,7 @@ func (cmd *filterer) RunCommand(prog string, args []string, stdin io.Reader, std
                }
        }
 
-       if *mincoverage < 1 {
+       if *mincoverage > 0 {
                mincov := int(*mincoverage * float64(len(cgs)*2))
                cov := make([]int, ntags)
                for _, cg := range cgs {