Package htsjdk.samtools.reference
Class BlockCompressedIndexedFastaSequenceFile
java.lang.Object
htsjdk.samtools.reference.BlockCompressedIndexedFastaSequenceFile
- All Implemented Interfaces:
ReferenceSequenceFile
,Closeable
,AutoCloseable
A block-compressed FASTA file driven by an index for fast lookups.
Supports two interfaces: the ReferenceSequenceFile for old-style, stateful lookups and a direct getter.
-
Constructor Summary
ConstructorDescriptionBlockCompressedIndexedFastaSequenceFile
(String source, SeekableStream in, FastaSequenceIndex index, SAMSequenceDictionary dictionary, GZIIndex gziIndex) Initialize the given indexed fasta sequence file stream.BlockCompressedIndexedFastaSequenceFile
(Path path, FastaSequenceIndex index, GZIIndex gziIndex) -
Method Summary
Modifier and TypeMethodDescriptionvoid
close()
protected SAMSequenceDictionary
Attempts to find and load the sequence dictionary if present.protected static Path
findFastaIndex
(Path fastaFile) protected static Path
findRequiredFastaIndexFile
(Path fastaFile) protected static File
findSequenceDictionary
(File file) Deprecated.use findSequenceDictionary(Path) instead.protected static Path
findSequenceDictionary
(Path fastaPath) Attempts to locate the sequence dictionary file adjacent to the reference fasta file.protected String
Returns the full path to the reference file.getIndex()
protected Path
getPath()
Returns the path to the reference file.getSequence
(String contig) Retrieves the complete sequence described by this contig.Returns the list of sequence records associated with the reference sequence if found otherwise null.protected String
Returns the named source of the reference file.getSubsequenceAt
(String contig, long start, long stop) Gets the subsequence of the contig in the range [start,stop]final boolean
default implementation -- override if index is supportedGets the next sequence if available, or null if not present.protected int
readFromPosition
(ByteBuffer buffer, long position) Reads a sequence of bytes from this sequence file into the given buffer, starting at the given file position.void
reset()
Reset the iterator over the index.protected static void
sanityCheckDictionaryAgainstIndex
(String fastaFile, SAMSequenceDictionary sequenceDictionary, FastaSequenceIndex index) Do some basic checking to make sure the dictionary and the index match.toString()
Returns the full path to the reference file, or the source if no path was specified.
-
Constructor Details
-
BlockCompressedIndexedFastaSequenceFile
- Throws:
FileNotFoundException
-
BlockCompressedIndexedFastaSequenceFile
-
BlockCompressedIndexedFastaSequenceFile
public BlockCompressedIndexedFastaSequenceFile(Path path, FastaSequenceIndex index, GZIIndex gziIndex) -
BlockCompressedIndexedFastaSequenceFile
public BlockCompressedIndexedFastaSequenceFile(String source, SeekableStream in, FastaSequenceIndex index, SAMSequenceDictionary dictionary, GZIIndex gziIndex) Initialize the given indexed fasta sequence file stream.- Parameters:
source
- The named source of the reference file (used in error messages).in
- The input stream to read the fasta file from; should not be decompressed already.index
- The fasta index.dictionary
- The sequence dictionary, or null if there isn't one.gziIndex
- The GZI index; may not be null.
-
-
Method Details
-
readFromPosition
Reads a sequence of bytes from this sequence file into the given buffer, starting at the given file position.- Parameters:
buffer
- the buffer into which bytes are to be transferredposition
- the position to start reading at- Returns:
- the number of bytes read
- Throws:
IOException
- if an I/O error occurs while reading
-
close
- Throws:
IOException
-
findRequiredFastaIndexFile
- Throws:
FileNotFoundException
-
findFastaIndex
-
sanityCheckDictionaryAgainstIndex
protected static void sanityCheckDictionaryAgainstIndex(String fastaFile, SAMSequenceDictionary sequenceDictionary, FastaSequenceIndex index) Do some basic checking to make sure the dictionary and the index match.- Parameters:
fastaFile
- Used for error reporting only.sequenceDictionary
- sequence dictionary to check against the index.index
- index file to check against the dictionary.
-
getIndex
-
nextSequence
Gets the next sequence if available, or null if not present.- Returns:
- next sequence if available, or null if not present.
-
reset
public void reset()Reset the iterator over the index. -
isIndexed
public final boolean isIndexed()default implementation -- override if index is supported- Specified by:
isIndexed
in interfaceReferenceSequenceFile
- Returns:
- true if getSequence and getSubsequenceAt methods are allowed.
-
getSequence
Retrieves the complete sequence described by this contig.- Specified by:
getSequence
in interfaceReferenceSequenceFile
- Parameters:
contig
- contig whose data should be returned.- Returns:
- The full sequence associated with this contig.
-
getSubsequenceAt
Gets the subsequence of the contig in the range [start,stop]- Specified by:
getSubsequenceAt
in interfaceReferenceSequenceFile
- Parameters:
contig
- Contig whose subsequence to retrieve.start
- inclusive, 1-based start of region.stop
- inclusive, 1-based stop of region.- Returns:
- The partial reference sequence associated with this range.
-
findAndLoadSequenceDictionary
Attempts to find and load the sequence dictionary if present. -
findSequenceDictionary
Deprecated.use findSequenceDictionary(Path) instead. -
findSequenceDictionary
Attempts to locate the sequence dictionary file adjacent to the reference fasta file. -
getPath
Returns the path to the reference file. -
getSource
Returns the named source of the reference file. -
getSequenceDictionary
Returns the list of sequence records associated with the reference sequence if found otherwise null.- Specified by:
getSequenceDictionary
in interfaceReferenceSequenceFile
- Returns:
- a list of sequence records representing the sequences in this reference file
-
getAbsolutePath
Returns the full path to the reference file. -
toString
Returns the full path to the reference file, or the source if no path was specified.- Specified by:
toString
in interfaceReferenceSequenceFile
- Overrides:
toString
in classObject
- Returns:
- Reference name, file name, or something other human-readable representation.
-