CENTIPEDE applies a hierarchical Bayesian mixture model to infer regions of the genome that are bound by particular transcription factors. It starts by identifying a set of candidate binding sites, and then aims to classify the sites according to whether each site is bound or not bound by a transcription factor. CENTIPEDE is an unsupervised learning algorithm that discriminates between two different types of motif instances using as much relevant information as possible.

This package provides basic routines for estimation of gene-specific transcriptional derivatives and visualization of the resulting velocity patterns.

BWA-PSSM is a probabilistic short genomic sequence read aligner based on the use of position specific scoring matrices (PSSM). Like many of the existing aligners it is fast and sensitive. Unlike most other aligners, however, it is also adaptible in the sense that one can direct the alignment based on known biases within the data set. It is coded as a modification of the original BWA alignment program and shares the genome index structure as well as many of the command line options.

This package provides a fast and accurate analysis toolkit for single cell ATAC-seq (Assay for transposase-accessible chromatin using sequencing). Single cell ATAC-seq can resolve the heterogeneity of a complex tissue and reveal cell-type specific regulatory landscapes. However, the exceeding data sparsity has posed unique challenges for the data analysis. This package r-snapatac is an end-to-end bioinformatics pipeline for analyzing large- scale single cell ATAC-seq data which includes quality control, normalization, clustering analysis, differential analysis, motif inference and exploration of single cell ATAC-seq sequencing data.

Express Beta Diversity (EBD) calculates ecological beta diversity (dissimilarity) measures between biological communities. EBD implements a variety of diversity measures including those that make use of phylogenetic similarity of community members.

Pybiomart provides a simple pythonic interface to biomart.

NanoSV is a software package that can be used to identify structural genomic variations in long-read sequencing data, such as data produced by Oxford Nanopore Technologies’ MinION, GridION or PromethION instruments, or Pacific Biosciences RSII or Sequel sequencers.

The alignment module of BioJava provides an API that contains

• implementations of dynamic programming algorithms for sequence alignment;

• reading and writing of popular alignment file formats;

• a single-, or multi- threaded multiple sequence alignment algorithm.

Change-O is a collection of tools for processing the output of V(D)J alignment tools, assigning clonal clusters to immunoglobulin (Ig) sequences, and reconstructing germline sequences.

This package provides a convenient interface to minimap2, a fast and accurate C program to align genomic and transcribe nucleotide sequences.

This package provides a package that makes it easy to implement sankey, alluvial and sankey bump plots in ggplot2.

This package is a client implementation of the GA4GH htsget protocol. It provides a simple and reliable way to retrieve genomic data from servers supporting the protocol.

Minimap2 is a versatile sequence alignment program that aligns DNA or mRNA sequences against a large reference database. Typical use cases include:

1. mapping PacBio or Oxford Nanopore genomic reads to the human genome;

2. finding overlaps between long reads with error rate up to ~15%;

3. splice-aware alignment of PacBio Iso-Seq or Nanopore cDNA or Direct RNA reads against a reference genome;

4. aligning Illumina single- or paired-end reads;

5. assembly-to-assembly alignment;

6. full-genome alignment between two closely related species with divergence below ~15%.

Circe is a Python package for inferring co-accessibility networks from single-cell ATAC-seq data, using skggm for the graphical lasso and python-scanpy for data processing.

Proteinortho is a tool to detect orthologous genes across different species. For doing so, it compares similarities of given gene sequences and clusters them to find significant groups. The algorithm was designed to handle large-scale data and can be applied to hundreds of species at once.

Biopython is a set of tools for biological computation including parsers for bioinformatics files into Python data structures; interfaces to common bioinformatics programs; a standard sequence class and tools for performing common operations on them; code to perform data classification; code for dealing with alignments; code making it easy to split up parallelizable tasks into separate processes; and more.

Kraken is a taxonomic sequence classifier that assigns taxonomic labels to DNA sequences. Kraken examines the k-mers within a query sequence and uses the information within those k-mers to query a database. That database maps k-mers to the lowest common ancestor (LCA) of all genomes known to contain a given k-mer.

This package provides a toolbox to process, analyze and visualize spatial single-cell expression data.

Bismark is a program to map bisulfite treated sequencing reads to a genome of interest and perform methylation calls in a single step. The output can be easily imported into a genome viewer, such as SeqMonk, and enables a researcher to analyse the methylation levels of their samples straight away. Its main features are:

• Bisulfite mapping and methylation calling in one single step

• Supports single-end and paired-end read alignments

• Supports ungapped and gapped alignments

• Alignment seed length, number of mismatches etc are adjustable

• Output discriminates between cytosine methylation in CpG, CHG and CHH context

PhenoGraph is a clustering method designed for high-dimensional single-cell data. It works by creating a graph representing phenotypic similarities between cells and then identifying communities in this graph.

Segemehl is software to map short sequencer reads to reference genomes. Segemehl implements a matching strategy based on enhanced suffix arrays (ESA). It accepts fasta and fastq queries (gzip'ed and bgzip'ed). In addition to the alignment of reads from standard DNA- and RNA-seq protocols, it also allows the mapping of bisulfite converted reads (Lister and Cokus) and implements a split read mapping strategy. The output of segemehl is a SAM or BAM formatted alignment file.

ctxcore is part of the SCENIC suite of tools. It provides core functions for pycisTarget and SCENIC.

Harmony is an algorithm for integrating multiple high-dimensional datasets with fuzzy k-means and locally linear adjustments.

This package stores information and scripts to work with, or load various different spatial datasets that can be used with the Giotto Suite workflow. It also contains example mini Giotto subobjects.