Package org.broadinstitute.gatk.utils.fasta

Examples of org.broadinstitute.gatk.utils.fasta.CachingIndexedFastaSequenceFile

1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
    }

    @Test(expectedExceptions=UserException.BadArgumentValue.class)
    public void testExceptionUponLegacyIntervalSyntax() throws Exception {
        GenomeAnalysisEngine toolkit = new GenomeAnalysisEngine();
        toolkit.setGenomeLocParser(new GenomeLocParser(new CachingIndexedFastaSequenceFile(new File(BaseTest.hg19Reference))));

        // Attempting to use the legacy -L "interval1;interval2" syntax should produce an exception:
        IntervalBinding<Feature> binding = new IntervalBinding<Feature>("1;2");
        binding.getIntervals(toolkit);
    }
View Full Code Here
220
221
222
223
224
225
226
227
228
229
230
231
    }

    @Test (enabled = true)
    public void testReadWithNsRefIndexInDeletion() throws FileNotFoundException {

        final IndexedFastaSequenceFile seq = new CachingIndexedFastaSequenceFile(new File(b37KGReference));
        final SAMFileHeader header = ArtificialSAMUtils.createArtificialSamHeader(seq.getSequenceDictionary());
        final int readLength = 76;

        final GATKSAMRecord read = ArtificialSAMUtils.createArtificialRead(header, "myRead", 0, 8975, readLength);
        read.setReadBases(Utils.dupBytes((byte) 'A', readLength));
        read.setBaseQualities(Utils.dupBytes((byte)30, readLength));
View Full Code Here
236
237
238
239
240
241
242
243
244
245
246
247
    }

    @Test (enabled = true)
    public void testReadWithNsRefAfterDeletion() throws FileNotFoundException {

        final IndexedFastaSequenceFile seq = new CachingIndexedFastaSequenceFile(new File(b37KGReference));
        final SAMFileHeader header = ArtificialSAMUtils.createArtificialSamHeader(seq.getSequenceDictionary());
        final int readLength = 76;

        final GATKSAMRecord read = ArtificialSAMUtils.createArtificialRead(header, "myRead", 0, 8975, readLength);
        read.setReadBases(Utils.dupBytes((byte) 'A', readLength));
        read.setBaseQualities(Utils.dupBytes((byte)30, readLength));
View Full Code Here
44
45
46
47
48
49
50
    @Param({"true", "false"})
    boolean useContigIndex; // set automatically by framework

    @Override protected void setUp() throws Exception {
        seq = new CachingIndexedFastaSequenceFile(new File("/Users/depristo/Desktop/broadLocal/localData/human_g1k_v37.fasta"));
    }
View Full Code Here
61
62
63
64
65
66
67
68
    private IndexedFastaSequenceFile seq;

    @BeforeClass
    public void setup() throws FileNotFoundException {
        // sequence
        seq = new CachingIndexedFastaSequenceFile(new File(b37KGReference));
        genomeLocParser = new GenomeLocParser(seq);
    }
View Full Code Here
53
54
55
56
57
58
59
60
    private static SAMSequenceDictionary dict;

    @BeforeClass
    public void init() throws FileNotFoundException {
        // sequence
        seq = new CachingIndexedFastaSequenceFile(new File(b37KGReference));
        dict = seq.getSequenceDictionary();
    }
View Full Code Here
62
63
64
65
66
67
68
        ATCATC = Allele.create("ATCATC");
        ATCATCT = Allele.create("ATCATCT");
        ATref = Allele.create("AT",true);
        Anoref = Allele.create("A",false);
        GT = Allele.create("GT",false);
        genomeLocParser = new GenomeLocParser(new CachingIndexedFastaSequenceFile(new File(hg18Reference)));
    }
View Full Code Here
54
55
56
57
58
59
60
61
    private GenomeLocParser genomeLocParser;

    @BeforeClass
    public void init() throws FileNotFoundException {
        // sequence
        seq = new CachingIndexedFastaSequenceFile(new File(hg18Reference));
        genomeLocParser = new GenomeLocParser(seq);
    }
View Full Code Here
56
57
58
59
60
61
62
63
    private IndexedFastaSequenceFile seq;

    @BeforeClass
    public void setup() throws FileNotFoundException {
        // sequence
        seq = new CachingIndexedFastaSequenceFile(new File(b37KGReference));
        genomeLocParser = new GenomeLocParser(seq);
    }
View Full Code Here
351
352
353
354
355
356
357
358
359
360
361
    //
    // -----------------------------------------------------------------------------------------------

    @DataProvider(name = "GetOverlapping")
    public Object[][] makeGetOverlappingTest() throws Exception {
        final GenomeLocParser genomeLocParser = new GenomeLocParser(new CachingIndexedFastaSequenceFile(new File(b37KGReference)));

        List<Object[]> tests = new ArrayList<Object[]>();

        final GenomeLoc prev1 = genomeLocParser.createGenomeLoc("19", 1, 10);
        final GenomeLoc prev2 = genomeLocParser.createGenomeLoc("19", 20, 50);
View Full Code Here
TOP

Related Classes of org.broadinstitute.gatk.utils.fasta.CachingIndexedFastaSequenceFile

  • htsjdk.samtools.reference.ReferenceSequence
  • htsjdk.samtools.SAMSequenceRecord
  • org.broadinstitute.gatk.engine.datasources.providers.ReferenceOrderedViewUnitTest
  • org.broadinstitute.gatk.engine.datasources.providers.ReferenceViewTemplate
  • org.broadinstitute.gatk.engine.datasources.reads.FilePointerUnitTest
  • org.broadinstitute.gatk.engine.datasources.reads.SAMDataSourceUnitTest
  • org.broadinstitute.gatk.engine.datasources.reference.ReferenceDataSource
  • org.broadinstitute.gatk.engine.datasources.rmd.ReferenceOrderedDataPoolUnitTest
  • org.broadinstitute.gatk.engine.datasources.rmd.ReferenceOrderedQueryDataPoolUnitTest
  • org.broadinstitute.gatk.engine.ReadMetricsUnitTest
    • 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.