diff --git a/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/cancer/m2/MuTect2.java b/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/cancer/m2/MuTect2.java index 4722e51f3..335e13487 100644 --- a/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/cancer/m2/MuTect2.java +++ b/protected/gatk-tools-protected/src/main/java/org/broadinstitute/gatk/tools/walkers/cancer/m2/MuTect2.java @@ -200,8 +200,6 @@ public class MuTect2 extends ActiveRegionWalker, Integer> i private double log10GlobalReadMismappingRate; private static final int REFERENCE_PADDING = 500; private static final byte MIN_TAIL_QUALITY_WITH_ERROR_CORRECTION = 6; - private final static int maxReadsInRegionPerSample = 1000; // TODO -- should be an argument - private final static int minReadsPerAlignmentStart = 5; // TODO -- should be an argument @ArgumentCollection protected M2ArgumentCollection MTAC = new M2ArgumentCollection(); @@ -293,6 +291,16 @@ public class MuTect2 extends ActiveRegionWalker, Integer> i @Argument(fullName="doNotRunPhysicalPhasing", shortName="doNotRunPhysicalPhasing", doc="Disable physical phasing", required = false) protected boolean doNotRunPhysicalPhasing = false; + /** + * When downsampling, level the coverage of the reads in each sample to no more than maxReadsInRegionPerSample reads, + * not reducing coverage at any read start to less than minReadsPerAlignmentStart + */ + @Argument(fullName = "maxReadsInRegionPerSample", shortName = "maxReadsInRegionPerSample", doc="Maximum reads in an active region", required = false) + protected int maxReadsInRegionPerSample = 1000; + + @Argument(fullName = "minReadsPerAlignmentStart", shortName = "minReadsPerAlignStart", doc="Minimum number of reads sharing the same alignment start for each genomic location in an active region", required = false) + protected int minReadsPerAlignmentStart = 5; + public RodBinding getDbsnpRodBinding() { return MTAC.dbsnp.dbsnp; } @Override diff --git a/public/gatk-engine/src/main/java/org/broadinstitute/gatk/engine/traversals/TraverseActiveRegions.java b/public/gatk-engine/src/main/java/org/broadinstitute/gatk/engine/traversals/TraverseActiveRegions.java index 33b2d9dff..af3558f7f 100644 --- a/public/gatk-engine/src/main/java/org/broadinstitute/gatk/engine/traversals/TraverseActiveRegions.java +++ b/public/gatk-engine/src/main/java/org/broadinstitute/gatk/engine/traversals/TraverseActiveRegions.java @@ -162,8 +162,8 @@ public final class TraverseActiveRegions extends TraversalEngine extends Walker