Search in sources :

Example 1 with AllelicCountCollector

use of org.broadinstitute.hellbender.tools.copynumber.allelic.alleliccount.AllelicCountCollector in project gatk-protected by broadinstitute.

the class CollectAllelicCounts method doWork.

@Override
protected Object doWork() {
    validateArguments();
    final File referenceFile = referenceArguments.getReferenceFile();
    final IntervalList siteIntervals = IntervalList.fromFile(inputSiteIntervalsFile);
    final AllelicCountCollector allelicCountCollector = new AllelicCountCollector(referenceFile, readValidationStringency);
    logger.info("Collecting allelic counts...");
    final AllelicCountCollection allelicCounts = allelicCountCollector.collect(inputBAMFile, siteIntervals, minimumMappingQuality, minimumBaseQuality);
    allelicCounts.write(outputAllelicCountsFile);
    logger.info("Allelic counts written to " + outputAllelicCountsFile.toString());
    return "SUCCESS";
}
Also used : AllelicCountCollector(org.broadinstitute.hellbender.tools.copynumber.allelic.alleliccount.AllelicCountCollector) IntervalList(htsjdk.samtools.util.IntervalList) AllelicCountCollection(org.broadinstitute.hellbender.tools.copynumber.allelic.alleliccount.AllelicCountCollection) File(java.io.File)

Example 2 with AllelicCountCollector

use of org.broadinstitute.hellbender.tools.copynumber.allelic.alleliccount.AllelicCountCollector in project gatk by broadinstitute.

the class CollectAllelicCounts method doWork.

@Override
protected Object doWork() {
    validateArguments();
    final File referenceFile = referenceArguments.getReferenceFile();
    final IntervalList siteIntervals = IntervalList.fromFile(inputSiteIntervalsFile);
    final AllelicCountCollector allelicCountCollector = new AllelicCountCollector(referenceFile, readValidationStringency);
    logger.info("Collecting allelic counts...");
    final AllelicCountCollection allelicCounts = allelicCountCollector.collect(inputBAMFile, siteIntervals, minimumMappingQuality, minimumBaseQuality);
    allelicCounts.write(outputAllelicCountsFile);
    logger.info("Allelic counts written to " + outputAllelicCountsFile.toString());
    return "SUCCESS";
}
Also used : AllelicCountCollector(org.broadinstitute.hellbender.tools.copynumber.allelic.alleliccount.AllelicCountCollector) IntervalList(htsjdk.samtools.util.IntervalList) AllelicCountCollection(org.broadinstitute.hellbender.tools.copynumber.allelic.alleliccount.AllelicCountCollection) File(java.io.File)

Aggregations

IntervalList (htsjdk.samtools.util.IntervalList)2 File (java.io.File)2 AllelicCountCollection (org.broadinstitute.hellbender.tools.copynumber.allelic.alleliccount.AllelicCountCollection)2 AllelicCountCollector (org.broadinstitute.hellbender.tools.copynumber.allelic.alleliccount.AllelicCountCollector)2