Search in sources :

Example 6 with BayesianHetPulldownCalculator

use of org.broadinstitute.hellbender.tools.exome.pulldown.BayesianHetPulldownCalculator in project gatk-protected by broadinstitute.

the class GetBayesianHetCoverage method runNormalOnly.

/**
     * The normal-only workflow
     */
private void runNormalOnly() {
    final BayesianHetPulldownCalculator normalHetPulldownCalculator;
    final Pulldown normalHetPulldown;
    normalHetPulldownCalculator = new BayesianHetPulldownCalculator(REFERENCE_ARGUMENTS.getReferenceFile(), IntervalList.fromFile(snpFile), minimumMappingQuality, minimumBaseQuality, readDepthThreshold, VALIDATION_STRINGENCY, errorProbabilityAdjustmentFactor, new BalancedHeterozygousPileupPriorModel());
    logger.info("Calculating the Het pulldown from the normal BAM file using the BALANCED prior...");
    normalHetPulldown = normalHetPulldownCalculator.getHetPulldown(normalBamFile, hetCallingStringency);
    logger.info("Writing Het pulldown from normal reads to " + normalHetOutputFile.toString());
    normalHetPulldown.write(normalHetOutputFile, AllelicCountTableColumn.AllelicCountTableVerbosity.FULL);
}
Also used : Pulldown(org.broadinstitute.hellbender.tools.exome.pulldown.Pulldown) BalancedHeterozygousPileupPriorModel(org.broadinstitute.hellbender.tools.exome.pulldown.BalancedHeterozygousPileupPriorModel) BayesianHetPulldownCalculator(org.broadinstitute.hellbender.tools.exome.pulldown.BayesianHetPulldownCalculator)

Aggregations

BayesianHetPulldownCalculator (org.broadinstitute.hellbender.tools.exome.pulldown.BayesianHetPulldownCalculator)6 Pulldown (org.broadinstitute.hellbender.tools.exome.pulldown.Pulldown)6 BalancedHeterozygousPileupPriorModel (org.broadinstitute.hellbender.tools.exome.pulldown.BalancedHeterozygousPileupPriorModel)4 HeterogeneousHeterozygousPileupPriorModel (org.broadinstitute.hellbender.tools.exome.pulldown.HeterogeneousHeterozygousPileupPriorModel)2