Examples of SAMFileReader

Examples of htsjdk.samtools.SAMFileReader

73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
    }

    private void loadHeader() {

        InputStream is = null;
        SAMFileReader reader = null;
        try {
            is = ParsingUtils.openInputStreamGZ(new ResourceLocator(samFile));
            BufferedInputStream bis = new BufferedInputStream(is);
            SAMFileReader.setDefaultValidationStringency(ValidationStringency.SILENT);
            reader = new SAMFileReader(bis);
            header = reader.getFileHeader();
        } catch (IOException e) {
            log.error("Error loading header", e);
        } finally {
            try {
                if (is != null) {
                    is.close();
                }
                if (reader != null) {
                    reader.close();
                }

            } catch (Exception e) {

            }
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

119
120
121
122
123
124
125
126
127
128
129
        FeatureIndex.TileDef seekPos = featureIndex.getTileDef(sequence, startTileNumber);

        if (seekPos != null) {
            // Skip to the start of the query interval and open a sam file reader
            SAMFileReader reader = getSAMFileReader(samFile, seekPos.getStartPosition());
            CloseableIterator<SAMRecord> iter = reader.iterator();
            return new SAMQueryIterator(sequence, start, end, contained, iter);
        }
        return EmptyAlignmentIterator.getInstance();
    }
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

159
160
161
162
163
164
165
166
167
        return sequenceNames;

    }

    public CloseableIterator<PicardAlignment> iterator() {
        SAMFileReader reader = getSAMFileReader(samFile, -1);
        CloseableIterator<SAMRecord> iter = reader.iterator();
        return new SAMQueryIterator(iter);
    }
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

170
171
172
173
174
175
176
177
178
179
180
181
        try {
            SeekableStream stream = IGVSeekableStreamFactory.getInstance().getStreamFor(samFile);
            if (startPosition >= 0) {
                stream.seek(startPosition);
            }
            SAMFileReader reader = new SAMFileReader(stream);
            reader.setValidationStringency(ValidationStringency.SILENT);

            //Need to keep the file source, if loading lazily
            //TODO Can't reload from SAM files. See SAMTextReader.getIterator
            //reader.enableFileSource(PicardAlignment.DEFAULT_LAZY_LOAD);
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

44
45
46
47
48
49
50
51
52
53
54
    public SAMFileHeader getHeader() {
        if (samFiles.length == 0) {
            return null;
        }
        SAMFileReader reader = new SAMFileReader(samFiles[0]);
        SAMFileHeader header = reader.getFileHeader();
        reader.close();
        ;
        return header;
    }
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

91
92
93
94
95
96
97
98
99
100
101
                }
                if (currentReader != null) {
                    currentReader.close();
                }
                if (idx < samFiles.length) {
                    currentReader = new SAMFileReader(samFiles[idx]);
                    currentIter = currentReader.iterator();
                    idx++;
                }
            }
            next = currentIter == null ? null : currentIter.next();
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

68
69
70
71
72
73
74
75
76
77
78
79
            SeekableStream ss = new IGVSeekableBufferedStream(IGVSeekableStreamFactory.getInstance().getStreamFor(url), 128000);
            //SeekableStream ss = getSeekableStream(url);
            log.debug("Initializing SAMFileReader");

            reader = new SAMFileReader(ss, indexFile, false);
        } else {
            InputStream is = HttpUtils.getInstance().openConnectionStream(url);
            reader = new SAMFileReader(new BufferedInputStream(is));
        }

    }
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

42
43
44
45
46
47
48
49
50
51
52
    SAMFileHeader header;

    public BAMFileReader(File bamFile) {
        try {
            File indexFile = findIndexFile(bamFile);
            reader = new SAMFileReader(bamFile, indexFile);
            reader.setValidationStringency(ValidationStringency.SILENT);
            loadHeader();
        } catch (Exception e) {
            MessageUtils.showMessage("Error loading SAM header: " + e.getMessage());
            e.printStackTrace();
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

120
121
122
123
124
125
126
127
128
129
130
    public CloseableIterator<PicardAlignment> iterator() {
        try {
            if (reader == null) {
                InputStream is = HttpUtils.getInstance().openConnectionStream(url);
                reader = new SAMFileReader(new BufferedInputStream(is));
            }
            return new WrappedIterator(reader.iterator());
        } catch (IOException e) {
            log.error("Error creating iterator", e);
            throw new RuntimeException(e);
View Full Code Here

Examples of htsjdk.samtools.SAMFileReader

134
135
136
137
138
139
140
141
142
143
144
    public CloseableIterator<PicardAlignment> query(String sequence, int start, int end, boolean contained) {
        try {
            if (reader == null) {
                SeekableStream ss = new IGVSeekableBufferedStream(IGVSeekableStreamFactory.getInstance().getStreamFor(url));
                reader = new SAMFileReader(ss, indexFile, false);
            }
            CloseableIterator<SAMRecord> iter = reader.query(sequence, start + 1, end, contained);
            return new WrappedIterator(iter);
        } catch (IOException e) {
            log.error("Error opening SAM reader", e);
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. 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.