Java FeatureManager類代碼示例

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
@Override
protected Object doWork() {
    final FeatureCodec<? extends Feature, ?> codec = FeatureManager.getCodecForFile(inputBedFile);
    final Class<? extends Feature> featureType = codec.getFeatureType();
    if (BEDFeature.class.isAssignableFrom(featureType)) {
        final FeatureDataSource<? extends BEDFeature> source = new FeatureDataSource<>(inputBedFile);
        try {
            final List<Target> targets = StreamSupport.stream(source.spliterator(), false).map(ConvertBedToTargetFile::createTargetFromBEDFeature)
                    .collect(Collectors.toList());
            TargetWriter.writeTargetsToFile(outFile, targets);
        } catch (final TribbleException e) {
            throw new UserException.BadInput(String.format("'%s' has a .bed extension but does not seem to be a valid BED file.", inputBedFile.getAbsolutePath()));
        }
    } else {
        throw new UserException.BadInput(String.format("'%s' does not seem to be a BED file.", inputBedFile.getAbsolutePath()));
    }
    return "SUCCESS";
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
@Override
protected Object doWork() {
    if (!featureFile.canRead()) {
        throw new UserException.CouldNotReadInputFile(featureFile);
    }

    // Get the right codec for the file to be indexed. This call will throw an appropriate exception
    // if featureFile is not in a supported format or is unreadable.
    final FeatureCodec<? extends Feature, ?> codec = new ProgressReportingDelegatingCodec<>(FeatureManager.getCodecForFile(featureFile), ProgressMeter.DEFAULT_SECONDS_BETWEEN_UPDATES);

    final Index index = createAppropriateIndexInMemory(codec);
    final File indexFile = determineFileName(index);

    try {
        index.write(indexFile);
    } catch (final IOException e) {
        throw new UserException.CouldNotCreateOutputFile("Could not write index to file " + indexFile.getAbsolutePath(), e);
    }

    logger.info("Successfully wrote index to " + indexFile.getAbsolutePath());
    return indexFile.getAbsolutePath();
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
private void verifyFileType(
        final File resultVCFFile,
        final String outputExtension) {
    final FeatureCodec<? extends Feature, ?> featureCodec = FeatureManager.getCodecForFile(resultVCFFile);

    if (outputExtension.equals(".vcf") ||
        outputExtension.equals(".vcf.bgz") ||
        outputExtension.equals(".vcf.gz") ||
        outputExtension.equals(".tmp"))
    {
        Assert.assertEquals(featureCodec.getClass(), VCFCodec.class,
                "Wrong codec selected for file " + resultVCFFile.getAbsolutePath());
    }
    else if (outputExtension.equals(".bcf")) {
        Assert.assertEquals(featureCodec.getClass(), BCF2Codec.class,
                "Wrong codec selected for file " + resultVCFFile.getAbsolutePath());
    }
    else {
        throw new IllegalArgumentException("Unknown file extension in createVCFWriter test validation");
    }
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
@SuppressWarnings("unchecked")
private static FeatureCodec<VariantContext, ?> getCodecForVariantSource( final String variantSource ) {
    final FeatureCodec<? extends Feature, ?> codec = FeatureManager.getCodecForFile(new File(variantSource));
    if ( !VariantContext.class.isAssignableFrom(codec.getFeatureType()) ) {
        throw new UserException(variantSource + " is not in a format that produces VariantContexts");
    }
    return (FeatureCodec<VariantContext, ?>)codec;
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
@Test(dataProvider = "AnnotateRsIDData")
public void testRsIDFeatureContext(final VariantContext toAnnotate, final List<VariantContext> dbSNPRecords, final String expectedID, final boolean expectOverlap) throws Exception {
    final File file= new File(publicTestDir + "Homo_sapiens_assembly19.dbsnp135.chr1_1M.exome_intervals.vcf");
    final FeatureContext ctx = new FeatureContext(new FeatureManager(new ArtificialFeatureContainingCommandLineProgram_ForVariantOverlap(file)), new SimpleInterval(toAnnotate));
    final VariantOverlapAnnotator annotator = makeAnnotator(file, "dbsnp");
    final VariantContext annotated = annotator.annotateRsID(ctx, toAnnotate);
    Assert.assertNotNull(annotated);
    Assert.assertEquals(annotated, toAnnotate);   //nothing at given position
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
@Test(dataProvider = "AnnotateRsIDData")
public void testRsIDFeatureContextWith2Sets(final VariantContext toAnnotate, final List<VariantContext> dbSNPRecords, final String expectedID, final boolean expectOverlap) throws Exception {
    final File file= new File(publicTestDir + "Homo_sapiens_assembly19.dbsnp135.chr1_1M.exome_intervals.vcf");
    final FeatureContext ctx = new FeatureContext(new FeatureManager(new ArtificialFeatureContainingCommandLineProgram_ForVariantOverlap(file)), new SimpleInterval(toAnnotate));
    final VariantOverlapAnnotator annotator = makeAnnotator(file, "dbsnp", "binding");
    final VariantContext annotated = annotator.annotateRsID(ctx, toAnnotate);
    Assert.assertNotNull(annotated);
    Assert.assertEquals(annotated, toAnnotate);   //nothing at given position
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
@Test(dataProvider = "AnnotateRsIDData")
public void testRsIDFeatureContextWithNoDBSnP(final VariantContext toAnnotate, final List<VariantContext> dbSNPRecords, final String expectedID, final boolean expectOverlap) throws Exception {
    final File file= new File(publicTestDir + "Homo_sapiens_assembly19.dbsnp135.chr1_1M.exome_intervals.vcf");
    final FeatureContext ctx = new FeatureContext(new FeatureManager(new ArtificialFeatureContainingCommandLineProgram_ForVariantOverlap(file)), new SimpleInterval(toAnnotate));
    final VariantOverlapAnnotator annotator = makeAnnotator(file, null, "binding");
    final VariantContext annotated = annotator.annotateRsID(ctx, toAnnotate);
    Assert.assertNotNull(annotated);
    Assert.assertEquals(annotated, toAnnotate);   //nothing at given position
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
@Test(dataProvider = "AnnotateRsIDData")
public void testAnnotateOverlapsFeatureContext(final VariantContext toAnnotate, final List<VariantContext> records, final String expectedID, final boolean expectOverlap) throws Exception {
    final File file= new File(publicTestDir + "Homo_sapiens_assembly19.dbsnp135.chr1_1M.exome_intervals.vcf");
    final FeatureContext ctx = new FeatureContext(new FeatureManager(new ArtificialFeatureContainingCommandLineProgram_ForVariantOverlap(file)), new SimpleInterval(toAnnotate));
    final VariantOverlapAnnotator annotator = makeAnnotator(file, "dbsnp");
    final VariantContext annotated = annotator.annotateOverlaps(ctx, toAnnotate);
    Assert.assertEquals(annotated, toAnnotate);   //nothing at given position
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
@Test(dataProvider = "AnnotateRsIDData")
public void testAnnotateOverlapsFeatureContextWith2Sets(final VariantContext toAnnotate, final List<VariantContext> records, final String expectedID, final boolean expectOverlap) throws Exception {
    final File file= new File(publicTestDir + "Homo_sapiens_assembly19.dbsnp135.chr1_1M.exome_intervals.vcf");
    final FeatureContext ctx = new FeatureContext(new FeatureManager(new ArtificialFeatureContainingCommandLineProgram_ForVariantOverlap(file)), new SimpleInterval(toAnnotate));
    final VariantOverlapAnnotator annotator = makeAnnotator(file, "dbsnp", "binding");
    final VariantContext annotated = annotator.annotateOverlaps(ctx, toAnnotate);
    Assert.assertEquals(annotated, toAnnotate);   //nothing at given position
}
 

import org.broadinstitute.hellbender.engine.FeatureManager; //導入依賴的package包/類
/**
 * Initialize our source of Feature data (or set it to null if no Feature argument(s) were provided).
 *
 * Package-private so that engine classes can access it, but concrete tool child classes cannot.
 * May be overridden by traversals that require custom initialization of Feature data sources.
 *
 * By default, this method initializes the FeatureManager to use the lookahead cache of {@link FeatureDataSource#DEFAULT_QUERY_LOOKAHEAD_BASES} bases.
 */
void initializeFeatures() {
    features = new FeatureManager(this);
    if ( features.isEmpty() ) {  // No available sources of Features discovered for this tool
        features = null;
    }
}