FiberHMM
Chromatin footprint calling on single DNA molecules from Fiber-seq (m6A) and
DAF-seq (deamination) data. For every read, FiberHMM calls nucleosomes,
transcription-factor and other sub-nucleosomal footprints, and the
methylase-sensitive patches (MSPs) between them, using sequence-context-aware
emission models; across molecules, it finds recurrent footprint classes and
measures how often each molecule carries them. Calls are written into the BAM
as fibertools-compatible ns/nl/as/al and Molecular-annotation MA/AQ
tags, readable by ft, FIRE and FiberBrowser.
Documentation: fiberseq.github.io/FiberHMM
(source in docs/).
FiberHMM 3.0.0 is the third generation of the FiberHMM family, released with FiberBrowser 3.0.0. Nanopore Hia5 and DddB users: the Nanopore Hia5 and DddB emission tables shipped in 2.x releases were context-swapped; re-run those calls made with 2.x.
fiberhmm-check <your outputs>lists which files need re-running. See Upgrading from 2.x and the CHANGELOG.
Install
pip install fiberhmm # Python >= 3.10
pip install "fiberhmm[all]" # plus QC figures (matplotlib) and HDF5 posteriors (h5py)
Quick start
# DAF-seq amplicons / plasmids: reads + reference in, FiberBrowser-ready calls out
fiberhmm-pipeline reads.fastq.gz --reference plasmid.dna --enzyme dddb -o run/
# Aligned BAMs
fiberhmm-call -i aligned.bam -o calls.bam --enzyme hia5 --seq pacbio -c 8 --region-parallel # Hia5 PacBio
fiberhmm-call -i aligned.bam -o calls.bam --enzyme hia5 --seq nanopore -c 8 --region-parallel # Hia5 Nanopore
fiberhmm-call -i aligned.bam -o calls.bam --enzyme dddb -c 8 --region-parallel # DddB DAF-seq (DddA: --enzyme ddda)
fiberhmm-extract -i calls.bam -o tracks/ # BED12 / bigBed tracks
A synthetic demo data set and a walk-through of every chemistry are in the Quick start.
Documentation
- Getting started:
installation, quick start, choosing
--enzymeand--seq - Concepts: how FiberHMM works, chemistries and platforms, coordinate frames, tags and scores
- Workflows:
calling, re-calling, DAF-seq, duplex pairing, tracks, QC, footprint classes
(
fiberhmm-consensus), transfer, strand rescue, training - Reference: every command and option, BAM tags, header declarations, bundled models, consensus output files, Python API
- Troubleshooting
License
MIT; see LICENSE. Citing FiberHMM.
Metadata
Release files for fiberhmm 3.0.0
For a detailed explanation of source distributions (sdists) and built distributions (wheels), please see the package formats documentation.
Source distribution (sdist)
| File | Size | Uploaded | |
|---|---|---|---|
| fiberhmm-3.0.0.tar.gz | 5.7 MB | Details |
Built distribution (wheel)
| File | Interpreter | ABI | Platform | Reset |
|---|---|---|---|---|
| fiberhmm-3.0.0-py3-none-any.whl | Python 3 | none | any | Details |
Total release size: 8.9 MB
Release files / fiberhmm-3.0.0.tar.gz
| Download URL | fiberhmm-3.0.0.tar.gz |
|---|---|
| Size | 5.7 MB |
| Tags | Source |
|
SHA-256 checksum How to use checksums |
ef512ce2a32b40e071efea26fcb410b157f5347226c94956dae1db0af6631c19
|
|
BLAKE2b-256 checksum How to use checksums |
671441ae5a8f0e900c3c24ebb1599e08334a21bc67acf800f05d788c124ad6f5
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/6.2.0 CPython/3.10.13
|
Release files / fiberhmm-3.0.0-py3-none-any.whl
| Download URL | fiberhmm-3.0.0-py3-none-any.whl |
|---|---|
| Size | 3.2 MB |
| Tags | Python 3 |
|
SHA-256 checksum How to use checksums |
3a3e8184b28c8569447be901c90b60348dc26b5eb2bde9110c330e154249b8c0
|
|
BLAKE2b-256 checksum How to use checksums |
0f43be709d916ca3774c06d75d6665117ccb6fa0bade3db4fcde2e18c071eb3b
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
twine/6.2.0 CPython/3.10.13
|