pyhmmer

Added

• Sequence.L property to get the length of a sequence consistently with HMM.M.
• Matrix.flatten method to get a Vector over the whole C-contiguous buffer of a Matrix.
• allocM property to Profile and OptimizedProfile to get the allocated capacity.
• HMM.emit_sequence, HMM.emit_alignment and Profile.emit_sequence to generate sequences from a model (#102).
• Missing type hints for gap_index, gap_symbol and type properties of Alphabet.

Removed

• Profile.transition_scores due to memory leak.

Fixed

• SequenceFile.__init__ crashing on opening SSI indices for gzip-compressed files.

Added

• Support for Windows AMD64 builds with MinGW (#92, see documentation).
• SequenceBlock.write method to write all sequences from a block to a file.
• SequenceBlock.total_length to compute the sum of lengths of sequences in a block.
• Subclasses DNA, RNA and AA of pyhmmer.easel.Alphabet to allow marking the alphabet type with type annotations.
• VectorI and MatrixI classes to pyhmmer.easel to store C int types (usually 32 bits).
• OptimizedProfile.ssv_filter method for NEON target in addition to SSE.
• msv_filter method to Profile and OptimizedProfile classes.
• Background.null1 to compute the null1 lod score for an arbitrary Sequence.
• HMMFile._file exposing the internal file-like object given to HMMFile constructor, if any (#89).
• SSIReader.primary_keys exposing a read-only sequence of primary keys inside a SSI file.
• SequenceFile.indexed and MSAFile.indexed to access data indexed with a SSI index (#85).

Changed

• Setup compilation in Stable ABI mode for Python 3.12 and later.
• Setup distribution to generate a CMake package file to facilitate usage with downstream scikit-build-core dependent packages.
• breaking: Make most textual attributes and properties str rather than bytes (#88).
  • Turn the following properties of Sequence into str instead of bytes: accession, description, name, source.
  • Turn the following properties of MSA into str instead of bytes: accession, description, name, source, author, names, reference, model_mask, secondary_structure, surface_accessibility, posterior_probabilities.
  • Turn the following properties of HMM, Profile and OptimizedProfile into str instead of bytes: name, accession, description.
  • Turn the following properties of Alignment into str instead of bytes: hmm_name, hmm_accession, hmm_sequence, target_name, target_sequence.
  • Make SequenceBlock.indexed and MSA.indexed mapping use str for keys instead of bytes.
  • Make SSIReader and SSIWriter use str for keys and aliases instead of bytes.
• Make DigitalSequence, DigitalSequenceBlock, DigitalMSA, SequenceFile, MSAFile, HMM, Profile, OptimizedProfileBlock, Background, Builder, Pipeline, HMMFile, HMMPressedFile generic over the alphabet type.
• Reorganize setup so that Easel and HMMER libraries, C headers and Cython headers are installed to the Python site path.
• Reorganize tests to skip tests in the absence of data files and remove larger test files from PyPI distributions.
• Enforce detection of invalid buffer sizes in _from_raw_bytes class constructors of various objects.
• Use PyUnicode_DecodeASCII to decode ASCII strings of known length instead of PyUnicode_FromString where applicable.
• Allow passing a buffer-protocol object interpreted as an ASCII string to the TextSequence constructor.

Fixed

• Remove hmmlogo.c from compiled HMMER sources to avoid an issue with the presence of a main function when linking.
• Avoid using multiprocessing.Value in pyhmmer.hmmer while in single-threaded mode for improved compatibility.
• Memory leak in Matrix caused by outdated allocation logic for zero dimensions.
• Make HMMFile constructor pretend to stream when given a file-like object to prevent HMMER from calling fseek and ftell.
• Incorrect declaration of Cython type declarations in various MSA property setters.
• Ensure correct background alphabet in Pipeline constructor when given a Background as argument.
• TopHits.merge potentially recomputing inclusion/reporting flags for hits obtained with bitscore cuttoffs.
• SSIReader missing methods in type stubs (#85).

Removed

• Outdated compatibility code for PyPy 3.6.
• Support for positional arguments other than alphabet in TextMSA, DigitalMSA, TextSequence and DigitalSequence constructors.

Added

• Support for Windows AMD64 builds with MinGW (#92, see documentation).
• SequenceBlock.write method to write all sequences from a block to a file.
• SequenceBlock.total_length to compute the sum of lengths of sequences in a block.
• Subclasses DNA, RNA and AA of pyhmmer.easel.Alphabet to allow marking the alphabet type with type annotations.
• VectorI and MatrixI classes to pyhmmer.easel to store C int types (usually 32 bits).
• OptimizedProfile.ssv_filter method for NEON target in addition to SSE.
• msv_filter method to Profile and OptimizedProfile classes.
• Background.null1 to compute the null1 lod score for an arbitrary Sequence.
• HMMFile._file exposing the internal file-like object given to HMMFile constructor, if any (#89).

Changed

• Setup compilation in Stable ABI mode for Python 3.12 and later.
• Setup distribution to generate a CMake package file to facilitate usage with downstream scikit-build-core dependent packages.
• breaking: Make most textual attributes and properties str rather than bytes (#88).
  • Turn the following properties of Sequence into str instead of bytes: accession, description, name, source.
  • Turn the following properties of MSA into str instead of bytes: accession, description, name, source, author, names, reference, model_mask, secondary_structure, surface_accessibility, posterior_probabilities.
  • Turn the following properties of HMM, Profile and OptimizedProfile into str instead of bytes: name, accession, description.
  • Turn the following properties of Alignment into str instead of bytes: hmm_name, hmm_accession, hmm_sequence, target_name, target_sequence.
  • Make SequenceBlock.indexed and MSA.indexed mapping use str for keys instead of bytes.
  • Make SSIReader and SSIWriter use str for keys and aliases instead of bytes.
• Make DigitalSequence, DigitalSequenceBlock, DigitalMSA, SequenceFile, MSAFile, HMM, Profile, OptimizedProfileBlock, Background, Builder, Pipeline, HMMFile, HMMPressedFile generic over the alphabet type.
• Reorganize setup so that Easel and HMMER libraries, C headers and Cython headers are installed to the Python site path.
• Reorganize tests to skip tests in the absence of data files and remove larger test files from PyPI distributions.
• Enforce detection of invalid buffer sizes in _from_raw_bytes class constructors of various objects.
• Use PyUnicode_DecodeASCII to decode ASCII strings of known length instead of PyUnicode_FromString where applicable.
• Allow passing a buffer-protocol object interpreted as an ASCII string to the TextSequence constructor.

Fixed

• Remove hmmlogo.c from compiled HMMER sources to avoid an issue with the presence of a main function when linking.
• Avoid using multiprocessing.Value in pyhmmer.hmmer while in single-threaded mode for improved compatibility.
• Memory leak in Matrix caused by outdated allocation logic for zero dimensions.
• Make HMMFile constructor pretend to stream when given a file-like object to prevent HMMER from calling fseek and ftell.
• Incorrect declaration of Cython type declarations in various MSA property setters.
• Ensure correct background alphabet in Pipeline constructor when given a Background as argument.
• TopHits.merge potentially recomputing inclusion/reporting flags for hits obtained with bitscore cuttoffs.
• SSIReader missing methods in type stubs (#85).

Removed

• Outdated compatibility code for PyPy 3.6.
• Support for positional arguments other than alphabet in TextMSA, DigitalMSA, TextSequence and DigitalSequence constructors.

Added

• Support for Windows AMD64 builds with MinGW (#92, see documentation).
• SequenceBlock.write method to write all sequences from a block to a file.
• SequenceBlock.total_length to compute the sum of lengths of sequences in a block.
• Subclasses DNA, RNA and AA of pyhmmer.easel.Alphabet to allow marking the alphabet type with type annotations.
• VectorI and MatrixI classes to pyhmmer.easel to store C int types (usually 32 bits).
• OptimizedProfile.ssv_filter method for NEON target in addition to SSE.
• msv_filter method to Profile and OptimizedProfile classes.
• Background.null1 to compute the null1 lod score for an arbitrary Sequence.
• HMMFile._file exposing the internal file-like object given to HMMFile constructor, if any (#89).

Changed

• Setup compilation in Stable ABI mode for Python 3.12 and later.
• Setup distribution to generate a CMake package file to facilitate usage with downstream scikit-build-core dependent packages.
• breaking: Make most textual attributes and properties str rather than bytes (#88).
  • Turn the following properties of Sequence into str instead of bytes: accession, description, name, source.
  • Turn the following properties of MSA into str instead of bytes: accession, description, name, source, author, names, reference, model_mask, secondary_structure, surface_accessibility, posterior_probabilities.
  • Turn the following properties of HMM, Profile and OptimizedProfile into str instead of bytes: name, accession, description.
  • Turn the following properties of Alignment into str instead of bytes: hmm_name, hmm_accession, hmm_sequence, target_name, target_sequence.
  • Make SequenceBlock.indexed and MSA.indexed mapping use str for keys instead of bytes.
  • Make SSIReader and SSIWriter use str for keys and aliases instead of bytes.
• Make DigitalSequence, DigitalSequenceBlock, DigitalMSA, SequenceFile, MSAFile, HMM, Profile, OptimizedProfileBlock, Background, Builder, Pipeline, HMMFile, HMMPressedFile generic over the alphabet type.
• Reorganize setup so that Easel and HMMER libraries, C headers and Cython headers are installed to the Python site path.
• Reorganize tests to skip tests in the absence of data files and remove larger test files from PyPI distributions.
• Enforce detection of invalid buffer sizes in _from_raw_bytes class constructors of various objects.
• Use PyUnicode_DecodeASCII to decode ASCII strings of known length instead of PyUnicode_FromString where applicable.
• Allow passing a buffer-protocol object interpreted as an ASCII string to the TextSequence constructor.

Fixed

• Remove hmmlogo.c from compiled HMMER sources to avoid an issue with the presence of a main function when linking.
• Avoid using multiprocessing.Value in pyhmmer.hmmer while in single-threaded mode for improved compatibility.
• Memory leak in Matrix caused by outdated allocation logic for zero dimensions.
• Make HMMFile constructor pretend to stream when given a file-like object to prevent HMMER from calling fseek and ftell.
• Incorrect declaration of Cython type declarations in various MSA property setters.
• Ensure correct background alphabet in Pipeline constructor when given a Background as argument.
• TopHits.merge potentially recomputing inclusion/reporting flags for hits obtained with bitscore cuttoffs.
• SSIReader missing methods in type stubs (#85).

Removed

• Outdated compatibility code for PyPy 3.6.
• Support for positional arguments other than alphabet in TextMSA, DigitalMSA, TextSequence and DigitalSequence constructors.

Added

• SequenceBlock.write method to write all sequences from a block to a file.
• SequenceBlock.total_length to compute the sum of lengths of sequences in a block.
• Subclasses DNA, RNA and AA of pyhmmer.easel.Alphabet to allow marking the alphabet type with type annotations.
• VectorI and MatrixI classes to pyhmmer.easel to store C int types (usually 32 bits).
• OptimizedProfile.ssv_filter method for NEON target in addition to SSE.
• msv_filter method to Profile and OptimizedProfile classes.
• Background.null1 to compute the null1 lod score for an arbitrary Sequence.

Changed

• Setup compilation in Stable ABI mode for Python 3.12 and later.
• Setup distribution to generate a CMake package file to facilitate usage with downstream scikit-build-core dependent packages.
• breaking: Make most textual attributes and properties str rather than bytes (#88).
  • Turn the following properties of Sequence into str instead of bytes: accession, description, name, source.
  • Turn the following properties of MSA into str instead of bytes: accession, description, name, source, author, names, reference, model_mask, secondary_structure, surface_accessibility, posterior_probabilities.
  • Turn the following properties of HMM, Profile and OptimizedProfile into str instead of bytes: name, accession, description.
  • Turn the following properties of Alignment into str instead of bytes: hmm_name, hmm_accession, hmm_sequence, target_name, target_sequence.
  • Make SequenceBlock.indexed and MSA.indexed mapping use str for keys instead of bytes.
  • Make SSIReader and SSIWriter use str for keys and aliases instead of bytes.
• Make DigitalSequence, DigitalSequenceBlock, DigitalMSA, SequenceFile, MSAFile, HMM, Profile, OptimizedProfileBlock, Background, Builder, Pipeline, HMMFile, HMMPressedFile generic over the alphabet type.
• Reorganize setup so that Easel and HMMER libraries, C headers and Cython headers are installed to the Python site path.
• Reorganize tests to skip tests in the absence of data files and remove larger test files from PyPI distributions.
• Enforce detection of invalid buffer sizes in _from_raw_bytes class constructors of various objects.
• Use PyUnicode_DecodeASCII to decode ASCII strings of known length instead of PyUnicode_FromString where applicable.
• Allow passing a buffer-protocol object interpreted as an ASCII string to the TextSequence constructor.

Fixed

• Remove hmmlogo.c from compiled HMMER sources to avoid an issue with the presence of a main function when linking.
• Avoid using multiprocessing.Value in pyhmmer.hmmer while in single-threaded mode for improved compatibility.
• Memory leak in Matrix caused by outdated allocation logic for zero dimensions.
• Make HMMFile constructor pretend to stream when given a file-like object to prevent HMMER from calling fseek and ftell.
• Incorrect declaration of Cython type declarations in various MSA property setters.
• Ensure correct background alphabet in Pipeline constructor when given a Background as argument.
• TopHits.merge potentially recomputing inclusion/reporting flags for hits obtained with bitscore cuttoffs.
• SSIReader missing methods in type stubs (#85).

Removed

• Outdated compatibility code for PyPy 3.6.
• Support for positional arguments other than alphabet in TextMSA, DigitalMSA, TextSequence and DigitalSequence constructors.

Added

• SequenceBlock.write method to write all sequences from a block to a file.
• SequenceBlock.total_length to compute the sum of lengths of sequences in a block.
• Subclasses DNA, RNA and AA of pyhmmer.easel.Alphabet to allow marking the alphabet type with type annotations.
• VectorI and MatrixI classes to pyhmmer.easel to store C int types (usually 32 bits).
• OptimizedProfile.ssv_filter method for NEON target in addition to SSE.

Changed

• Setup compilation in Stable ABI mode for Python 3.12 and later.
• breaking: Make most textual attributes and properties str rather than bytes (#88).
  • Turn the following properties of Sequence into str instead of bytes: accession, description, name, source.
  • Turn the following properties of MSA into str instead of bytes: accession, description, name, source, author, names, reference, model_mask, secondary_structure, surface_accessibility, posterior_probabilities.
  • Turn the following properties of HMM, Profile and OptimizedProfile into str instead of bytes: name, accession, description.
  • Turn the following properties of Alignment into str instead of bytes: hmm_name, hmm_accession, hmm_sequence, target_name, target_sequence.
  • Make SequenceBlock.indexed and MSA.indexed mapping use str for keys instead of bytes.
  • Make SSIReader and SSIWriter use str for keys and aliases instead of bytes.
• Make DigitalSequence, DigitalSequenceBlock, DigitalMSA, SequenceFile, MSAFile, HMM, Profile, OptimizedProfileBlock, Background, Builder, Pipeline, HMMFile, HMMPressedFile generic over the alphabet type.
• Reorganize setup so that Easel and HMMER libraries, C headers and Cython headers are installed to the Python site path.
• Reorganize tests to skip tests in the absence of data files and remove larger test files from PyPI distributions.
• Enforce detection of invalid buffer sizes in _from_raw_bytes class constructors of various objects.
• Use PyUnicode_DecodeASCII to decode ASCII strings of known length instead of PyUnicode_FromString where applicable.
• Allow passing a buffer-protocol object interpreted as an ASCII string to the TextSequence constructor.

Fixed

• Remove hmmlogo.c from compiled HMMER sources to avoid an issue with the presence of a main function when linking.
• Avoid using multiprocessing.Value in pyhmmer.hmmer while in single-threaded mode for improved compatibility.
• Memory leak in Matrix caused by outdated allocation logic for zero dimensions.
• Make HMMFile constructor pretend to stream when given a file-like object to prevent HMMER from calling fseek and ftell.
• Incorrect declaration of Cython type declarations in various MSA property setters.
• Ensure correct background alphabet in Pipeline constructor when given a Background as argument.

Removed

• Outdated compatibility code for PyPy 3.6.
• Support for positional arguments other than alphabet in TextMSA, DigitalMSA, TextSequence and DigitalSequence constructors.

Added

• SequenceBlock.write method to write all sequences from a block to a file.
• Subclasses DNA, RNA and AA of pyhmmer.easel.Alphabet to allow marking the alphabet type with type annotations.

Changed

• Setup compilation in Stable ABI mode for Python 3.12 and later.
• breaking: Make most textual attributes and properties str rather than bytes (#88).
  • Turn the following properties of Sequence into str instead of bytes: accession, description, name, source.
  • Turn the following properties of MSA into str instead of bytes: accession, description, name, source, author, names, reference, model_mask, secondary_structure, surface_accessibility, posterior_probabilities.
  • Turn the following properties of HMM, Profile and OptimizedProfile into str instead of bytes: name, accession, description.
  • Turn the following properties of Alignment into str instead of bytes: hmm_name, hmm_accession, hmm_sequence, target_name, target_sequence.
  • Make SequenceBlock.indexed and MSA.indexed mapping use str for keys instead of bytes.
  • Make SSIReader and SSIWriter use str for keys and aliases instead of bytes.
• Make DigitalSequence, DigitalSequenceBlock, DigitalMSA, SequenceFile, MSAFile, HMM, Profile, OptimizedProfileBlock, Background, Builder, Pipeline, HMMFile, HMMPressedFile generic over the alphabet type.
• Reorganize setup so that Easel and HMMER libraries, C headers and Cython headers are installed to the Python site path.
• Reorganize tests to skip tests in the absence of data files and remove larger test files from PyPI distributions.
• Enforce detection of invalid buffer sizes in _from_raw_bytes class constructors of various objects.
• Use PyUnicode_DecodeASCII to decode ASCII strings of known length instead of PyUnicode_FromString where applicable.

Fixed

• Remove hmmlogo.c from compiled HMMER sources to avoid an issue with the presence of a main function when linking.
• Avoid using multiprocessing.Value in pyhmmer.hmmer while in single-threaded mode for improved compatibility.
• Memory leak in Matrix caused by outdated allocation logic for zero dimensions.
• Make HMMFile constructor pretend to stream when given a file-like object to prevent HMMER from calling fseek and ftell.
• Incorrect declaration of Cython type declarations in various MSA property setters.

Removed

• Outdated compatibility code for PyPy 3.6.
• Support for positional arguments other than alphabet in TextMSA, DigitalMSA, TextSequence and DigitalSequence constructors.

Added

• SequenceBlock.write method to write all sequences from a block to a file.
• Subclasses DNA, RNA and AA of pyhmmer.easel.Alphabet to allow marking the alphabet type with type annotations.

Changed

• breaking: Make most textual attributes and properties str rather than bytes (#88).
  • Turn the following properties of Sequence into str instead of bytes: accession, description, name, source.
  • Turn the following properties of MSA into str instead of bytes: accession, description, name, source, author, names, reference, model_mask, secondary_structure, surface_accessibility, posterior_probabilities.
  • Turn the following properties of HMM, Profile and OptimizedProfile into str instead of bytes: name, accession, description.
  • Turn the following properties of Alignment into str instead of bytes: hmm_name, hmm_accession, hmm_sequence, target_name, target_sequence.
  • Make SequenceBlock.indexed and MSA.indexed mapping use str for keys instead of bytes.
  • Make SSIReader and SSIWriter use str for keys and aliases instead of bytes.
• Make DigitalSequence, DigitalSequenceBlock, DigitalMSA, SequenceFile, MSAFile, HMM, Profile, OptimizedProfileBlock, Background, Builder, Pipeline, HMMFile, HMMPressedFile generic over the alphabet type.
• Reorganize setup so that Easel and HMMER libraries, C headers and Cython headers are installed to the Python site path.
• Reorganize tests to skip tests in the absence of data files and remove larger test files from PyPI distributions.
• Enforce detection of invalid buffer sizes in _from_raw_bytes class constructors of various objects.
• Use PyUnicode_DecodeASCII to decode ASCII strings of known length instead of PyUnicode_FromString where applicable.

Fixed

• Remove hmmlogo.c from compiled HMMER sources to avoid an issue with the presence of a main function when linking.
• Avoid using multiprocessing.Value in pyhmmer.hmmer while in single-threaded mode for improved compatibility.
• Memory leak in Matrix caused by outdated allocation logic for zero dimensions.
• Make HMMFile constructor pretend to stream when given a file-like object to prevent HMMER from calling fseek and ftell.
• Incorrect declaration of Cython type declarations in various MSA property setters.

Removed

• Outdated compatibility code for PyPy 3.6.
• Support for positional arguments other than alphabet in TextMSA, DigitalMSA, TextSequence and DigitalSequence constructors.

Added

• Support for multithreading in hmmalign.

Fixed

• Formatting issues in TestMSA.compute_weights documentation (#96).
• Trace.__eq__ exiting the interpreter on inequality (EddyRivasLab/hmmer#344).

Changed

• Make Traces more similar to SequenceBlock and implement list-like interface (append, extend, clear, ...).

Added

• VectorD and MatrixD classes for double-precision data.
• __abs__ implementations for floating-point Vector and Matrix types.
• MSA.sequence_weights to access the sequence weights of an alignment.
• MSA.compute_weights to compute the sequence weights of an alignment like esl-weight (#94)
• MSA.mark_fragments method to identify fragmented sequences in an alignment.
• DigitalMSA.reverse_complement method to reverse-complement a nucleotide alignment.

Fixed

• Erroneous categories in documentation (#95, by @mtrsl).
• TextMSA and DigitalMSA constructors now default-initialize the sequence weights of the Easel alignment.

Changed

• Relax psutil dependencies to allow versions past 5.0.
• Replace some cases of ValueError exceptions in pyhmmer.easel with more specific InvalidParameter errors.

Added

• Support for Python 3.14.

Fixed

• Issue with HMMER manually settings denormals-are-zero flags on x86-64 with SSE4.1 enabled (#71).
• Missing documentation of multiple easel class constructors.

Changed

• Improve documentation of supported MSA formats (#93).

Added

• DigitalMSA.identity_filter method to remove sequences too similar from a multiple sequence alignment (#84).
• reference, model_mask, secondary_structure, surface_accessibility, posterior_probabilities properties of MSA to get and set additional column annotation for a sequence alignment.
• MSA.select to select a subset of columns and rows of a sequence alignment given an iterable of indices.
• New parallelization strategy in hmmsearch suitable for low query counts to parallelize on target chunks rather than individual queries (more similar to original HMMER).
• TextSequence.sample and DigitalSequence.sample constructors to generate random sequences for testing.
• DigitalMSA.sample constructor to generate a random multiple sequence alignment for testing.
• MSA.indexed and SequenceBlock.indexed to get the sequences of an alignment or block by name.
• DigitalMSA.alignment to get the rows of an alignment in digital mode as VectorU8 objects.

Changed

• Remove some test data from pyhmmer.tests.data to reduce distribution size.
• Use fused types to specialize queries of Pipeline search methods.
• Improve error reporting message in Builder.build_msa on input format errors.
• Deprecate positional arguments in TextSequence, DigitalSequence, TextMSA and DigitalMSA constructors.
• Ensure to poison workers before yielding the last result in pyhmmer.hmmer dispatchers.
• Allow passing a seed integer or None to HMM.sample.
• Use pure-Python collections.abc.Sequence classes to expose MSA.sequences instead of Cython class.
• Use pure-Python collections.abc.Sequence classes to expose TextMSA.alignment instead of generating a tuple on demand.
• Relax psutil dependency to support both 6.0 and 7.0.

Fixed

• Deployment to Arch User Repository.
• Use of deprecated API in documentation.
• Pipeline methods not raising errors on unsupported targets (#87).
• Matrix.__getitem__ raising IndexError when slicing with a slice without upper bound.
• TopHits.merge not properly merging book-keeping attributes when merging TopHits with zero hits.
• Relax equality comparison in TopHits.merge when the query was an OptimizedProfile.
• Some uncaught connection errors in pyhmmer.hmmer background workers.
• Broken assertion in TextMSA.sequences due to Easel bug in esl_sq_FetchFromMSA (EddyRivasLab/easel#80).
• Fused Cython methods not being properly documented in Sphinx documentation.

Added

• Missing type annotations to specific options of LongTargetsPipeline.
• __len__ implementation for HMMPressedFile using the entries in the SSI index.
• Setters for residue-wise annotation properties in HMM.
• Support for process-based parallelism in pyhmmer.hmmer in addition to thread-based.
• Read-only buffer protocol implementation for TextSequence and DigitalSequence classes.

Changed

• Drop support for Python 3.6.
• Use CMake and scikit-build-core to build the package instead of setuptools.
• Use TypedDict API to mark allowed keyword arguments in pyhmmer.hmmer functions.
• Reorganize detection and handling of alphabets for arbitrary queries in pyhmmer.hmmer.
• Allow passing a SequenceFile to most pyhmmer.hmmer functions.
• Use faster PyUnicode_FromStringAndSize function to decode strings of known lengths in several plan7 classes.
• Make SequenceFile and MSAFile generic on the individual sequence and MSA types.
• Reorganize pyhmmer.hmmer into different submodules.

Fixed

• Type annotations not using typing-extensions for Python3.8 to 3.10.
• Logic error causing out-of-bounds memory access in TopHits.__getstate__.
• Avoid creating a new Pipeline object when running single-threaded searches in pyhmmer.hmmer.
• Detect the appropriate SequenceFile type based on the digital flag value (#72).

Removed

• Deprecated properties of TopHits (query_name, query_length, query_accession).

Added

• query property TopHits referencing the original object used to create the TopHits #76.

Changed

• Require the query object to create a TopHits object.
• Make TopHits generic over its query property.
• Deprecate old query properties of TopHits (query_name, query_length, query_accession).

Removed

• Detection of SSE flush from setup.py (#71).

Added

• Detection of SSE flush modes to setup.py for possible performance gains on x86 platforms.

Changed

• Migrate documentation to pydata-sphinx-theme.

Fixed

• Documentation examples not using permanent resource links.

Changed

• Allow AlphabetMismatch error to allow for an unknown actual alphabet.
• Make HMMFile and HMMPressedFile raise AlphabetMismatch on files with mixed alphabets.

Fixed

• Avoid calling fclose with null pointers in Sequence.write and MSA.write.

Fixed

• HMM.__setstate__ not properly extracting the cutoff from pickle state for some HMMs (#67).

Changed

• Update and remove some test files to reduce size of distributed package data.

Fixed

• Compilation of Easel and HMMER code not using SSE4.1 extensions.

Fixed

• Implement write function for fopencookie with off_t instead of off64_t for compatibility.
• Fix handling of NULL buffers passed to read and write methods of fopencookie.

Fixed

• Reallocation issue causing segmentation faults in nhmmer with more than 64 sequences (#62).

Added

• Getter to access the strand of a Domain produced by a LongTargetsPipeline.

Changed

• Display model and cutoff names in MissingCutoffs error message, if any.
• Allow LongTargetsPipeline to be configured with window length and beta parameters.
• Make nhmmer use the window length and beta from the options when creating a Builder.

Fixed

• nhmmer not computing E-values for non-default window lengths (moshi4/pybarrnap#2).
• SequenceFile and MSAFile crashing with a segmentation fault when given the path to a folder rather than a file.

Added

• Pre-compiled wheels for PyPy 3.10.

Fixed

• Invalid pointer cast in __getbuffer__ method of Matrix and Vector objects.
• Remaining tests failing to run on missing importlib-resources.
• pyhmmer.hmmer dispatchers possibly dead-locking on background thread errors (#60).

Added

• armv7 and aarch64 to the PKGBUILD architectures.

Changed

• SSIReader and SSIWriter constructors now accept path-like objects.
• Skip tests dependending on importlib.resources.files when it is not available on the host machine.

Fixed

• Memory leak caused by alphabet allocation in Pipeline._scan_loop_file.

Added

• Alignment properties to get the original lengths of the sequence and HMM being stored.
• Hit.length property storing the length of the hit sequence (or HMM).
• TopHits.query_length storing the length of the hit HMM (or query).
• Alignment.posterior_probabilities property showing an encoded representation of posteriors (#59, by @arajkovic).
• Trace.score method to compute a trace score from a given profile and sequence.
• Alignment.__sizeof__ implementation leveraing p7_alidisplay_SizeOf.

Fixed

• Cutoffs proxy objects not recording their owner to prevent deallocation.
• Avoid GIL re-acquisition in GeneticCode.translate.
• Query metadata not being recorded in Hits obtained from daemon.Client.
• Empty MatrixU8 creation attempting zero-allocation.
• VectorU8.zeros allocating 4x more memory than required.
• Memory leak caused by string duplication in __getbuffer__ methods of Matrix and Vector types.

Added

• residue_markups argument to TextSequence and DigitalSequence constructors.
• __reduce__ implementation to TextSequence, DigitalSequence, TextSequenceBlock and DigitalSequenceBlock.

Changed

• Handling of easel I/O methods to avoid implicit GIL acquisition for error checking.

Fixed

• Syntax errors in type annotation files.

Added

• Out-of-band pickle serialization of Bitfield objects.
• Getters for float attributes and forward/backward parameters of OptimizedProfile.
• InvalidHMM error raised by HMM.validate.

Changed

• Mark HMM.zero method as noexcept.
• Increase size of buffer for the query queue in the hmmer dispatcher.

Fixed

• Unneeded semaphore in pyhmmer.hmmer message passing implementation.
• Broken assertion in Bitfield._from_raw_bytes.
• Relax tolerance of HMM validation in TraceAligner.align_traces.

Fixed

• Invalid buffer write in DigitalSequenceBlock.translate (#50).

Added

• HMM.set_consensus method to set the consensus for a method or compute it from the emission probabilities.

Fixed

• Platform detection for MacOS and Armv7 platforms in setup.py.
• pyhmmer.plan7.HMM constructor setting a consensus string forcefully.

Added

• Support for compiling wheels for Aarch64 and NEON-enabled Arm platforms.

Changed

• Updated HMMER to v3.4.
• Updated Easel to v0.49.
• Use cibuildwheel to build wheel distributions.

Fixed

• Patch missing PyInterpreterState_GetID preventing the package from working on PyPy 3.9.

Added

• TopHits.mode property showing from which pipeline mode (search or scan) the hits were obtained.

Changed

• Updated the code for Cython v3.0.

Fixed

• TopHits.merge not properly handling inclusion and reporting for domains (#46, #47, by @zdk123).

Added

• Bracket-style repr implementation to HMM, Profile and OptimizedProfile showing model alphabet, length and name.
• MissingCutoffs and InvalidParameter exceptions inheriting ValueError.

Changed

• Replace pthread locks with PyThread API for synchronizing models in OptimizedProfileBlock.

Fixed

• Sequence length extraction in LongTargetsPipeline.search_hmm (#42).
• LongTargetsPipeline.search_msa not building a HMM with Builder.build_msa.

Added

• HMM.validate method to ensure a HMM holds HMMER structural constraints.
• plan7.Transitions enum with transition names for indexing HMM.transition_probabilities.

PyHMMER has been accepted for publication in Bioinformatics. Paper can be reached at doi:10.1093/bioinformatics/btad214.

Added

• pyhmmer.hmmer.jackhmmer function to run several JackHMMER iterative searches in parallel using multithreading (#35, by @zdk123).
• HMM.to_profile shortcut method to allocate and configure a new Profile object.

Fixed

• Type annotations of Pipeline.iterate_seq and Pipeline.iterate_hmm.
• Potential memory leak on exceptions raised by HMMPressedFile.read.
• Offsets.profile not recording offsets properly, causing pyhmmer.hmmer.hmmpress to produce invalid pressed files (#37).

Changed

• HMM.__init__ and HMM.sample now take the Alphabet as the first argument, for consistency with the rest of the API.
• HMM now require a name argument.

Removed

• Deprecated ignore_gaps argument in SequenceFile.__init__.
• Deprecated Sequence.taxonomy_id property.

Added

• Recipes page to the documentation with code example for loading multiple HMM files (#24, by @zdk123).

Fixed

• TraceAligner methods causing a segfault when passed an uninitialized HMM (#36).

Changed

• HMM default constructor now always creates a valid HMM (with respects to probability arrays).
• TraceAligner now validates the input HMM before calling the HMMER code.
• Use stack allocation for all error buffers instead of creating empty bytearray objects where applicable.

Fixed

• Wrong argument type in IterativeSearch.iterate_hmm method (#34, by @zdk123).

Added

• easel.GeneticCode class wrapping an ESL_GENCODE struct for configuring translation.
• DigitalSequence.translate method to translate a nucleotide sequence to a protein sequence. Metadata is copied from the source sequence to its translation (#31, by @valentynbez).

Deprecated

• Sequence.taxonomy_id property, as it is not used by Easel and implementation is not consistent (see EddyRivasLab/easel#68).

Added

• Missing __reduce__ method to TopHits.

Fixed

• Build detection of available platform functions in setup.py.

Added

• Bitfield.zeros and Bitfield.ones classmethods for constructing an empty bitfield of known size.
• Bitfield.copy method to copy a bitfield object.
• SequenceBlock and OptimizedProfileBlock classes to store Python objects next to a contiguous array of pointers for iterating with the GIL released.
• SequenceFile.read_block method to read a whole sequence block from a file.
• HMM.sample class method to generate a HMM at random given a Randomness source.
• hmmscan function to scan a profile database with sequence queries.
• deepcopy implementations to HMM, Profile and OptimizedProfile classes of plan7.
• rewind method to HMMFile, HMMPressedFile and SequenceFile to reset a file back to its initial position.
• name attribute to HMMFile, HMMPressedFile, MSAFile and SequenceFile to expose the path of a file (when it was created from path).
• local property to Profile and OptimizedProfile, indicating whether a profile is in local or global mode.
• multihit property to Profile and OptimizedProfile, indicating whether a profile is in unihit or multihit mode, with a setter taking care of the reconfiguration.
• Domain.included and Domain.reported settable properties to report the inclusion and reporting status of a single domain.
• TopHits.included and TopHits.reported sized iterator to iterate only on included and reported hits.
• Domains.included and Domains.reported sized iterator to iterate only on included and reported domains.

Changed

• Bitfield, Vector and Matrix can now be created from an iterable.
• Pipeline search methods now expect a DigitalSequenceBlock or a SequenceFile for the target sequence database.
• Pipeline scan methods now expect an OptimizedProfileBlock or a HMMPressedFile for the target profile database.
• TraceAligner now expect a DigitalSequenceBlock for the sequences to align to the HMM.
• Profile.configure now uses a default value of 400 for the L argument.
• hmmsearch, nhmmer and phmmer support being given a single query instead of requiring an iterable.
• HMMPressedFile can now be created, closed and used as a context manager directly without having to manage the source HMMFile.
• Renamed Profile.optimized method to Profile.to_optimized.
• Replaced Randomness.is_fast method with the Randomness.fast property.
• Rewrite handling of Hit flags using settable properties (Hit.included, Hit.reported, Hit.new, Hit.dropped, Hit.duplicate) instead of methods.

Fixed

• Memory leak in the LongTargetsPipeline search loop.
• PyPy behaviour change of readinto methods now expecting unsigned char* instead of char* memoryview.
• NULL-pointer dereference in Pipeline.search_hmm when given a query without name.
• LongTargetsPipeline not recording the query name and accession.
• Memory leak caused by using a non-default prior scheme when constructing a Builder.

Removed

• PipelineSearchTargets, replaced in functionality with easel.DigitalSequenceBlock.
• is_local and is_multihit methods of Profile and OptimizedProfile, replaced with equivalent properties.
• Hit.manually_drop and Hit.manually_include methods, replaced with the different Hit properties.

Fixed

• Error not being raised on alphabet detection failure in SequenceFile or MSAFile.
• Add check in DigitalSequence constructor to make sure encoded characters are in valid range (#25).

Added

• SequenceFile.guess_alphabet and MSAFile.guess_alphabet to guess the alphabet from an open file.
• Alphabet.encode and Alphabet.decode to convert raw sequences between digital and text format.

Changed

• hmmsearch, phmmer and nhmmer functions will reduce the requested number of threads to the number of queries, if it can be detected using operator.length_hint.

Added

• Documentation for loading all HMMs from an HMMFile object at once (#23).
• List of projects depending on PyHMMER to the Examples page of the documentation.

Added

• pickle protocol support for TopHits objects, using the HMMER network serialization.
• TopHits.write method to write hits to a file in tabular format.
• query_name and query_accession properties to TopHits objects to access the name and accession of the query that produced the hits.

Fixed

• Extraction of filename from file-like objects in the HMMFile constructor.
• Use os.cpu_count instead of multiprocessing.cpu_count where applicable to preserve OS scheduling.
• Wrong return type in docstring of HMM.insert_emissions.
• TopHits.searched_nodes returning the searched number of residues instead of the searched number of model nodes.
• Unsound decoding of pickled MatrixF or VectorF when data comes from a source of different endianness.

Changed

• Rewrite pyhmmer.hmmer threading code using Deque instead of collections.Queue to store the queries and results.
• Reduce memory consumption of pyhmmer.hmmer by reducing the number of semaphores and event flags used concurrently.
• Make pyhmmer.hmmer main threads block on query insertion rather than result retrieval to make sure worker threads are never idling.

Added

• pyhmmer.daemon module with an client implementation to communicate to a hmmpgmd server.
• Pipeline.arguments methods to get a list of CLI arguments from the parameters used to initialize the Pipeline.
• Setters for name, accession and description properties of plan7.Hit.
• Constructor for individual plan7.Trace objects outside a plan7.Traces list.
• plan7.Trace.from_sequence constructor to create a faux trace from a single sequence.
• manually_include and manually_drop methods to plan7.Hit for manually selecting the inclusion status of a Hit in a TopHits instance.
• compare_ranking method to plan7.TopHits for comparing the order of the hits compared to a previous run on the same targets stored in an easel.KeyHash object.
• Pipeline.iterate_seq and Pipeline.iterate_hmm to run iterative queries like JackHMMER.
• repr implementations for easel.MSAFile, easel.SequenceFile and easel.HMMFile showing the path or file object they were created from.
• repr implementation for easel.Randomness showing the seed and the RNG algorithm in use.
• str implementation for plan7.Alignment using HMMER original code to display a domain alignment like in search/scan results.

Changed

• plan7.Trace.posterior_probabilities property may now be None in case no memory is allocated for the posteriors in the P7_TRACE struct.
• TopHits.to_msa can now add additional sequences passed as arguments to the alignment.
• plan7.HMMPressedFile now raises an exception on attempts to create a new instance manually.
• ignore_gaps argument of easel.SequenceFile is now deprecated.
• repr implementations for easel types now use the fully qualified class name.

Fixed

• easel.SequenceFile.readinto docstring not rendering properly in documentation.
• Type annotations of hits_included and hits_reported of plan7.TopHits marking these properties as bool instead of int.
• Setters of name, accession, description and author properties of easel.MSA crashing when given None values.
• Exception value raised from Easel code not being properly extracted.
• Plain strings being used in example for easel.TextSequence and easel.TextMSA constructors where byte strings are expected (#20).

Added

• plan7.PipelineSearchTargets to reduce the overhead when searching the same sequences several times with different. query profiles.
• TopHits.copy method to duplicate a TopHits instance.
• TopHits.merge method to merge hits obtained with the same query on different targets.
• Buffer protocol implementation for pyhmmer.easel.Bitfield.

Changed

• Renamed TopHits.included and TopHits.reported properties to TopHits.hits_included and TopHits.hits_included.
• MSAFile and SequenceFile are now directly in digital mode if they are instantiated with digital=True.
• SequenceFile.parse can now return a sequence in digital mode.
• Reorganized tests to make then runnable from a site install.

Fixed

• Usage of memcpy in contexts where it may have had undefined behaviour.
• VectorF.__eq__ crashing when comparing two empty objects.
• SequenceFile and MSAFile not closing file handles when raising an error in __init__.

Added

• plan7.HMMFile.read method to read a single plan7.HMM from an plan7.HMMFile (instead of using next).
• closed property on easel.SequenceFile, easel.MSAFile and plan7.HMMFile to mark whether a file object is closed.
• plan7.HMMFile.is_pressed method to check whether a HMM file has associated pressed data.
• plan7.HMMFile.optimized_profiles methods to read the plan7.OptimizedProfile entries in an plan7.HMMFile is there are associated pressed data available.
• Getters for the name, accession, description, consensus, consensus_structure, evalue_parameters and cutoffs properties of a plan7.OptimizedProfile.
• plan7.OptimizedProfile.__eq__ implementation to compare two optimized profiles.
• __sizeof__ implementations for plan7.OptimizedProfile and plan7.Profile to get the allocated size of a profile.

Fixed

• Double-free caused by the Cython cycle breaking feature on several view types (easel.Randomness, easel.Vector, easel.Matrix, plan7.Cutoffs, plan7.EvalueParameters, plan7.Offsets, plan7.Trace)
• plan7.Hit.description using the pointer to the accession string erroneously, causing occasional NULL dereference.
• plan7.OptimizedProfile.copy performing a shallow copy instead of a deep copy as expected.

Changed

• pyhmmer.hmmer type annotations now explicit support for plan7.Profile or plan7.OptimizedProfile inputs where applicable.

Added

• entropy and relative_entropy methods to easel.VectorF to compute the Shannon entropy of a vector and the Kullback-Leibler divergence of two vectors.
• mean_match_entropy, mean_match_information and mean_match_relative_entropy methods to plan7.HMM to get information statistics of an HMM model.
• match_occupancy method to plan7.HMM to compute the occupancy for each match state as an easel.VectorF.

Fixed

• plan7.Builder.build_msa using the gap-open and gap-extend probabilities instead of the MSA itself to compute the transition probabilities for the new HMM.

Changed

• plan7.Builder.build will now only load the score system once and reuse it unless a different score system is requested between calls.

Added

• plan7.ScoreData class to store the substitution scores and maximal extensions for a long target search.
• plan7.LongTargetsPipeline to run searches on targets longer than 100,000 residues.
• Alphabet methods to check whether an Alphabet object is a DNA, RNA, nucleotide or protein alphabet.
• window_length and window_beta arguments to plan7.Builder to set the max length of nucleotide HMM created by builder objects.

Changed

• pyhmmer.hmmer.nhmmer now uses a LongTargetsPipeline instead of a Pipeline to search the target sequences.
• pyhmmer.hmmer.nhmmer now supports HMM queries in addition to DigitalSequence and DigitalMSA queries.
• pyhmmer.hmmer.phmmer now always assumes protein queries.
• Z and domZ attributes of plan7.TopHits objects is now read-only.

Fixed

• nhmmer now uses DNA as the default alphabet instead of amino acid alphabet like it did before (#12).

Added

• Constructor arguments and properties to plan7.Pipeline to support bit score thresholds instead to filter top hits.
• Support for creating a SequenceFile and an MSAFile using a Python file-like object instead of only supporting filenames.
• Support for reading individual sequences from an MSA file with SequenceFile.
• TextMSA.alignment to access the actual alignment as a tuple of strings.
• Subtraction and division support for easel.Vector subclasses

Changed

• plan7.Cutoffs now support setting the bit score cutoffs, but requires both to be set or cleared at the same time.
• easel.Vector will always allocate some memory when created manually to avoid having a special empty case in every vector method.
• pyhmmer.easel.AllocationError now stores the size it failed to allocate, and the number of elements when allocating an array.

Fixed

• TextSequence.digitize will not raise a ValueError when the sequence contains invalid characters for the alphabet (previously was an UnexpectedError).

Added

• TraceAligner, Trace and Traces classes to pyhmmer.plan7 to get tracebacks after aligning several sequences against an HMM.
• pyhmmer.hmmalign function with the same features as the hmmalign binary from HMMER3.
• Support for out-of-band pickling in easel.Vector and easel.Matrix.

Changed

• Allow creating an empty Vector or Matrix by calling their constructor without arguments.

Fixed

• Potential unreported exceptions in plan7.OptimizedProfile.write and several plan7.SSIWriter methods.

Added

• pickle protocol for easel.Alphabet, easel.Bitfield, easel.KeyHash, easel.Vector, easel.Matrix and plan7.HMM.
• taxonomy_id and residue_markups properties to easel.Sequence.
• sum_score property to plan7.Hit.
• plan7.EvalueParameters class to expose the e-value parameters of a plan7.HMM or a plan7.Profile.
• Equality checks and slicing for easel.Matrix and easel.Vector.
• Support for creating and manipulating zero-sized easel matrices and vectors.
• plan7.Cutoffs class to expose the Pfam score cutoffs of a plan7.HMM or a plan7.Profile.
• Keyword arguments to configure E-value thresholds when creating a plan7.Pipeline object.
• Support for using model-specific thresholding options in plan7.Pipeline.

Changed

• Use the replace error handler when decoding error messages to skip potential decoding issues when already building an exception.
• Improve pyhmmer.hmmer to ensure background threads exit on a KeyboardInterrupt.
• easel.VectorU8.__eq__ accepts any object implementing the buffer protocol.
• plan7.HMM.creation_time now takes and returns a datetime.datetime object, assuming the field is only ever set with asctime.
• Refactor easel.Vector and easel.Matrix and mark exposed memory as C-contiguous.

Fixed

• easel.Alphabet not reporting potential allocation errors.
• Potential buffer overflow in easel.Matrix and easel.Vector when calling __init__ more than once.

Added

• OptimizedProfile.convert method to configure an optimized profile from a Profile without reallocating a new P7_OPROFILE struct.

Changed

• Rewrite the plan7.Pipeline search loop to avoid reacquiring the GIL between reference sequences.
• Require the reference sequences to be stored in a collection (instead of an iterable) when passing them to the search_hmm, search_msa and search_seq methods of plan7.Pipeline.
• Avoid reallocating a new OptimizedProfile every time a new HMM is passed to Pipeline.search_hmm.
• Relax the GIL while sorting and thresholding TopHits in Pipeline search methods.

Added

• ignore_gaps parameter to pyhmmer.plan7.SequenceFile, allowing to skip the gap characters when reading a sequence from an ungapped format.
• __sizeof__ implementation for some
• Dedicated check for sequence length before running the platform-specific code in pyhmmer.plan7.Pipeline.

Fixed

• Score system not being set in pyhmmer.plan7.Builder.build_msa.
• Alphabet not being checked after the first sequence in Pipeline search and scan methods.

Fixed

• File object wrappers not reporting exceptions raised when seeking on OSX/BSD platforms.

Added

• pyhmmer.easel.Randomness class exposing a deterministic random number generator.
• pyhmmer.plan7.Builder.randomness and pyhmmer.plan7.Pipeline.randomness attributes exposing the internal random number generator used by each object.
• pyhmmer.plan7.Hit.best_domain property mapping to the highest scoring domain of a hit.
• pyhmmer.plan7.OptimizedProfile.rbv property exposing match scores.
• pyhmmer.plan7.Domain.pvalue and pyhmmer.plan7.Hit.pvalue reporting the p-value for a domain or hit bitscore.

Fixed

• Dimensions of the pyhmmer.plan7.OptimizedProfile.sbv matrix not being properly set.

Fixed

• Main buffer not being freed in MatrixF.__dealloc__ and MatrixU8.__dealloc__ when created without owner.

Added

• Additional configuration values for pyhmmer.plan7.Pipeline as both constructor arguments and mutable properties.
• consensus, consensus_structure and offsets properties to pyhmmer.plan7.Profile objects.

Changed

• Make OptimizedProfile.ssv_filter check the alphabet of the given sequence.

Added

• Linear algebra primitives to expose 1D (Vector) and 2D (Matrix) contiguous buffers containing numerical values to pyhmmer.easel.
• Documentation for the Z and domZ parameters of the pyhmmer.plan7.Pipeline constructor.
• pyhmmer.errors.AlphabetMismatch exception deriving from ValueError to specifically report mismatching Easel alphabets where applicable.
• scale and normalize methods to pyhmmer.plan7.HMM objects.
• Property to access pyhmmer.plan7.Background residue frequencies as a VectorF object.
• Property to access pyhmmer.plan7.HMM mean residue composition as a VectorF object.
• Property to access pyhmmer.plan7.HMM probabilities and emissions as MatrixF objects.
• ssv_filter methods to pyhmmer.plan7.OptimizedProfile to get the SSV filter score of the profile for a given sequence.
• Several additional properties to access the pyhmmer.plan7.OptimizedProfile internals.

Removed

• Unused report_e parameter of pyhmmer.plan7.Pipeline constructor.
• pyhmmer.plan7.TopHits.clear method which could lead to segfault if it was called while a Hit is being held.

Changed

• Multithreaded loop in pyhmmer.hmmer to reduce memory consumption while still yielding hits in order.
• pyhmmer.easel.DigitalSequence.sequence property is now a VectorU8.

Fixed

• Type annotations in pyhmmer.hmmer.
• Potential double free in pyhmmer.plan7.HMM.command_line property setter.
• Minor floating-point precision issues in pyhmmer.plan7.Builder constructor.
• Segfault in TextMSA.digitize caused by esl_msa_Copy not digitizing on-the-fly like esl_sq_Copy.
• Exceptions not being raised in some methods of pyhmmer.plan7.Profile and pyhmmer.plan7.TopHits.

Added

• Pipeline.scan_seq method to query a database of profiles with one or more sequences.
• transition_probabilities, match_emissions, insert_emissions properties to the HMM class, providing access to the numerical parameters of the HMM.
• consensus_structure and consensus_accessibility properties to the HMM class to get consensus lines from the source alignment if the HMM was created from a MSA.
• nseq and nseq_effective properties to the HMM class to get the number of training sequences and effective sequences used to build the HMM.

Changed

• HMM.checksum is now None if the p7H_CHKSUM flag is not set.
• Builder methods will now record sys.argv when creating a HMM.

Fixed

• HMM.write(..., binary=False) crashing on HMMs without a consensus line. (#5). Fixed upstream in (EddyRivasLab/HMMER#236).
• Pipeline.reset mishandling the Z and domZ values if those were detected from the number of targets.
• pyhmmer.hmmer functions will not block until all results have been collected anymore when run in multithreaded mode.

Added

• easel.MSAFile to read from a file containing
• accession, author, name and description properties to easel.MSA objects.
• plan7.Builder.build_msa to build a pHMM from a sequence alignment.
• Additional methods to easel.KeyHash, allowing to use it as a dict/set hybrid.
• Sequence.write and MSA.write methods to format a sequence or an alignment to a file handle.
• plan7.TopHits.to_msa method to convert all the top hits of a query against a database into a multiple sequence alignment.
• easel.MSA.sequences attribute to access individual sequences of an alignment using the collections.abc.Sequence interface.
• easel.DigitalMSA.textize method to convert a multiple sequence alignment in digital mode to its text-mode counterpart.
• Read-only name, accession and description properties to plan7.Profile showing attributes inherited from the HMM it was configured with.
• plan7.HMM.consensus property, allowing to access the consensus sequence of a pHMM.
• plan7.HMM equality implementation, using zero tolerance.
• plan7.Pipeline.search_msa to query a MSA against a sequence database.
• easel.Sequence.reverse_complement method allowing to reverse-complement inplace or to build a copy.
• errors.AlphabetMismatch exception for use in cases where an alphabet is expected but not matched by the input.
• hmmer.nhmmer function with the same behaviour as hmmer.phmmer, except it expects inputs with a DNA alphabet.

Fixed

• plan7.Builder.copy not copying some parameters correctly, causing pyhmmer.hmmer.phmmer to give inconsistent results in multithreaded mode.
• easel.Bitfield not properly handling index overflows.
• Documentation not rendering for the __init__ method of all classes.

Changed

• plan7.Builder gap-open and gap-extend probabilities are now set on instantiation and depend on the alphabet type.
• Constructors for easel.TextMSA and easel.DigitalMSA, which can now be given an iterable of easel.Sequence objects to store in the alignment.

Removed

• Unimplemented easel.SequenceFile.fetch and easel.SequenceFile.fetchinto methods.