Merge pull request #58 from broadinstitute/mc_callset_assesment_GSATDG-52

AGBT scripts, tool updates and misc
This commit is contained in:
depristo 2013-02-22 07:23:59 -08:00
commit 50612ac981
1 changed files with 6 additions and 3 deletions

View File

@ -70,6 +70,9 @@ public class FindCoveredIntervals extends ActiveRegionWalker<GenomeLoc, Long> {
@Output(required = true)
private PrintStream out;
@Argument(fullName = "uncovered", shortName = "u", required = false, doc = "output intervals that fail the coverage threshold instead")
private boolean outputUncovered = false;
@Argument(fullName = "coverage_threshold", shortName = "cov", doc = "The minimum allowable coverage to be considered covered", required = false)
private int coverageThreshold = 20;
@ -86,10 +89,10 @@ public class FindCoveredIntervals extends ActiveRegionWalker<GenomeLoc, Long> {
@Override
public GenomeLoc map(final org.broadinstitute.sting.utils.activeregion.ActiveRegion activeRegion, final RefMetaDataTracker tracker) {
if (activeRegion.isActive())
if ((!outputUncovered && activeRegion.isActive()) || (outputUncovered && !activeRegion.isActive()))
return activeRegion.getLocation();
else
return null;
return null;
}
@Override