make varcall more tweakable

ExpressionAnalysis · 5d0cad0268d14cd9c111de3509264fceea93f3f9 · earonesty committed Sep 25, 2014 · Sep 25, 2014 · 5d0cad0
1 parent fb9245c
commit 5d0cad0268d14cd9c111de3509264fceea93f3f9
Unified Split

Showing with 2 additions and 4 deletions.

+2 −4 clipper/varcall.cpp
diff --git a/clipper/varcall.cpp b/clipper/varcall.cpp
@@ -1019,7 +1019,7 @@ inline void PileupSummary::Parse(char *line, PileupReads &rds, tidx *adex, char
     }
 
     int meanreadlen = rds.MeanReadLen();
-    int maxdepthbypos = meanreadlen <= 0 ? 10 : min(10, round(10 * artifact_filter * (Depth/meanreadlen)));
+    int maxdepthbypos = meanreadlen <= 0 ? 10 : max(10, round(10.0 * artifact_filter * (Depth/(double)meanreadlen)));
 
     Calls.clear();
 
@@ -1569,12 +1569,12 @@ void VarCallVisitor::Visit(PileupSummary &p) {
 		    fprintf(stderr,"xpos-skip-amp\t%d\n",p.SkipAmp);
 		    fprintf(stderr,"xpos-skip-mapq\t%d\n",p.SkipMinMapq);
 		    fprintf(stderr,"xpos-skip-qual\t%d\n",p.SkipMinQual);
+			exit(0);
         }
 		++SkippedDepth;
 		return;
 	}
 
-
 	int ins_fwd = p.Calls.size() > 6 ? p.Calls[6].fwd : 0;
 	int ins_rev = p.Calls.size() > 6 ? p.Calls[6].rev : 0;
 
@@ -1974,8 +1974,6 @@ void VarCallVisitor::Visit(PileupSummary &p) {
             }
 			exit(0);
 		}
-
-
 	}
 }