zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added -forceActive argument to ActiveRegionWalkers 

-- Causes the ART tool to treat all bases as active.   Useful for debugging
This commit is contained in:
Mark DePristo 2013-03-28 17:35:00 -04:00
parent 8656bd5e29
commit 167cd49e71
2 changed files with 9 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
public/java/src/org/broadinstitute/sting/gatk/traversals/TraverseActiveRegions.java
View File

@ -487,6 +487,7 @@ public class TraverseActiveRegions<M, T> extends TraversalEngine<M,T,ActiveRegio
final AlignmentContext locus) {
// must be called, even if we won't use the result, to satisfy walker contract
final ActivityProfileState state = walker.isActive( tracker, refContext, locus );
if ( walker.forceActive) state.isActiveProb = 1.0;
if ( ! walkerHasPresetRegions ) {
activityProfile.add(state);
}

8
public/java/src/org/broadinstitute/sting/gatk/walkers/ActiveRegionWalker.java
View File

@ -88,6 +88,14 @@ public abstract class ActiveRegionWalker<MapType, ReduceType> extends Walker<Map
@Argument(fullName="activeRegionExtension", shortName="activeRegionExtension", doc="The active region extension; if not provided defaults to Walker annotated default", required = false)
public Integer activeRegionExtension = null;
/**
* For the active region walker to treat all bases as active. Useful for debugging when you want to force something like
* the HaplotypeCaller to process a specific interval you provide the GATK
*/
@Advanced
@Argument(fullName="forceActive", shortName="forceActive", doc="If provided, all bases will be tagged as active", required = false)
public boolean forceActive = false;
@Advanced
@Argument(fullName="activeRegionMaxSize", shortName="activeRegionMaxSize", doc="The active region maximum size; if not provided defaults to Walker annotated default", required = false)
public Integer activeRegionMaxSize = null;