Package org.broadinstitute.gatk.utils.activeregion

Examples of org.broadinstitute.gatk.utils.activeregion.ActiveRegion


        final DummyActiveRegionWalker walker = new DummyActiveRegionWalker(activeRegions, false);

        final Map<GenomeLoc, ActiveRegion> activeRegionsMap = getActiveRegions(traversal, walker, intervals, bamBuilder.makeTemporarilyBAMFile());

        final ActiveRegion region = activeRegionsMap.values().iterator().next();
        int nReadsExpectedInRegion = 0;

        final Set<String> readNamesInRegion = readNamesInRegion(region);
        for ( final GATKSAMRecord read : bamBuilder.makeReads() ) {
            Assert.assertTrue(readNamesInRegion.contains(read.getReadName()),
                    "Region " + region + " should contain read " + read + " with cigar " + read.getCigarString() + " but it wasn't");
            nReadsExpectedInRegion++;
        }

        Assert.assertEquals(region.size(), nReadsExpectedInRegion, "There are more reads in active region " + region + "than expected");
    }
View Full Code Here

TOP

Related Classes of org.broadinstitute.gatk.utils.activeregion.ActiveRegion

Copyright © 2018 www.massapicom. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.