ACEMD is a high performance molecular dynamics code for biomolecular systems designed specifically for NVIDIA GPUs. Simple and fast, ACEMD uses very similar commands and input files of NAMD and output files as NAMD or Gromacs.

AMBER (Assisted Model Building with Energy Refinement) is a package of molecular simulation programs.

The Accelerating Therapeutics for Opportunites in Medicine (ATOM) Consortium Modeling PipeLine for Drug Discovery. AMPL is an open-source, modular, extensible software pipeline for building and sharing models to advance in silico drug discovery.

APBS (Adaptive Poisson-Boltzmann Solver) is a software package for the numerical solution of the Poisson-Boltzmann equation (PBE), one of the most popular continuum models for describing electrostatic interactions between molecular solutes in salty, aqueous media.

Autodock is a suite of automated docking tools. It is designed to predict how small molecules, such as substrates or drug candidates, bind to a receptor of known 3D structure.

AutoDock Vina is a program for drug discovery, molecular docking and virtual screening, offering multi-core capability, high performance and enhanced accuracy and ease of use. It is closely tied to Autodock.

CHARMM is a general and flexible software application for modeling the structure and behavior of molecular systems.

GAMESS is a general ab initio quantum chemistry package.

Gaussian is a connected system of programs for performing semiempirical and ab initio molecular orbital (MO) calculations.

Gromacs is a versatile package to perform molecular dynamics, i.e. simulate the Newtonian equations of motion for systems with hundreds to millions of particles. It is primarily designed for biochemical molecules like proteins and lipids that have a lot of complicated bonded interactions, but since GROMACS is extremely fast at calculating the nonbonded interactions (that usually dominate simulations) many groups are also using it for research on non-biological systems, e.g. polymers.

NAMD is a parallel molecular dynamics program for UNIX platforms designed for high-performance simulations in structural biology. VMD, the associated molecular visualization program, is also available.

A high-performance toolkit for molecular simulation. Use it as an application, a library, or a flexible programming environment. It includes extensive language bindings for Python, C, C++, and even Fortran.

Cross-program parameter and topology file editor and molecular mechanical simulator engine.

Psi4 is an ab-initio electronic structure code that supports various methods for calculating energies and gradients of molecular systems.

Q-Chem is a comprehensive ab initio quantum chemistry package for accurate predictions of molecular structures, reactivities, and vibrational, electronic and NMR spectra.

A limited number of Schrödinger applications are available on the Biowulf cluster through the Molecular Modeling Interest Group. Most are available through the Maestro GUI.

VMD is a molecular visualization program for displaying, animating, and analyzing large biomolecular systems using 3-D graphics and built-in scripting. To use, type vmd at the prompt.

Simulates the behavior of spins under a specific # MR sequence within a microvascular network.

Finite element solver for computational modeling of the cardio vascular system.

A multi-GPU accelerated software package that fully automates motion-corrected marker-free tomographic alignment and reconstruction. It is improved version of AreTomo that allow multi-gpu parallelization and better tilt-series alignment.

AreTomo3 is a multi-GPU accelerated software package that enables real-time fully automated reconstruction of cryoET tomograms in parallel with cryoET data collection. Integrating MotionCor3, AreTomo2, and GCtfFind in a single application, AreTomo3 has established an autonomous preprocessing pipeline that, whenever a new tilt series is collected, can activate and repeat itself from correction of beam induced motion and assembling tilt series to CTF estimation and correction, tomographic alignment, and 3D reconstruction throughout entire session of data collection without human intervention.

The Collaborative Computational Project for electron cryo-microscopy (CCP-EM) supports users and developers in biological cryogenic EM.

cisTEM is user-friendly software to process cryo-EM images of macromolecular complexes and obtain high-resolution 3D reconstructions from them.

Content-Aware Image Restoration for Cryo-Transmission Electron Microscopy Data: utilizes advanced image processing algorithms to enhance the quality of images obtained from Cryo-TEM (Cryo-Transmission Electron Microscopy) by selectively improving the contrast and detail of the important features within the image while minimizing noise, effectively making the data more interpretable for analysis, particularly in structural biology applications.

CryoDRGN is an algorithm that leverages the representation power of deep neural networks to directly reconstruct continuous distributions of 3D density maps and map per-particle heterogeneity of single-particle cryo-EM datasets. It contains interactive tools to visualize a dataset’s distribution of per-particle variability, generate density maps for exploratory analysis, extract particle subsets for use with other tools and generate trajectories to visualize molecular motions.

CryoPROS is a computational framework specifically designed to tackle misalignment errors caused by preferred orientation issues in single-particle cryo-EM.

CryoSieve is an advanced software solution designed for particle sorting/sieving in single particle analysis (SPA) for Cryogenic Electron Microscopy (cryo-EM).

CryoSPARC (Cryo-EM Single Particle Ab-Initio Reconstruction and Classification) is a state of the art HPC software solution for complete processing of single-particle cryo-electron microscopy (cryo-EM) data. CryoSPARC is useful for solving cryo-EM structures of membrane proteins, viruses, complexes, flexible molecules, small particles, phase plate data and negative stain data.

Programs for finding CTFs of electron micrographs

DeepDeWedge is A Deep Learning Method for Simultaneous Denoising and Missing Wedge Reconstruction in Cryogenic Electron Tomography".

A deep learning solution for cryo-EM map.

EMAN2 is a broadly based greyscale scientific image processing suite with a primary focus on processing data from transmission electron microscopes.

An application improving the quality and interpretability of cryo-EM maps by local and non-local deep learning.

Frealign is a program for high-resolution refinement of 3D reconstructions from cryo-EM images of single particles.

Frealix is a program for the refinement of helical filament reconstructions from cryo electron micrographs. It is primarily used to process images of amyloid fibrils, though it has also been tested on TMV and actin filaments.

Fully automatic acccurate, convenient and extremely fast particle picking for EM

Gctf provides accurate estimation of the contrast transfer function (CTF) for near-atomic resolution cryo electron microscopy (cryoEM) reconstruction using GPUs.

ModelAngelo is an automatic atomic model building program for cryo-EM maps.

MotionCor2 is a multi-GPU accelerated program that provides iterative, patch-based motion detection combining spatial and temporal constraints and dose weighting for both single particle and tomographic cryo-electon microscopy images.

An improved implementation of MotionCor2 with addition of CTF estimation, Multi-GPU accelerated software package that enables single-pixel level correction of anisotropic beam induced sample motion for cryoEM and cryET images.

Installed with one user's request for their cryo-ET classification of the nuclear pore complex C++ version of subtomogram averaging (SA) scripts from TOM/AV3 package https://doi.org/10.1073/pnas.0409178102. Both CPU and GPU parallelization is supported although the latter performs significantly worse, build against CUDA10. Recommend the user to only cpus.

Protofilament Refinement is a software package to further refine microtubule structures, by aligning individual protofilaments rather than full microtubule segements.

UCSF pyem is a collection of Python modules and command-line utilities for electron microscopy of biological samples.

RELION (for REgularised LIkelihood OptimisatioN) is a stand-alone computer program for Maximum A Posteriori refinement of (multiple) 3D reconstructions or 2D class averages in cryo-electron microscopy.

ResMap (Resolution Map) is a Python (NumPy/SciPy) application with a Tkinter GUI and a command-line interface. It is a software package for computing the local resolution of 3D density maps studied in structural biology, primarily electron cryo-microscopy (cryo-EM).

Surface_morphometrics: A pipeline of tools to generate robust open mesh surfaces from voxel segmentations of biological membranes using the Screened Poisson algorithm, calculate morphological features including curvature and membrane-membrane distance using pycurv's vector voting framework, and tools to convert these orphological quantities into morphometric insights.

TomoSegMemTV is a software package for segmenting membranes in tomograms. It is based on a Gaussian-like model of membrane profile, a local differential structure approach and anisotropic propagation of the local structural information using the tensor voting algorithm. Standalone Linux version.

topaz is a pipeline for particle detection in cryo-electron microscopy images using convolutional neural networks trained from positive and unlabeled examples. Topaz also includes methods for micrograph and tomogram denoising using deep denoising models.

Warp is a set of tools for cryo-EM and cryo-ET data processing including: Warp, M, WarpTools, MTools, MCore, and Noise2Map.

B-SOID (Behavioral Segmentation in Deeplabcut) is an unsupervised learning algorithm that serves to discover and classify behaviors that are not pre-defined by users. It segregates statistically different, sub-second rodent behaviors with a single bottom-up perspective video cameraR by performing a novel expectation maximization fitting of Gaussian mixture models on t-Distributed Stochastic Neighbor Embedding (t-SNE).

BioBERT is a biomedical language representation model designed for biomedical text mining tasks such as biomedical named entity recognition, relation extraction, question answering, etc.

BioGANs is a novel application of Generative Adversarial Networks (GAN) to the synthesis of cells imaged by fluorescence microscopy. It allows to infer the correlation between the spatial pattern of different fluorescent proteins that reflects important biological functions. The synthesized images capture these relationships, which are relevant for biological applications.

The NVIDIA CUDA Deep Neural Network library (cuDNN) is a GPU-accelerated library of primitives for deep neural networks. cuDNN provides highly tuned implementations for standard routines such as forward and backward convolution, pooling, normalization, and activation layers.

DanQ is a hybrid convolutional and recurrent deep neural network for quantifying the function of DNA sequences

Real-time denoising of fluorescence time-lapse imaging using deep self-supervised learning

The DeepCell-tf library allows users to apply pre-existing models to imaging data as well as to develop new deep learning models for single-cell analysis. The library specializes in models for cell segmentation (whole-cell and nuclear) in 2D and 3D images as well as cell tracking in 2D time-lapse datasets. The models are applicable to data ranging from multiplexed images of tissues to dynamic live-cell imaging movies.

DeepLabCut is an open source toolbox that builds on a state-of-the-art human pose estimation algorithm. It allows training of a deep neural network by using limited training data to precisely track user-defined features, so that the human labeling accuracy will be matched.

This project aims to offer easy access to Deep Learning for segmentation of structures of interest in biomedical 3D scans. It is a system that allows the easy creation of a 3D Convolutional Neural Network, which can be trained to detect and segment structures if corresponding ground truth labels are provided for training. The system processes NIFTI images, making its use straightforward for many biomedical tasks.

DeepMM implements fully automated de novo structure modeling method, MAINMAST, which builds three-dimensional models of a protein from a near-atomic resolution EM map. The method directly traces the protein’s main-chain and identifies Cα positions as tree-graph structures in the EM map.

Deep-Learning framework for predicting structural variations from bulk and single-cell contact maps.

The GCN_Cancer application employs graph convolutional network (GCN) models to classify the gene expression data samples from The Cancer Genonme Atlas (TCAG) as 33 designated tumor types or as normal. It has been trained on 10,340 cancer samples and 731 normal tissue samples from TCGA dataset.

RoseTTAFold All-Atom (RFAA), a deep network capable of modeling full biological assemblies containing proteins, nucleic acids, small molecules, metals, and covalent modifications given the sequences of the polymers and the atomic bonded geometry of the small molecules and covalent modifications. The heme_binder_diffusion pipeline employs RoseTTAFold All-Atom to perform de novo heme binding protein design.

IsoNet is a deep learning-based software package that iteratively reconstructs the missing-wedge information and increases signal-to-noise ratio, using the knowledge learned from raw tomograms. Without the need for sub-tomogram averaging, IsoNet generates tomograms with significantly reduced resolution anisotropy.

Keras is a deep learning API written in Python, running on top of the machine learning platform TensorFlow. It was developed with a focus on enabling fast experimentation. Being able to go from idea to result as fast as possible is key to doing good research.

Kilosort is a state-of-the-art neural signal processing algorithm designed for the analysis of large-scale neural recordings. It excels in sorting spikes in high-density multielectrode array recordings, making it an essential tool for neuroscientists.

The MPSAN_MP application employs neural networks implementing the Message Passing or Self-Attention mechanism for Molecular Property prediction. It takes as input a set of molecules represented by either graphs or SMILES strings. For each molecule, it predicts a propery value (PV), which may be either a discrete/binary label (classification task) or a continuous value (regression task). of molecular property prediction:
This application is been be used as a biological example by the class #7 of the course "Deep Learning by Example on Biowulf".

NeST-VNN is an interpretable neural network-based model that predicts cell response to a drug. This framework integrates information across multiple levels of cancer cell biology to understand drug response, and can serve to identify and explain biomarkers for clinical application.

nnU-Net is a PyTorch-based semantic segmentation application. It will analyze the provided training cases and automatically configure a matching U-Net-based segmentation pipeline.

PyTom is a software package for the analysis of volumetric data obtained by cryo electron tomography (cryo-ET). It covers a complete pipeline of processing steps for tomogram reconstruction, localization of macromolecular complexes in tomograms, fine alignment of subtomograms extracted at these locations, and their classification.

ReLeaSE (Reinforcement Learning for Structural Evolution) is an application for de-novo Drug Design based on Reinforcement Learning. It integrates two deep neural networks: generative and predictive, that are trained separately but are used jointly to generate novel targeted chemical libraries. ReLeaSE uses simple representation of molecules by their simplified molecular input line entry specification (SMILES) strings only.

Rosetta Fold (RF) dissusion is an open source method for structure generation, with or without conditional information (a motif, target etc). It can perform motif scaffolding, unconditional protein generation, and other tasks.

Saturn is a deep learning approach that couples gene expression with protein representations learnt using large protein language models for cross-species integration. The key idea in SATURN is to map cells from all datasets to a shared space of functionally related genes that we name macrogenes. Using macrogenes, SATURN is uniquely able to detect functionally related genes co-expressed across species.

SpliceAI is a deep neural network that accurately predicts splice junctions from an arbitrary pre-mRNA transcript sequence, enabling precise prediction of noncoding genetic variants that cause cryptic splicing.

Talos is a hyperparameter optimization package for deep learning. It works with any Keras, TensorFlow (tf.keras) or PyTorch model, takes minutes to implement, involves no new syntax to learn and adds zero new overhead to your workflow.

ensoorQTL leverages general-purpose libraries and graphics processing units (GPUs) to achieve high efficiency of computations at low costR. Using PyTorch or TensorFlow it allows > 200-fold decreases in runtime and ~ 5–10-fold reductions in cost when running on GPUs relative to CPUs.

Tybalt implements a Variational EutoEncoder (VAE), a deep neural network approach capable of generating meaningful latent spaces for image and text data. Tybalt has been trained on The Cancer Genome Atlas (TCGA) pan-cancer RNA-seq data and used to identify specific patterns in the VAE encoded features.

U-Net is an image segmentation tool. It relies on the strong use of data augmentation to use the available annotated samples more efficiently. The architecture consists of a contracting path to capture context and a symmetric expanding path that enables precise localization.

Emacs is a text and source code editor for text terminals and X. It has a vast set of features and is well suited for doing everything from reading mail and simple text editing to managing and editing large programming projects. It has its own help and tutorial which can be accessed by typing Ctrl-h i and Ctrl-h t respectively. Type emacs [filename] to edit a file. For more info, see here.

Emacs Speaks Statistics is an Emacs mode for interactive statistical programming and data analysis. Languages supported: the S family (S, S-PLUS and R), SAS, BUGS/JAGS, Stata and XLispStat. First load an R module per our R applications page. Putting the line (load "/usr/local/share/emacs/site-lisp/ess-17.11/lisp/ess-site") in your .emacs file, or its one-time equivalent M-x load-library /usr/local/share/emacs/site-lisp/ess-17.11/lisp/ess-site, will make an *ESS* buffer available.

Nano is a simple, user-friendly text editor derived from the editor in the Pine email client. Type nano [filename] to edit a file and use the key commands listed at the bottom of the screen to access various functions. It is equivalent to Pico.

NEdit is an GUI style editor for plain text and source code files. It provides mouse based editing and a streamlined editing style, based on popular Macintosh and MS Windows editors, using the X-window system. NEdit requires an X-based workstation or X-Terminal. Type nedit [filename] to edit a file.

Neovim is a refactor, and sometimes redactor, in the tradition of Vim (which itself derives from Stevie). It is not a rewrite but a continuation and extension of Vim.

SciTE or SCIntilla based Text Editor is a cross-platform text editor. Lightweight and built for speed, it is designed mainly for source editing, and performs syntax highlighting and inline function reference for many different languages.

Tex + Latex + associated packages for high-quality text formatting.

Texinfo is the official documentation format of the GNU project

a text editor that is upwards compatible to Vi. It can be used to edit all kinds of plain text. It is especially useful for editing programs with syntactical coloring.

Visual Studio Code is a lightweight but powerful source code editor which runs on your desktop and is available for Windows, macOS and Linux. It comes with built-in support for JavaScript, TypeScript and Node.js and has a rich ecosystem of extensions for other languages (such as C++, C#, Java, Python, PHP, Go) and runtimes (such as .NET and Unity).

Abyss represents Assembly By Short Sequences - a de novo, parallel, paired-end sequence assembler. The parallel version is implemented using MPI and is capable of assembling larger genomes.

Automatic Filtering, Trimming, Error Removing and Quality Control for fastq data.

ANGSD is a software for analyzing next generation sequencing data. The software can handle a number of different input types from mapped reads to imputed genotype probabilities. Most methods take genotype uncertainty into account instead of basing the analysis on called genotypes. This is especially useful for low and medium depth data. The software is written in C++ and has been used on large sample sizes.

AnnotSV is a program designed for annotating Structural Variations (SV). This tool compiles functionally, regulatory and clinically relevant information and aims at providing annotations useful to i) interpret SV potential pathogenicity and ii) filter out SV potential false positives.

ANNOVAR is an efficient software tool to utilize update-to-date information to functionally annotate genetic variants detected from diverse genomes.

AscatNGS contains the Cancer Genome Projects workflow implementation of the ASCAT copy number algorithm for paired end sequencing.

AUGUSTUS is a program that predicts genes in eukaryotic genomic sequences.

This tool is used to extract raw sequences (with qualities) from bam files.

bamliquidator is a set of tools for analyzing the density of short DNA sequence read alignments in the BAM file format.

Bam-readcount generates metrics at single nucleotide positions. There are number of metrics generated which can be useful for filtering out false positive calls.

BamTools provides a fast, flexible C++ API & toolkit for reading, writing, and manipulating BAM files.

bamUtil is a repository that contains several programs that perform operations on SAM/BAM files. All of these programs are built into a single executable, bam.

Bartender is an accurate clustering algorithm to detect barcodes and their abundances from raw next-generation sequencing data. In contrast with existing methods that cluster based on sequence similarity alone, Bartender uses a modified two-sample proportion test that also considers cluster size. This modification results in higher accuracy and lower rates of under- and over-clustering artifacts.

Command line interface for Illumina's BaseSpace

A tool to extract paired reads in FASTQ format from coordinate sorted BAM files

An extensive set of bioinformatics tools including bbmap (short read aligner), bbnorm (kmer based normalization), dedupe (deduplication and clustering of unaligned reads), reformat (formatting and trimming reads) and many more.

The Illumina BCL Convert is a standalone local software app that converts the Binary Base Call (BCL) files produced by Illumina sequencing systems to FASTQ files. BCL Convert also provides adapter handling (through masking and trimming) and UMI trimming and produces metric outputs.

a tool to handle bcl conversion and demultiplexing

Bedops is a suite of tools to address common questions raised in genomic studies - mostly with regard to overlap and proximity relationships between data sets - BEDOPS aims to be scalable and flexible, facilitating the efficient and accurate analysis and management of large-scale genomic data.

The BEDTools utilities allow one to address common genomics tasks such finding feature overlaps and computing coverage. In addition, one can develop sophisticated pipelines that answer complicated research questions by "streaming" several BEDTools together.

Binding and expression target analysis (BETA) is a software package that integrates ChIP-seq of TFs or chromatin regulators with differential gene expression data to infer direct target genes. The combination of ChIP-seq and transcriptome analysis is a compelling approach to unravel the regulation of gene expression.

bigSCale is a complete framework for the analysis and visualization of single cell data. It allows to cluster, phenotype, perform pseudotime analysis, infer gene regulatory networks and reduce large datasets in smaller datasets with higher quality.

Regular awk with support for several common biological data formats, including optionally gzip'ed BED, GFF, SAM, VCF, FASTA/Q and TAB-delimited formats with column names.

bioBakery is a meta’omic analysis environment and collection of individual software tools with the capacity to process raw shotgun sequencing data into actionable microbial community feature profiles, summary reports, and publication-ready figures. It includes a collection of preconfigured analysis modules also joined into workflows for reproducibility. Each individual module has been developed to perform a particular task, e.g. quantitative taxonomic profiling or statistical analysis.

Tools for early stage alignment file processing.

tool (and library) to manipulate Biological Observation Matrix (BIOM) Format files

biomodal duet multiomics solution bioinformatics pipeline

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.

A PyTorch Basecaller for Oxford Nanopore Reads

bowtie is an ultrafast, memory-efficient short read aligner geared toward quickly aligning large sets of short DNA sequences (reads) to large genomes.

A version of bowtie that's particularly good at aligning reads of about 50 up to 100s or 1,000s of characters, and particularly good at aligning to relatively long (e.g. mammalian) genomes

BRASS analyses one or more related BAM files of paired-end sequencing to determine potential rearrangement breakpoints.

provides genome-wide detection of structural variants from next generation paired-end sequencing reads.

breseq is a computational pipeline for finding mutations relative to a reference sequence in short-read DNA re-sequencing data. It is intended for haploid microbial genomes (<20 Mb).

BSMAP is a short reads mapping software for bisulfite sequencing reads. Bisulfite treatment converts unmethylated Cytosines into Uracils (sequenced as Thymine) and leave methylated Cytosines unchanged, hence provides a way to study DNA cytosine methylation at single nucleotide resolution. BSMAP aligns the Ts in the reads to both Cs and Ts in the reference.

BUSCO completeness assessments employ sets of Benchmarking Universal Single-Copy Orthologs from OrthoDB (www.orthodb.org) to provide quantitative measures of the completeness of genome assemblies, annotated gene sets, and transcriptomes in terms of expected gene content.

BWA is a fast light-weighted tool that aligns short sequences to a sequence database, such as the human reference genome.

The next version of the bwa-mem algorithm in bwa.

CADD (Combined Annotation Dependent Depletion) is a tool for scoring the deleteriousness of single nucleotide variants as well as insertion/deletions variants in the human genome. Currently, it supports the builds: GRCh37/hg19 and GRCh38/hg38.

Canu is a fork of the Celera Assembler designed for high-noise single-molecule sequencing (such as the PacBio RSII or Oxford Nanopore MinION). Canu will correct the reads, then trim suspicious regions (such as remaining SMRTbell adapter), then assemble the corrected and cleaned reads into unitigs.

Canvas is a tool for calling copy number variants (CNVs) from human DNA sequencing data.

CCBR Pipeliner provides access to a set of best-practices NGS pipelines developed, tested, and benchmarked by experts at CCBR and NCBR for Biowulf. Contact CCBR_Pipeliner@mail.nih.gov with questions

Cis-regulatory Element Annotation System is a tool designed to characterize genome-wide protein-DNA interaction patterns from ChIP-chip and ChIP-Seq of both sharp and broad binding factors.

CellBender is a software package for eliminating technical artifacts from high-throughput single-cell omics data, including scRNA-seq, snRNA-seq, and CITE-seq.

An open-source application for biological image analysis

Cell Ranger is a set of analysis pipelines that processes Chromium single cell 3’ RNA-seq output to align reads, generate gene-cell matrices and perform clustering and gene expression analysis.

Cell Ranger ARC is a set of analysis pipelines that process Chromium Single Cell Multiome ATAC + Gene Expression sequencing data to generate a variety of analyses pertaining to gene expression, chromatin accessibility and their linkage. Furthermore, since the ATAC and gene expression measurements are on the very same cell, we are able to perform analyses that link chromatin accessibility and gene expression.

Cell Ranger ATAC is a set of analysis pipelines that process Chromium Single Cell ATAC data.

Efficient genotyping bi-allelic SNPs on single cells

Detect subclonality and copy number in matched NGS data

Rapid assessment of genome bin quality using machine learning.

The CHESS (Comparison of Hi-C Experiments using Structural Similarity) application implements an algorithm for the comparison of chromatin contact maps and automatic differential feature extraction.

Filtering and trimming for long-read sequencing data (PacBio/ONT).

ChromHMM is software for learning and characterizing chromatin states.

A combined strategy to identify circular RNAs (circRNAs and ciRNAs)

Circos is a program for the generation of publication-quality, circularly composited renditions of genomic data and related annotations. Circos is particularly suited for visualizing alignments, conservation and intra and inter-chromosomal relationships. Also, Circos is useful to visualize any type of information that benefits from a circular layout. Thus, although it has been designed for the field of genomics, it is sufficiently flexible to be used in other data domains.

Clair3 is a small variant caller for Illumina, PacBio and ONT long reads. Compare to PEPPER (r0.4), Clair3 (v0.1) shows a better SNP F1-score with ≤30-fold of ONT data (precisionFDA Truth Challenge V2), and a better Indel F1-score, while runs generally four times faster.

A method based on a supervised sequence classification using discriminative k-mers

Robust detection of clinically relevant structural and copy number variation from whole genome sequencing data

Copy number variant detection from targeted DNA sequencing

CNVnator is a tool for CNV discovery and genotyping from depth of read mapping.

Cogent NGS Analysis Pipeline (CogentAP) is bioinformatic software for analyzing RNA-seq NGS data generated using various takara kits.

A library to combine, analyze, group and correct p-values in BED files. Unique tools involve correction for spatial autocorrelation. This is useful for ChIP-Seq probes and Tiling arrays, or any data with spatial correlation.

Concordance and contamination estimator for tumor–normal pairs

Software pipeline for the analysis of CRISPR-Cas9 genome editing outcomes from sequencing data

CrossMap is a program for convenient conversion of genome coordinates between different assemblies (e.g. mm9->mm10). It can convert SAM, BAM, bed, GTF, GFF, wig/bigWig, and VCF files

csvkit is a suite of command-line tools for converting to and working with CSV, the king of tabular file formats.

Cufflinks assembles transcripts and estimates their abundances in RNA-Seq samples. It accepts aligned RNA-Seq reads and assembles the alignments into a parsimonious set of transcripts. Cufflinks then estimates the relative abundances of these transcripts based on how many reads support each one.

cutadapt removes adapter sequences from DNA high-throughput sequencing data. This is usually necessary when the read length of the machine is longer than the molecule that is sequenced, such as in microRNA data.

cutruntools2 is a major update of CutRunTools, including a set of new features specially designed for CUT&RUN and CUT&Tag experiments. Both of the bulk and single-cell data can be processed, analyzed and interpreted.

A toolkit for Dynamic Analysis of Nucleosome and Protein Occupancy by Sequencing, version 2

The DeepCell-tf library allows users to apply pre-existing models to imaging data as well as to develop new deep learning models for single-cell analysis. The library specializes in models for cell segmentation (whole-cell and nuclear) in 2D and 3D images as well as cell tracking in 2D time-lapse datasets. The models are applicable to data ranging from multiplexed images of tissues to dynamic live-cell imaging movies.

DeepConsensus uses gap-aware sequence transformers to correct errors in Pacific Biosciences (PacBio) Circular Consensus Sequencing (CCS) data

A deep-learning method for detecting DNA methylation state from Oxford Nanopore sequencing reads.

deepTools is a suite of user-friendly tools for the visualization, quality control and normalization of data from deep-sequencing DNA sequencing experiments.

DeepVariant is an analysis pipeline that uses a deep neural network to call genetic variants from next-generation DNA sequencing data.

deFuse is a software package for gene fusion discovery using RNA-Seq data. The software uses clusters of discordant paired end alignments to inform a split read alignment analysis for finding fusion boundaries. The software also employs a number of heuristic filters in an attempt to reduce the number of false positives and produces a fully annotated output for each predicted fusion.

DELLY is an integrated structural variant prediction method that can detect deletions, tandem duplications, inversions and translocations at single-nucleotide resolution in short-read massively parallel sequencing data. It uses paired-ends and split-reads to sensitively and accurately delineate genomic rearrangements throughout the genome.

DETONATE is a tool for evaluation of de novo transcriptome assemblies from RNA-Seq data. It consists of two component packages, RSEM-EVAL and REF-EVAL. RSEM-EVAL is a reference-free evaluation method based on a novel probabilistic model that depends only on an assembly and the RNA-Seq reads used for its construction. REF-EVAL is a toolkit of reference-based measures.

Genome-wide quantification of differential transcription factor activity

A modular Hi-C mapping pipeline for reproducible data analysis, it was used for Micro-C analysis too.

Dorado is a high-performance, easy-to-use, open source basecaller for Oxford Nanopore reads.

A pipeline to find aberrant events in RNA-Seq data.

Command-line tools for processing biological sequencing data. Barcode demultiplexing, adapter trimming, etc.

EDD is a ChIP-seq peak caller for detection of megabase domains of enrichment.

The EDTA package was designed to filter out false discoveries in raw TE candidates and generate a high-quality non-redundant TE library for whole-genome TE annotations.

This pipeline is designed for automated end-to-end quality control and processing of ATAC-seq or DNase-seq data.

EPACTS (Efficient and Parallelizable Association Container Toolbox) is a versatile software pipeline to perform various statistical tests for identifying genome-wide association from sequence data through a user-friendly interface, both to scientific analysts and to method developers.

ESPRESSO is a novel method for processing alignment of long read RNA-seq data, which can effectively improve splice junction accuracy and isoform quantification. ESPRESSO jointly considers alignments of all long reads aligned to a gene and uses error profiles of individual reads to improve the identification of splice junctions and the discovery of their corresponding transcript isoforms.

EukDetect: Detect eukaryotes from shotgun metagenomic data

The Exomiser is a Java program that functionally annotates variants from whole-exome sequencing data starting from a VCF file.

eXpress is a streaming tool for quantifying the abundances of a set of target sequences from sampled subsequences.

Fragmentase Artifact Detection and Elimination

FAN-C is a toolkit for the analysis and visualization of Hi-C data. Beyond objects generated within FAN-C, the toolkit is largely compatible with Hi-C files from Cooler and Juicer.

It provide quality control functions to next gen sequencing data.

fastq-tools a collection of small and efficient programs for performing some common and uncommon tasks with FASTQ files.

FastQ Screen allows you to screen a library of sequences in FastQ format against a set of sequence databases so you can see if the composition of the library matches with what you expect.

The FASTX-Toolkit is a collection of command line tools for Short-Reads FASTA/FASTQ files preprocessing.

FCS is a toolset to remove contaminant sequences from a genome assembly.

The Fulcrum Genomics tools are a set of utilities for working with BAM files, VCF files, and Unique Molecular IDs. Theey are accessed as subprograms from a Java jar, like GATK or Picard.

Fit-Hi-C is a tool for assigning statistical confidence estimates to intra-chromosomal contact maps produced by genome-wide genome architecture assays such as Hi-C.

Gene-flank analysis tool

FlashPCA performs fast principal component analysis (PCA) of single nucleotide polymorphism (SNP) data, similar to smartpca from EIGENSOFT (http://www.hsph.harvard.edu/alkes-price/software/) and shellfish (https://github.com/dandavison/shellfish). FlashPCA is based on the https://github.com/yixuan/spectra/ library.

Flexbar preprocesses high-throughput sequencing data efficiently. It demultiplexes barcoded runs and removes adapter sequences. Moreover, trimming and filtering features are provided. Flexbar increases read mapping rates and improves genome and transcriptome assemblies. It supports next-generation sequencing data in fasta and fastq format, e.g. from Illumina and the Roche 454 platform

Fast and accurate de novo assembler for single molecule sequencing reads

Tools for manipulating fastq files

Bayesian haplotype-based polymorphism discovery and genotyping

Control-FREEC is a tool for detection of copy-number changes and allelic imbalances (including LOH) using deep-sequencing data

Fusion Gene Detection Using Whole-Exome Sequencing Data in Cancer Patients

FusionCatcher searches for novel/known somatic fusion genes, translocations, and chimeras in RNA-seq data (paired-end or single-end reads from Illumina NGS platforms like Solexa/HiSeq/NextSeq/MiSeq) from diseased samples.

GATK, from the Broad Institute, is a structured software library that makes writing efficient analysis tools using next-generation sequencing data very easy, and second it's a suite of tools for working with human medical resequencing projects such as 1000 Genomes and The Cancer Genome Atlas. These tools include things like a depth of coverage analyzers, a quality score recalibrator, a SNP/indel caller and a local realigner.

High resolution peak calling and motif discovery for ChIP-seq and ChIP-exo data

GEMINI (GEnome MINIng) is designed to be a flexible framework for exploring genetic variation in the context of the wealth of genome annotations available for the human genome. By placing genetic variants, sample genotypes, and useful genome annotations into an integrated database framework, GEMINI provides a simple, flexible, yet very powerful system for exploring genetic variation for for disease and population genetics.

GEMMA is the software implementing the Genome-wide Efficient Mixed Model Association algorithm for a standard linear mixed model and some of its close relatives for genome-wide association studies (GWAS).

geNomad's primary goal is to identify viruses and plasmids in sequencing data (isolates, metagenomes, and metatranscriptomes).

The Genome Browser Mirror Fragments is a mirror of the UCSC Genome Browser. The URL is https://hpcnihapps.cit.nih.gov/genome. Users can also access the MySQL databases, supporting files directly, and a huge number of associated executables.

Genrich is a peak-caller for genomic enrichment assays (e.g. ChIP-seq, ATAC-seq). It analyzes alignment files generated following the assay and produces a file detailing peaks of significant enrichment.

The GeoMx NGS Pipeline, processes RNA-sequencing files (FASTQ files) from Illumina sequencers according to parameters defined in the Configuration File (which is generated from the GeoMx DSP run). The Pipeline processes information from these files and outputs .dcc files, which can then be uploaded to the GeoMx DSP system for data analysis.

Facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers.

A Genomic Mapping and Alignment Programs

GoPeaks is a peak caller designed for CUT&TAG/CUT&RUN sequencing data. GoPeaks by default works best with narrow peaks such as H3K4me3 and transcription factors. However, broad epigenetic marks like H3K27Ac/H3K4me1 require different the step, slide, and minwidth parameters.

Seed-and-extend program for aligning long error-prone reads to genome graphs.

GRIDSS is a module software suite containing tools useful for the detection of genomic rearrangements. GRIDSS includes a genome-wide break-end assembler, as well as a structural variation caller for Illumina sequencing data. GRIDSS calls variants based on alignment-guided positional de Bruijn graph genome-wide break-end assembly, split read, and read pair evidence.

GT-Pro utilizes an exact matching algorithm to perform ultra-rapid and accurate genotyping of known SNPs from metagenomes.

This module makes available the tools used in the GTEX RNA-Seq pipeline.

GTOOL is a program for transforming sets of genotype data for use with the programs SNPTEST and IMPUTE. GTOOL can be used to (a) generate subsets of genotype data, (b) to convert genotype data between the PED file format and the FILE FORMAT used by SNPTEST and IMPUTE.

Genome UNClutterer (GUNC) is a tool for detection of chimerism and contamination in prokaryotic genomes resulting from mis-binning of genomic contigs from unrelated lineages

Local accelerated basecalling for Nanopore data

Hail is an open-source, scalable framework for exploring and analyzing genomic data.

A set of programs based on htslib to benchmark variant calls against gold standard truth datasets.

Tools to process, normalize and visualize Hi-C data

hichipper is a preprocessing and QC pipeline for HiChIP data. This package takes output from a HiC-Pro run and a sample manifest file (.yaml) that coordinates optional high-quality peaks (identified through ChIP-Seq) and restriction fragment locations (see folder here) as input and produces output that can be used to 1) determine library quality, 2) identify and characterize DNA loops and 3) interactively visualize loops.

HiC-Pro: An optimized and flexible pipeline for Hi-C data processing

A framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines.

Pacific Biosciences workflow for analyzing human PacBio whole genome sequencing data.

Hifiasm is a fast haplotype-resolved de novo assembler initially designed for PacBio HiFi reads. Its latest release supports telomere-to-telomere assembly by utilizing ultralong Oxford Nanopore reads. It can produce better haplotype-resolved assemblies when given parental short reads or Hi-C data.

a computational method to detect CNVs and Translocations from Hi-C data.

Tool for genotyping short tandom repeats from Illumina sequencing data

HISAT is a fast and sensitive spliced alignment program which uses Hierarchical Indexing for Spliced Alignment of Transcripts.

Fast HLA type inference from whole-genome data. Previously known as HLA-PRG-LA.

HMMRATAC peak caller for ATAC-seq data

HOMER (Hypergeometric Optimization of Motif EnRichment) is a suite of tools for Motif Discovery and ChIP-Seq analysis.

High-Throughput Genome-Wide Translocation Sequencing pipeline

HTSbox is a fork of early HTSlib. It is a collection of small experimental tools manipulating HTS-related files. While some of these tools are already part of the official SAMtools package, others are for niche use cases.

HTSeq is a Python package that provides infrastructure to process data from high-throughput sequencing assays.

HUMAnN is a pipeline for efficiently and accurately profiling the presence/absence and abundance of microbial pathways in a community from metagenomic or metatranscriptomic sequencing data (typically millions of short DNA/RNA reads).

IDBA is a practical iterative De Bruijn Graph De Novo Assembler for sequence assembly in bioinfomatics. Most assemblers based on de Bruijn graph build a de Bruijn graph with a specific k to perform the assembling task. For all of them, it is very crucial to find a specific value of k. If k is too large, there will be a lot of gap problems in the graph. If k is too small, there will a lot of branch problems. IDBA uses not only one specific k but a range of k values to build the iterative de Bruijn graph. It can keep all the information in graphs with different k values.

iDEP (integrated Differential Expression and Pathway analysis) is shiny application for analyzing RNA-seq data

iDiffIR is a tool for identifying differential IR from RNA-seq data. It accepts any sorted, indexed BAM file for single- or paired-end reads.

The IDR (Irreproducible Discovery Rate) framework is a unified approach to measure the reproducibility of findings identified from replicate experiments and provide highly stable thresholds based on reproducibility. The IDR method compares a pair of ranked lists of identifications (such as ChIP-seq peaks).

The Integrative Genomics Viewer is a high-performance visualization tool for interactive exploration of large, integrated genomic datasets.

IGVTools provides utilities for working with ascii file formats used by the Integrated Genome Viewer. The files can be sorted, tiled, indexed, and counted.

The Immcantation framework provide a start-to-finish analytical ecosystem for high-throughput AIRR-seq datasets. Immcantation includes Alakazam, Change-O, Dowser, pRESTO, SCOPer, SHazaM, and TIgGER.

Impute is a program for estimating ("imputing") unobserved genotypes in SNP association studies.

In 2015, the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) published updated standards and guidelines for the clinical interpretation of sequence variants with respect to human diseases on the basis of 28 criteria. However, variability between individual interpreters can be extensive because of reasons such as the different understandings of these guidelines and the lack of standard algorithms for implementing them, yet computational tools for semi-automated variant interpretation are not available. To address these problems, InterVar implements these criteria to help human reviewers interpret the clinical significance of variants. InterVar can take a pre-annotated or VCF file as input and generate automated interpretation on 18 criteria.

a tool for intersection and visualization of multiple genomic region sets

iPHoP stands for Integrated Phage Host Prediction. It is an automated command-line pipeline designed to predict the host genus of novel bacteriophages and archaeoviruses based on their genome sequences.

IsoQuant is used to analyze long read RNA sequencing data

IVA is a de novo assembler designed to assemble virus genomes that have no repeat sequences, using Illumina read pairs sequenced from mixed populations at extremely high and variable depth.

iVar is a computational package that contains functions broadly useful for viral amplicon-based sequencing. Additional tools for metagenomic sequencing are actively being incorporated into iVar. While each of these functions can be accomplished using existing tools, iVar contains an intersection of functionality from multiple tools that are required to call iSNVs and consensus sequences from viral sequencing data across multiple replicates.

A One-Click System for Analyzing Loop-Resolution Hi-C Experiments

Java tools for bioinformatics

KAT (K-mer Analysis Toolkit) is a suite of tools that analyse Jellyfish hashes or sequence files (fasta or fastq) using kmer counts.

kb-python is a python package for processing single-cell RNA-sequencing. It wraps the kallisto | bustools single-cell RNA-seq command line tools in order to unify multiple processing workflows.

KMC is a disk-based programm for counting k-mers from (possibly gzipped) FASTQ/FASTA files

KmerGO is a user-friendly tool to identify the group-specific sequences on two groups of high throughput sequencing datasets.

KneadData is a tool designed to perform quality control on metagenomic and metatranscriptomic sequencing data, especially data from microbiome experiments.

Kraken is a system for assigning taxonomic labels to short DNA sequences, usually obtained through metagenomic studies

KrakenUniq is a taxonomic sequence classifier that assigns taxonomic labels to short DNA reads. It combines the fast k-mer-based classification of Kraken (Wood and Salzberg, Genome Biology 2014) with an efficient algorithm for assessing the coverage of unique k-mers found in each species in a dataset. On various test datasets, KrakenUniq gives better recall and precision than other methods and effectively classifies and distinguishes pathogens with low abundance from false positives in infectious disease samples. By using the probabilistic cardinality estimator HyperLogLog, KrakenUniq runs as fast as Kraken and requires little additional memory.

LEfSe (Linear discriminant analysis Effect Size) determines the features (organisms, clades, operational taxonomic units, genes, or functions) most likely to explain differences between classes by coupling standard tests for statistical significance with additional tests encoding biological consistency and effect relevance.

The La Jolla Assembler (LJA) is a tool for assemblies of long and accurate reads. It reduces the error rate in these reads by three orders of magnitude (making them nearly error-free) and constructs the de Bruijn graph for large genomes and large k-mer sizes. Since the de Bruijn graph constructed for a fixed k-mer size is typically either too tangled or too fragmented, LJA uses a new concept of a multiplex de Bruijn graph with varying k-mer sizes.

LocusZoom is designed to facilitate viewing of local association results together with useful information about a locus, such as the location and orientation of the genes it includes, linkage disequilibrium coefficients and local estimates of recombination rates

LoFreq is a fast and sensitive variant-caller for inferring SNVs and indels from next-generation sequencing data.

Longshot is a variant calling tool for diploid genomes using long error prone reads such as Pacific Biosciences (PacBio) SMRT and Oxford Nanopore Technologies (ONT). It takes as input an aligned BAM file and outputs a phased VCF file with variants and haplotype information. It can also output haplotype-separated BAM files that can be used for downstream analysis. Currently, it only calls single nucleotide variants (SNVs).

A probabilistic framework for structural variant discovery.

mach2qtl uses dosages/posterior probabilities inferred with MACH as predictors in a linear regression to test association with a quantitative trait

Model-based Analysis of ChIP-Seq (MACS) on short reads sequencers such as Genome Analyzer (Illumina / Solexa). MACS empirically models the length of the sequenced ChIP fragments, which tends to be shorter than sonication or library construction size estimates, and uses it to improve the spatial resolution of predicted binding sites. MACS also uses a dynamic Poisson distribution to effectively capture local biases in the genome sequence, allowing for more sensitive and robust prediction.

Multiple alignment program for amino acid or nucleotide sequences

MAGeCK is Model-based Analysis of Genome-wide CRISPR/Cas9 Knockout (MAGeCK) method for prioritizing single-guide RNAs, genes and pathways in genome-scale CRISPR/Cas9 knockout screens. It demonstrates better performance compared with other methods, identifies both positively and negatively selected genes simultaneously, and reports robust results across different experimental conditions.

MAGeCK-VISPR is a comprehensive quality control, analysis and visualization workflow for CRISPR/Cas9 screens.

Magic-BLAST is a tool for mapping large next-generation RNA or DNA sequencing runs against a whole genome or transcriptome. Each alignment optimizes a composite score, taking into account simultaneously the two reads of a pair, and in case of RNA-seq, locating the candidate introns and adding up the score of all exons. This is very different from other versions of BLAST, where each exon is scored as a separate hit and read-pairing is ignored.

Modeling Alternative Junction Inclusion Quantification. MAJIQ and Voila are two software packages that together define, quantify, and visualize local splicing variations (LSV) from RNA-Seq data.

MAnorm is for quantitative comparison of ChIP-Seq data sets describing transcription factor binding sites and epigenetic modifications. The quantitative binding differences inferred by MAnorm showed strong correlation with both the changes in expression of target genes and the binding of cell type-specific regulators.

Structural variant and indel caller for mapped sequencing data

a set of multiple scripts used to analyze PLAC-Seq and HiChIP data.

mash is a command line tool and library to provide fast genome and metagenome distance estimation using MinHash. Only command line tool is installed

Minimizer based sparse de Bruijn Graph constructor.

medaka is a tool to create a consensus sequence from nanopore sequencing data. This task is performed using neural networks applied from a pileup of individual sequencing reads against a draft assembly.

A draft genome scaffolder that uses multiple reference genomes in a graph-based approach.

MEGAHIT is a single node assembler for large and complex metagenomics NGS reads, such as soil. It makes use of succinct de Bruijn graph (SdBG) to achieve low memory assembly. MEGAHIT can optionally utilize a CUDA-enabled GPU to accelerate its SdBG construction.

Megalodon is a research command line tool to extract high accuracy modified base and sequence variant calls from raw nanopore reads by anchoring the information rich basecalling neural network output to a reference genome/transcriptome.

MEtaGenome ANalyzer that takes a file of reads and a Blast output from comparison against a reference genome, and automatically calculate a taxonomic classification of the reads and if desired, a functional classification.

Improved variant filtering and polishing via k-mer validation

Evaluate genome assemblies with k-mers and more

Meryl: a genomic k-mer counter (and sequence utility) with nice features. It is built into the Celera Assembler and is also available as a stand-alone application. Meryl uses a sorting-based approach that sorts the k-mers in lexicographical order.

MetaBAT: A robust statistical framework for reconstructing genomes from metagenomic data

The METAL software is designed to facilitate meta-analysis of large datasets (such as several whole genome scans) in a convenient, rapid and memory efficient manner.

MetaPhlAn is a computational tool for profiling the composition of microbial communities (Bacteria, Archaea, Eukaryotes and Viruses) from metagenomic shotgun sequencing data (i.e. not 16S) with species-level. With the newly added StrainPhlAn module, it is now possible to perform accurate strain-level microbial profiling.

Ultrafast de novo assembly for long noisy reads (though having no consensus step)

minimac is a low memory, computationally efficient implementation of the MaCH algorithm for genotype imputation. It is designed to work on phased genotypes and can handle very large reference panels with hundreds or thousands of haplotypes. 'mini' refers to the low amount of computational resources it needs.

Minimap2 is a fast sequence mapping and alignment program that can find overlaps between long noisy reads, or map long reads or their assemblies to a reference genome optionally with detailed alignment (i.e. CIGAR).

Miniprot aligns protein sequences to whole genomes with splicing and frameshift. It is primarily intended for annotating protein-coding genes in a new species using known genes from other species.

miRDeep2 is a completely overhauled tool which discovers microRNA genes by analyzing sequenced RNAs.

MISO (Mixture-of-Isoforms) is a probabilistic framework that quantitates the expression level of alternatively spliced genes from RNA-Seq data

mitosuite is a graphical tool for human mitochondrial genome profiling in massively parallel sequencing

MiXCR is a universal software for fast and accurate analysis of T- and B- cell receptor repertoire sequencing data.

MiXeR is Causal Mixture Model for GWAS summary statistics. The version(1.3) contains a Python port of MiXeR, wrapping the C/C++ core. Also data preprocessing code sumstats.py is included too.

A set of tools to manipulate and visualize DNA/RNA base modification data that are stored in bam format.

A bioinformatics tool for working with modified bases from Oxford Nanopore. Specifically for converting modBAM to bedMethyl files using best practices, but also manipulating modBAM files and generating summary statistics.

modPhred is a pipeline for detection, annotation and visualisation of DNA/RNA modifications.

SNV calling from single cell sequencing

a machine learning method that leverages read-based phasing and read-level features to accurately detect mosaic SNVs (SNPs, small indels) from NGS data.It builds on existing algorithms to achieve a multifold increase in specificity.

Fast BAM/CRAM depth calculation for WGS, exome, or targeted sequencing.

A bioinformatics pipeline aimed at the analysis of mitochondrial DNA (mtDNA) in high throughput sequencing studies.

aggregates results for various frequently used bioinformatics tools across multiple samples into a nice visual report

MuSE is an approach to somatic variant calling based on the F81 Markov substitution model for molecular evolution, which models the evolution of the reference allele to the allelic composition of the matched tumor and normal tissue at each genomic locus.

MuTect is a method developed at the Broad Institute for the reliable and accurate identification of somatic point mutations in next generation sequencing data of cancer genomes.

Tools for analyzing and processing long reads and alignments

nanopolish is a software package for signal-level analysis of Oxford Nanopore sequencing data. Nanopolish can calculate an improved consensus sequence for a draft genome assembly, detect base modifications, call SNPs and indels with respect to a reference genome and more (see Nanopolish modules, below).

NeuSomatic is based on deep convolutional neural networks for accurate somatic mutation detection. With properly trained models, it can robustly perform across sequencing platforms, strategies, and conditions. NeuSomatic summarizes and augments sequence alignments in a novel way and incorporates multi-dimensional features to capture variant signals effectively. It is not only a universal but also accurate somatic mutation detection method.

NGMLR is a long-read mapper designed to align PacBio or Oxford Nanopore (standard and ultra-long) to a reference genome with a focus on reads that span structural variations.

ngsplot is an easy-to-use global visualization tool for next-generation sequencing data.

NgsRelate can be used to infer relatedness, inbreeding coefficients and many other summary statistics for pairs of individuals from low coverage Next Generation Sequencing (NGS) data by using genotype likelihoods instead of called genotypes.

Package includes aligner for single-ended and paired-end reads from the Illumina Genome Analyser. Novoalign finds global optimum alignments using full Needleman-Wunsch algorithm with affine gap penalties.

NucleoATAC is a software for nucleosome calling using ATAC-seq. It can identify the rotational and translational positions of nucleosomes with up to base-pair resolution and provide quantitative measures of nucleosome occupancy.

The PacBio obc2fastq software is designed for use with Onso sequencing data. You can demultiplex samples and generate FASTQ files for use with downstream data pipelines through a command-line interface.

Octopus is a mapping-based variant caller that implements several calling models within a unified haplotype-aware framework. Octopus takes inspiration from particle filtering by constructing a tree of haplotypes and dynamically pruning and extending the tree based on haplotype posterior probabilities in a sequential manner. This allows octopus to implicitly consider all possible haplotypes at a given loci in reasonable time.

Tools to manipulate HDF5 files of the Oxford Nanopore .fast5 file format

Pairtools is a simple and fast command-line framework to process sequencing data from a Hi-C experiment. Pairtools perform various operations on Hi-C pairs and occupy the middle position in a typical Hi-C data processing pipeline. Pairtools aim to be an all-in-one tool for processing Hi-C pairs.

the PALEOMIX pipelines are a set of pipelines and tools designed to aid the rapid processing of High-Throughput Sequencing (HTS) data

An updated pipeline for pangenome investigation

The Clara Parabricks toolkit is a set of GPU-accelerated genome analysis tools for secondary analysis of next generation sequencing data.

Web interface designed specifically for the analysis needs of next generation sequencing applications including RNA, small RNA, and DNA sequencing.

Tools for analyzing CpG/5mC data from PacBio HiFi reads aligned to a reference genome

Improved Phased Assembler (IPA) is the official PacBio software for HiFi genome assembly. IPA was designed to utilize the accuracy of PacBio HiFi reads to produce high-quality phased genome assemblies.

The PBSuite contains two projects created for analysis of Pacific Biosciences long-read sequencing data: PBHoney and PBJelly. PBHoney is an implementation of two variant-identification approaches designed to exploit the high mappability of long reads (i.e., greater than 10,000 bp). PBHoney considers both intra-read discordance and soft-clipped tails of long reads to identify structural variants. PBJelly is a highly automated pipeline that aligns long sequencing reads (such as PacBio RS reads or long 454 reads in fasta format) to high-confidence draft assembles. PBJelly fills or reduces as many captured gaps as possible to produce upgraded draft genomes.

A supervised learning framework for chromatin loop detection in genome-wide contact maps.

PEAR is an ultrafast, memory-efficient and highly accurate pair-end read merger. It is fully parallelized and can run with as low as just a few kilobytes of memory.

peddy is used to compare sex and familial relationships given in a PED file with those inferred from a VCF file

PennCNV is a free software tool for Copy Number Variation (CNV) detection from SNP genotyping arrays. Currently it can handle signal intensity data from Illumina and Affymetrix arrays. With appropriate preparation of file format, it can also handle other types of SNP arrays and oligonucleotide arrays.

PEPATAC is a robust pipeline for Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) built on a loosely coupled modular framework. It may be easily applied to ATAC-seq projects of any size, from one-off experiments to large-scale sequencing projects. It is optimized on unique features of ATAC-seq data to be fast and accurate and provides several unique analytical approaches.

PhaseFinder detects DNA inversion mediated phase variation in bacterial genomes using genomic or metagenomic sequencing data. It works by identifying regions flanked by inverted repeats, mimicking their inversion in silico, and identifying regions where sequencing reads support both orientations.

Picard comprises Java-based command-line utilities that manipulate SAM files, and a Java API (SAM-JDK) for creating new programs that read and write SAM files. Both SAM text format and SAM binary (BAM) format are supported.

Pindel can detect breakpoints of large deletions, medium sized insertions, inversions, tandem duplications and other structural variants at single-based resolution from next-gen sequence data. It uses a pattern growth approach to identify the breakpoints of these variants from paired-end short reads.

tool for variant-detection in high-throughput sequencing data.

PLINK is whole genome association analysis toolset, designed to perform a range of basic, large-scale analyses in a computationally efficient manner.

Tool for manipulating the pod5 format of nanopore reads

Pomoxis comprises a set of basic bioinformatic tools tailored to nanopore sequencing. Notably tools are included for generating and analysing draft assemblies. Many of these tools are used by the research data analysis group at Oxford Nanopore Technologies.

Trim/demultiplex Oxford Nanopore reads

Poretools is a toolkit for manipulating and exploring nanopore sequencing data sets. Poretools operates on individual FAST5 files, directory of FAST5 files, and tar archives of FAST5 files.

predicting library complexity and genome coverage in high-throughput sequencing

PRINSEQ is a tool that generates summary statistics of sequence and quality data and that is used to filter, reformat and trim next-generation sequence data. It is particular designed for 454/Roche data, but can also be used for other types of sequence data.

proseq-2.0 is a pipeline for preprocesses and alignment of run-on sequencing (PRO/GRO/ChRO-seq) data from Single-Read or Paired-End Illumina Sequencing

pVACtools is a cancer immunotherapy suite consisting of pVACseq, pVACfuse, pVACvector

Pychopper v2 is a tool to identify, orient and trim full-length Nanopore cDNA reads. The tool is also able to rescue fused reads.

PyClone is statistical model and software tool designed to infer the prevalence of point mutations in heterogeneous cancer samples.

pycoQC is a new tool to generate interactive quality control metrics and plots from basecalled nanopore reads or summary files generated by the basecallers Albacore, Guppy or MinKNOW. pycoQC has several novel features, including: 1) python support for creation of dynamic D3.js visualizations and interactive data exploration in Jupyter Notebooks; 2) simple command line interface to generate customizable interactive HTML reports; and 3) multiprocessing FAST5 feature extraction program to generate a summary file directly from FAST5 files.

A download client for the European Genome-phenome Archive (EGA). The EGA is designed to be a repository for all types of sequence and genotype experiments, including case-control, population, and family studies.

QCTOOL is a command-line utility program for basic quality control of gwas datasets.

A platform-independednt application written in Java and R that provides both a GUI and a co mmand-line interface to facilitate the quality control of alignment sequencing data.

QUAST stands for QUality ASsessment Tool. The tool evaluates genome assemblies by computing various metrics. The package includes the general QUAST tool for genome assemblies, MetaQUAST, the extension for metagenomic datasets, and Icarus, interactive visualizer for these tools.

Ultrafast consensus module for raw de novo genome assembly of long uncorrected reads.

RagTag is a collection of software tools for scaffolding and improving modern genome assemblies.

Rapids-singlecell offers enhanced single-cell data analysis as a near drop-in replacement predominantly for scanpy.Utilizing GPU computing with cupy and Nvidia’s RAPIDS, it emphasizes high computational efficiency.

RAREMETAL is a computationally efficient tool for meta-analysis of rare variants using sequencing or genotyping array data.

Rcorrector implements a k-mer based method to correct random sequencing errors in Illumina RNA-seq reads. Rcorrector uses a De Bruijn graph to compactly represent all trusted k-mers in the input reads. Unlike WGS read correctors, which use a global threshold to determine trusted k-mers, Rcorrector computes a local threshold at every position in a read.

regtools:Tools that integrate DNA-seq and RNA-seq data to help interpret mutations in a regulatory and splicing context

RepeatMasker is a program that screens DNA sequences for interspersed repeats and low complexity DNA sequences. The output of the program is a detailed annotation of the repeats that are present in the query sequence as well as a modified version of the query sequence in which all the annotated repeats have been masked (default: replaced by Ns). On average, almost 50% of a human genomic DNA sequence currently will be masked by the program.

RepeatModeler is a de novo transposable element (TE) family identification and modeling package. RepeatModeler assists in automating the runs of the various algorithms given a genomic database, clustering redundant results, refining and classifying the families and producing a high quality library of TE families suitable for use with RepeatMasker and ultimately for submission to the Dfam database (http://dfam.org).

REViewer is a tool for visualizing alignments of reads in regions containing tandem repeats. REViewer requires a BAMlet with graph-realigned reads generated by ExpansionHunter and the corresponding variant catalog.

Regulatory Genomics Toolbox: Python library and set of tools for the integrative analysis of high throughput regulatory genomics data. http://www.regulatory-genomics.org

RILseq computational protocol

R package: rliger. Linked Inference of Genomic Experimental Relationships. Use with rliger-R, rliger-Rscript.

RNA-SeQC is a java program which computes a series of quality control metrics for RNA-seq data.

Rockhopper is a comprehensive and user-friendly system for computational analysis of bacterial RNA-seq data. As input, Rockhopper takes RNA sequencing reads output by high-throughput sequencing technology (FASTQ, QSEQ, FASTA, SAM, or BAM files)

ROSE (Rank Ordering of Super-Enhancers) is tool for (1) creating stitched enhancers, and (2) separating super-enhancers from typical enhancers. given sequencing data (.bam) and a file of previously identified constituent enhancers (.gff)

Reciprocal Smallest Distance (RSD) is a pairwise orthology algorithm that uses global sequence alignment and maximum likelihood evolutionary distance between sequences to accurately detects orthologs between genomes.

RSEM is a software package for estimating gene and isoform expression levels from RNA-Seq data.

Rseqc comprehensively evaluate RNA-seq datasets generated from clinical tissues or other well annotated organisms such as mouse, fly and yeast.

variant detection for singletons, families, large pedigrees and populations, cancer, structural variant and CNV analysis, and microbial and metagenomic analysis

Rare Variant tests is a flexible software package for genetic association studies. It is designed to support unrelated individual or related (family-based) individuals

a tool for quantifying the expression of transcripts using RNA-seq data.

SalmonTE is an ultra-Fast and Scalable Quantification Pipeline of Transpose Element (TE) Abundances.

Sambamba is a high performance modern robust and fast tool (and library), written in the D programming language, for working with SAM and BAM files. Current parallelised functionality is an important subset of samtools functionality, including view, index, sort, markdup, and depth.

samblaster is a program for marking duplicates and finding discordant/split read pairs in read-id grouped paired-end SAM files. When marking duplicates, samblaster will use about 20MB per 1M read pairs. In a read-id grouped SAM file all alignments for a read-id (QNAME) are continuous. Aligners naturally produce such files. They can also be created by sorting a SAM file by read-id.

The samtools package now provides samtools, bcftools, tabix, and the underlying htslib library.

Scallop is a reference-based transcript assembler.

Bioinformatics pipeline for discovery of genetic variants from NGS reads.

Scanpy is a scalable toolkit for analyzing single-cell gene expression data. It includes preprocessing, visualization, clustering, pseudotime and trajectory inference and differential expression testing. The Python-based implementation efficiently deals with datasets of more than one million cells.

A graph neural network model to estimate cell-wise metabolic using single cell RNA-seq data

scHiCExplorer is a set of programs to process, normalize, analyse and visualize single-cell Hi-C data

Scramble is a mobile element insertion (MEI) detection tool. It identifies clusters of soft clipped reads in a BAM file, builds consensus sequences, aligns to representative L1Ta, AluYa5, and SVA-E sequences, and outputs MEI calls.

scvi-tools (single-cell variational inference tools) is a package for end-to-end analysis of single-cell omics data primarily developed and maintained by the Yosef Lab at UC Berkeley.

SeqAn is an open source C++ library of efficient algorithms and data structures for the analysis of sequences with the focus on biological data. It applies a unique generic design that guarantees high performance, generality, extensibility, and integration with other libraries. This package also contains a suite of apps, including Fiona, Gustaf, Mason, RazerS 3, Yara, SeqAn T-Coffee, Stellar, and searchjoin.

A cross-platform toolkit for FASTA/Q file manipulation

seqtk is a toolkit for processing sequences in FASTA/Q formats

Sequenza-utils is The supporting python library for the sequenza R package. Sequenza is a project the estimate purity/ploidy and copy number alteration from tumor sequencing experiments. Sequenza-utils provide command lines programs to transform common NGS file type, such as BAM, pileup and VCF, to input files for the R package

SHAPEIT is a fast and accurate haplotype inference software

De novo assembly from Oxford Nanopore reads

shmlast is a reimplementation of the Conditional Reciprocal Best Hits algorithm for finding potential orthologs between a transcriptome and a species-specific protein database. It uses the LAST aligner and the pydata stack to achieve much better performance while staying in the Python ecosystem.

SHRiMP is a software package for aligning genomic reads against a target genome. It was primarily developed with the multitudinous short reads of next generation sequencing machines in mind, as well as Applied Biosystem's colourspace genomic representation.

A clustering approach for identification of enriched domains from histone modification ChIP-Seq data

A windowed adaptive trimming tool for FASTQ files using quality

SIFT predicts whether an amino acid substitution affects protein function based on sequence homology and the physical properties of amino acids.

SlamDunk is a novel, fully automated software tool for automated, robust, scalable and reproducible SLAMseq data analysis.

The slapnap container is a tool for using the Compile, Analyze and Tally NAb Panels (CATNAP; Yoon et al. 2015) database to develop predictive models of HIV-1 neutralization sensitivity to one or several broadly neutralizing antibodies (bnAbs).

SMRT® Analysis is a bioinformatics software suite available for analysis of DNA sequencing data from Pacific Biosciences’ SMRT technology. Users can choose from a variety of analysis protocols that utilize PacBio® and third-party tools. Analysis protocols include de novo genome assembly, cDNA mapping, DNA base-modification detection, and long-amplicon analysis to determine phased consensus sequences.

Sniffles is a structural variation caller using third generation sequencing (PacBio or Oxford Nanopore). It detects all types of SVs (10bp+) using evidence from split-read alignments, high-mismatch regions, and coverage analysis.

snippy is a tool for rapid haploid variant calling and core genome alignment

SNP2HLA is a tool to impute amino acid polymorphisms and single nucleotide polymorphisms in human luekocyte antigenes (HLA) within the major histocompatibility complex (MHC) region in chromosome 6.

snpEff is a variant annotation and effect prediction tool. It annotates and predicts the effects of variants on genes (such as amino acid changes).

SNPTEST is a program for the analysis of single SNP association in genome-wide studies. The tests implemented include * Binary (case-control) phenotypes, single and multiple quantitative phenotypes * Bayesian and Frequentist tests * Ability to condition upon an arbitrary set of covariates * Various different methods for the dealing with imputed SNPs. The program is designed to work seamlessly with the output of both the genotype calling program CHIAMO, the genotype imputation program IMPUTE and the program GTOOL.

The purpose of this program is to identify single nucleotide positions that are different between tumor and normal (or, in theory, any two bam files). It takes a tumor bam and a normal bam and compares the two to determine the differences. It outputs a file in a format very similar to Samtools consensus format.

SortMeRNA is a local sequence alignment tool for filtering, mapping and clustering. The core algorithm is based on approximate seeds and allows for sensitive analysis of NGS reads. The main application of SortMeRNA is filtering rRNA from metatranscriptomic data. SortMeRNA takes as input a file of reads (fasta or fastq format) and one or multiple rRNA database file(s), and sorts apart aligned and rejected reads into two files specified by the user.

10x pipeline for processing Visium spatial RNA-seq data

SPAdes – St. Petersburg genome assembler – is intended for both standard isolates and single-cell MDA bacteria assemblies.

SpeedSeq is a genome analysis platform designed for rapid whole-genome variant detection and interpretation

Split-pool combinatorial barcoding makes it possible to scale projects to hundreds of samples and millions of cells, overcoming limitations of previous droplet based technologies. Spipe (split-pipe) implements combinatorial barcoding method for single cell RNA sequencing (scRNA-seq) with dramatically improved sensitivity.

The NCBI SRA Toolkit enables reading ("dumping") of sequencing files from the SRA database and writing ("loading") files into the .sra format.

SRST2 is a a read mapping-based tool for rapid molecular typing of bacterial pathogens. It allows fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment.

Short read aligner

Spliced Transcripts Alignment to a Reference

Transcript fusion detection

Tandem repeat expansion detection or genotyping from long-read alignments

STREAM stands for Single-cell Trajectories Reconstruction, Exploration And Mapping ofomics data. It is an interactive pipeline capable of disentangling and visualizing complex branching trajectories from both single-cell transcriptomic and epigenomic data.

Strelka is an analysis package designed to detect somatic SNVs and small indels from the aligned sequencing reads of matched tumor-normal samples.

StringTie is a fast and highly efficient assembler of RNA-Seq alignments into potential transcripts. It is primarily a genome-guided transcriptome assembler, although it can borrow algorithmic techniques from de novo genome assembly to help with transcript assembly.

STRipy-pipeline is non-graphical command-line version of STRipy that can be integrated into pipelines and analyse multiple STR or VNTR loci in parallel.

High-performance read alignment, quantification and mutation discovery

Fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions

SURVIVOR is a tool set for simulating/evaluating SVs, merging and comparing SVs within and among samples, and includes various methods to reformat or summarize SVs.

The svanna is an efficient and accurate pathogenicity prediction for coding and regulatory structural variants in long-read genome sequencing.

Svtyper is a Bayesian genotyper for structural variants.

svviz visualizes high-throughput sequencing data relevant to a structural variant. Only reads supporting the variant or the reference allele will be shown. svviz can operate in both an interactive web browser view to closely inspect individual variants, or in batch mode, allowing multiple variants (annotated in a VCF file) to be analyzed simultaneously.

Syri compares alignments between two chromosome-level assemblies and identifies synteny and structural rearrangements.

TALON is a Python package for identifying and quantifying known and novel genes/isoforms in long-read transcriptome data sets.

A Cross-platform and Efficient NCBI Taxonomy Toolkit

Single locus resolution of Transposable ELEment expression. Telescope estimates transposable element expression (retrotranscriptome) resolved to specific genomic locations. It directly addresses uncertainty in fragment assignment by reassigning ambiguously mapped fragments to the most probable source transcript as determined within a Bayesian statistical model.

TelSeq is a software that estimates telomere length from whole genome sequencing data (BAMs).

A package for including transposable elements in differential enrichment analysis of sequencing datasets.

Tumor Heterogeneity Analysis (THetA) is an algorithm used to estimate tumor purity and clonal/subclonal copy number aberrations simultaneously from high-throughput DNA sequencing data.

a suite of tools primarily for the identification of modified nucleotides from nanopore sequencing data.

TopHat is a fast splice junction mapper for RNA-Seq reads. It aligns RNA-Seq reads to mammalian-sized genomes using the ultra high-throughput short read aligner Bowtie, and then analyzes the mapping results to identify splice junctions between exons.

TPMCalculator quantifies mRNA abundance directly from the alignments by parsing BAM files

TranscriptClean is a Python program that corrects mismatches, microindels, and noncanonical splice junctions in long reads that have been mapped to the genome. It is designed for use with sam files from the PacBio Iso-seq and Oxford Nanopore transcriptome sequencing technologies

TransDecoder identifies candidate coding regions within transcript sequences, such as those generated by de novo RNA-Seq transcript assembly using Trinity, or constructed based on RNA-Seq alignments to the genome using Tophat and Cufflinks.

TransVar is a versatile annotator for 3-way conversion and annotation among genomic characterization(s) of mutations and transcript-dependent annotation(s).

trimAl is a tool for the automated removal of spurious sequences or poorly aligned regions from a multiple sequence alignment. It can consider several parameters, alone or in multiple combinations, in order to select the most-reliable positions in the alignment. These include the proportion of sequences with a gap, the level of residue similarity and, if several alignments for the same set of sequences are provided, the consistency level of columns among alignments. Moreover, trimAl is able to manually select a set of columns to be removed from the alignment.

Consistent quality and adapter trimming for RRBS or standard FastQ files.

Trimmomatic performs a variety of useful trimming tasks for illumina paired-end and single ended data.

Trinity, developed at the Broad Institute and the Hebrew University of Jerusalem, represents a novel method for the efficient and robust de novo reconstruction of transcriptomes from RNA-seq data.

Trinotate is a comprehensive annotation suite designed for automatic functional annotation of transcriptomes, particularly de novo assembled transcriptomes, from model or non-model organisms.

tRNAscan-SE 2.0 has advanced the state-of-the-art methodology in tRNA gene detection and functional prediction, captured by rich new content of the companion Genomic tRNA Database

Ultraplex is primarily designed for the demultiplexing of sequencing data generated using in-house library preparation protocols with custom adaptors

tools for dealing with Unique Molecular Identifiers (UMIs)/Random Molecular Tags (RMTs) and single cell RNA-Seq cell barcodes

VAtools is a python package that includes several tools to annotate VCF files with data from other tools.

A smarter, more reproducible, and more configurable tool for converting a VCF to a MAF.

annotate a VCF with other VCFs/BEDs/tabixed files

a simple C++ library for parsing and manipulating VCF files, + many command-line utilities

VCFtools contains a Perl API (Vcf.pm) and a number of Perl scripts that can be used to perform common tasks with VCF files such as file validation, file merging, intersecting, complements, etc.

Variant filtering using python expressions

VEP (Variant Effect Predictor) determines the effect of your variants (SNPs, insertions, deletions, CNVs or structural variants) on genes, transcripts, and protein sequence, as well as regulatory regions.

verifyBamID is a software that verifies whether the reads in particular file match previously known genotypes for an individual (or group of individuals), and checks whether the reads are contaminated as a mixture of two samples. verifyBamID can detect sample contamination and swaps when external genotypes are available. When external genotypes are not available, verifyBamID still robustly detects sample swaps.

Verkko is a hybrid genome assembly pipeline developed for telomere-to-telomere assembly of PacBio HiFi and Oxford Nanopore reads.

VIPER combines the use of several dozen RNA-seq tools, suites, and packages to create a complete pipeline that takes RNA-seq analysis from raw sequencing data all the way through alignment, quality control, unsupervised analyses, differential expression, and downstream pathway analysis

VirSorter2 is a DNA and RNA virus identification tool. It leverages genome-informed database advances across a collection of customized automatic classifiers to improve the accuracy and range of virus sequence detection.

Bioinformatics pipeline for viral transcriptome detection.

VSEARCH supports de novo and reference based chimera detection, clustering, full-length and prefix dereplication, rereplication, reverse complementation, masking, all-vs-all pairwise global alignment, exact and global alignment searching, shuffling, subsampling and sorting. It also supports FASTQ file analysis, filtering, conversion and merging of paired-end reads.

vt is a variant tool set that discovers short variants from Next Generation Sequencing data.

winnowmap is used for mapping ONT and PacBio reads to repetitive reference sequences.

Pipeline to process Xenium In Situ Gene Expression data

The HLA gene complex on human chromosome 6 is one of the most polymorphic regions in the human genome and contributes in large part to the diversity of the immune system. Accurate typing of HLA genes with short-read sequencing data has historically been difficult due to the sequence similarity between the polymorphic alleles. xHLA iteratively refines the mapping results at the amino acid level to achieve high typing accuracy for both class I and II HLA genes.

XHMM uses principal component analysis (PCA) normalization and a hidden Markov model (HMM) to detect and genotype copy number variation (CNV) from normalized read-depth data from targeted sequencing experiments.

Yet another K-mer analyzer. Yak is initially developed for two specific use cases: 1) to robustly estimate the base accuracy of CCS reads and assembly contigs, and 2) to investigate the systematic error rate of CCS reads.

A software platform for the analysis (including registration and interactive segmentation) and visualization (including volume rendering) of medical images and for research in image guided therapy.

AFNI (Analysis of Functional NeuroImages) is a set of C programs for processing, analyzing, and displaying functional MRI (FMRI) data - a technique for mapping human brain activity.

Advanced Normalization Tools (ANTs) extracts information from complex datasets that include imaging. Paired with ANTsR (answer), ANTs is useful for managing, interpreting and visualizing multidimensional data.

Alignment and Reconstruction for Electron Tomography

B-SOID (Behavioral Segmentation in Deeplabcut) is an unsupervised learning algorithm that serves to discover and classify behaviors that are not pre-defined by users. It segregates statistically different, sub-second rodent behaviors with a single bottom-up perspective video cameraR by performing a novel expectation maximization fitting of Gaussian mixture models on t-Distributed Stochastic Neighbor Embedding (t-SNE).

Baracus predicts brain age, based on data from Freesurfer. It combines data from cortical thickness, cortical surface area, and subcortical information

The ‘BrkRaw’ is a python module designed to provide a comprehensive tool to access raw data acquired from Bruker Biospin preclinical MRI scanner. This module is also compatible with the zip compressed data to enable use of the archived data directly. The module is comprised of four components, including graphical user interface (GUI), command-line tools, high-level and low-level python APIs.

Bsoft is a collection of programs and a platform for development of software for image and molecular processing in structural biology. Problems in structural biology are approached with a highly modular design, allowing fast development of new algorithms without the burden of issues such as file I/O. It provides an easily accessible interface, a resource that can be and has been used in other packages.

C3D is a command-line tool for converting 3D images between common file formats. The tool also includes a growing list of commands for image manipulation, such as thresholding and resampling.

A generalist algorithm for cellular segmentation with human-in-the-loop capabilities.

civet is a brain-imaging pipeline for analysis of large MR data sets. civet extracts and analyses cortical surfaces from MR images, as well as many other volumetric and corticometric functions.

CMTK is a Software toolkit for computational morphometry of biomedical images. CMTK provides a set of command line tools for processing and I/O.

CompuCell3D is a multiscale multicellular virtual tissue modeling and simulation environment. CompuCell3D is written in C++ and provides Python bindings for model and simulation development in Python. CompuCell3D is supported on Windows, Mac and Linux.

Tools to browse, download, explore, and analyze data from the Human Connectome Project (HCP). Allows users to compare their own data to that of the HCP.

CoolBox is an open-source, user-friendly toolkit for visual analysis of genomics data. It is highly compatible with the Python ecosystem and customizable with a well-designed user interface. It can bed used, for example, to produce high-quality genome track plots or fetch commonly used genomic data files with a Python script or command line, to explore genomic data interactively within Jupyter environment or web browser.

A configurable, open-source, Nipype-based, automated processing pipeline for resting state fMRI data.

The CTF MEG software has two main roles: - Provide a human-machine interface to the CTF MEG elec- tronics to collect MEG and/or EEG data. - Provide a tool for reviewing and (to a limited extent) ana- lyzing the MEG and/or EEG data acquired by the CTF MEG system.

Tool to convert dcm images to the neuroimaging BIDS format

DICOM to NIfTI converter

DeepCAD is a self-supervised deep-learning method for spatiotemporal enhancement of calcium imaging data that does not require any high signal-to-noise ratio (SNR) observations. DeepCAD suppresses detection noise and improves the SNR more than tenfold, which reinforces the accuracy of neuron extraction and spike inference and facilitates the functional analysis of neural circuits.

DeepLabCut is an open source toolbox that builds on a state-of-the-art human pose estimation algorithm. It allows training of a deep neural network by using limited training data to precisely track user-defined features, so that the human labeling accuracy will be matched.

This project aims to offer easy access to Deep Learning for segmentation of structures of interest in biomedical 3D scans. It is a system that allows the easy creation of a 3D Convolutional Neural Network, which can be trained to detect and segment structures if corresponding ground truth labels are provided for training. The system processes NIFTI images, making its use straightforward for many biomedical tasks.

Pipeline for preprocessing diffusion MRI datasets.

Dynamo is a software environment for subtomogram averaging of cryo-EM data.

a toolbox for rigid and nonrigid registration of images.

Fastsurfer is a neuroimaging pipeline based on deep learning.

Fiji is a distribution of the popular open-source software ImageJ focused on biological-image analysis. Fiji uses modern software engineering practices to combine powerful software libraries with a broad range of scripting languages to enable rapid prototyping of image-processing algorithms. Fiji facilitates the transformation of new algorithms into ImageJ plugins that can be shared with end users through an integrated update system.

Fitlins fits linear models to BIDS neuroimaging datasets.

A Robust Preprocessing Pipeline for fMRI Data

Freesurfer is a set of automated tools for reconstruction of the brain's cortical surface from structural MRI data, and overlay of functional MRI data onto the reconstructed surface.

FSL is a comprehensive library of image analysis and statistical tools for FMRI, MRI and DTI brain imaging data.

HiTIPS: High-Throughput Image Processing Software for the Study of Nuclear Architecture and Gene Expression. Documentations: https://hitips.readthedocs.io/en/latest/

Huygens is an image restoration, deconvolution, resolution and noise reduction. It can process images from all current optical microscopes, including wide-field, confocal, Nipkow (scanning disk confocal), multiple-photon, and 4Pi microscopes.

Imaris provides scientists with solutions for processing, visualizing and analyzing multi-dimensional microscopic images. It reads images in many of the most commonly used proprietary formats.

IMOD is a set of image processing, modeling and display programs used for tomographic reconstruction and for 3D reconstruction of EM serial sections and optical sections.

IsoNet is a deep learning-based software package that iteratively reconstructs the missing-wedge information and increases signal-to-noise ratio, using the knowledge learned from raw tomograms. Without the need for sub-tomogram averaging, IsoNet generates tomograms with significantly reduced resolution anisotropy.

ITK-SNAP is a tool for segmentation of 3D biomedical images. It requires a graphical connection to run on the cluster.

ITK-SNAP brings active contour segmentation to the fingertips of clinical researchers. This application fulfills a specific and pressing need of biomedical imaging research by providing a combination of manual and semiautomatic tools for extracting structures in 3D image data of different modalities and from different anatomical regions.

Tools to analyze layer fMRI datasets

libvips is a demand-driven, horizontally threaded image processing library. Compared to similar libraries, libvips runs quickly and uses little memory.

Given a set of labelled MR images (atlases) and unlabelled images (subjects), MAGeT produces a segmentation for each subject using a multi-atlas voting procedure based on a template library made up of images from the subject set.

Mango (Multi-image Analysis GUI) is a viewer for medical research images. It provides analysis tools and a user interface to navigate image volumes.

3D membrane segmentation for cryo-electron tomography.

This metaproject bundles multiple MINC-based packages that historically have been developed somewhat independently

The MIPAV (Medical Image Processing, Analysis, and Visualization) application enables quantitative analysis and visualization of medical images of numerous modalities such as PET, MRI, CT, or microscopy.

MONAI Label is a free and open-source platform that facilitates the development of AI-based applications that aim at reducing the time required to annotate 3D medical image datasets. It allows researchers to readily deploy their apps as services, which can be made available to clinicians via their preferred user-interface. Currently, MONAI Label readily supports locally installed (3DSlicer) and web-based (OHIF) frontends, and offers two Active learning strategies to facilitate and speed up the training of segmentation algorithms. MONAI Label allows researchers to make incremental improvements to their labeling apps by making them available to other researchers and clinicians alike.

Viewer for several types of brain scan formats

MRIQC is an MRI quality control tool

MRtrix provides a large suite of tools for image processing, analysis and visualisation, with a focus on the analysis of white matter using diffusion-weighted MRI.

MRtrix3Tissue is a fork of the MRtrix3 project. It aims to add capabilities for 3-Tissue CSD modelling and analysis to a complete version of the MRtrix3 software.

napari is an image viewer.Can be used to annotate images.

Preprocessing pipeline for neonate and infant MRI.

NiftyReg is an open-source software for efficient medical image registration. It has been mainly developed by members of the Centre for Medical Image Computing at University College London, UK.

nnU-Net is a PyTorch-based semantic segmentation application. It will analyze the provided training cases and automatically configure a matching U-Net-based segmentation pipeline.

OpenSlide is a C library for reading and manipulating digital slides of diverse vendor formats. It provides a simple interface to read whole-slide images (also known as virtual slides). OpenSlide has been used in the digital pathology projects.

PEET (Particle Estimation for Electron Tomography) is an open-source package for aligning and averaging particles in 3-D subvolumes extracted from tomograms. It seeks the optimal alignment of each particle against a reference volume through several iterations. If PEET and IMOD are both installed, most PEET operations are available from the eTomo graphical user interface in IMOD.

Positron Emission Tomography (PET) is a state-of-the-art neuroimaging tool for quantification of the in vivo spatial distribution of specific molecules in the brain. It is affected by various kinds of patient movement during a scan. The PETprep_HMC application allows for correction of the PET results in the presence of head motions.

Application for registration of medical images such as X-rays, CT, MRI and PET

PyTom is a software package for the analysis of volumetric data obtained by cryo electron tomography (cryo-ET). It covers a complete pipeline of processing steps for tomogram reconstruction, localization of macromolecular complexes in tomograms, fine alignment of subtomograms extracted at these locations, and their classification.

qsiprep configures pipelines for processing diffusion-weighted MRI (dMRI) data.

QuPath is open source software for bioimage analysis. It is often used for digital pathology applications because it offers a powerful set of tools for working with whole slide images - but it can be applied to lots of other kinds of image as well.

Rapidtide is a suite of Python programs used to model, characterize, visualize, and remove time varying, physiological blood signals from fMRI and fNIRS datasets. The primary workhorses of the package are the rapidtide program, which characterizes bulk blood flow, and happy, which focusses on the cardiac band.

Synthetic Aperture Magnetometry - The SANsrcV3 suite implements the latest advances in MEG source localization.

SimNIBS is a free software package for the Simulation of Non-invasive Brain Stimulation. It allows for realistic calculations of the electric field induced by transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS).

Cardiovascular Modeling and Simulation

Structural MRI PREProcessing (sMRIPrep) workflows for NIPreps (NeuroImaging PREProcessing tools).

SPHIRE (SPARX for High-Resolution Electron Microscopy) is an open-source, user-friendly software suite for the semi-automated processing of single particle electron cryo-microscopy (cryo-EM) data. It allows fast and reproducible structure determination from cryo-EM images.

The (S)tatistical (P)ara(M)etric application analyzes brain imaging data.

steinbock is a toolkit for processing multiplexed tissue images

Parallel computational fluid dynamics (CFD) finite element solver incorporating boundary conditions useful for blood flow simulations

Tedana is an application to denoise multi-echo fMRI datasets

TomoTwin - a deep metric learning based particle picking procedure for cryo-ET

(Tolerably Obsessive Registration and Tensor Optimization Indolent Software Ensemble) The TORTOISE software package is for processing diffusion MRI data.

TORTOISE (Tolerably Obsessive registration and Tensor Optimization Indolent Software Ensemble) is a suite of programs for for pre-processing, post-processing and analyzing diffusion MRI data

Tool for white matter bundle segmentation from Diffusion MRI.

xpcEngine performs denoising and estimation of Functional Connectivity on fMRI datasets

xcp_d is a postprocessing and noise regression pipeline for fMRI datasets (can use output from fmriprep and nibabies).

ADMIXTOOLS is a software package that supports formal tests of whether admixture occurred, and makes it possible to infer admixture proportions and dates.

ADMIXTURE is a software tool for maximum likelihood estimation of individual ancestries from multilocus SNP genotype datasets. It uses the same statistical model as STRUCTURE but calculates estimates much more rapidly using a fast numerical optimization algorithm.

high resolution HLA typing from RNA seq

Beagle is a package for imputing genotypes, inferring haplotype phase, and performing genetic association analysis. BEAGLE is designed to analyze large-scale data sets with hundreds of thousands of markers genotyped on thousands of samples.

BEAGLE is a high-performance library that can perform the core calculations at the heart of most Bayesian and Maximum Likelihood phylogenetics packages. It can make use of highly-parallel processors such as those in graphics cards (GPUs) found in many PCs.
module name: beagle-lib

BEAST (Bayesian Evolutionary Analysis Sampling Trees) is a cross-platform program for Bayesian MCMC analysis of molecular sequences.

bioBakery is a meta’omic analysis environment and collection of individual software tools with the capacity to process raw shotgun sequencing data into actionable microbial community feature profiles, summary reports, and publication-ready figures. It includes a collection of preconfigured analysis modules also joined into workflows for reproducibility. Each individual module has been developed to perform a particular task, e.g. quantitative taxonomic profiling or statistical analysis.

CD-HIT is a very widely used program for clustering and comparing protein or nucleotide sequences.

The EIGENSOFT package combines functionality from population genetics methods and EIGENSTRAT stratification correction method.

Progressive algorithm for large-scale multiple sequence alignments

In order to discover quantitative trait loci (QTLs), multi-dimensional genomic datasets combining DNA-seq and ChiP-/RNA-seq require methods that rapidly correlate tens of thousands of molecular phenotypes with millions of genetic variants while appropriately controlling for multiple testing. FastQTL implements a popular cis-QTL mapping strategy in a user- and cluster-friendly tool. FastQTL also proposes an efficient permutation procedure to control for multiple testing.

FastTree infers approximately-maximum-likelihood phylogenetic trees from alignments of nucleotide or protein sequences. FastTree can handle alignments with up to a million of sequences in a reasonable amount of time and memory. For large alignments, FastTree is 100-1,000 times faster than PhyML 3.0 or RAxML 7.

GCTA (Genome-wide Complex Trait Analysis) is designed to estimate the proportion of phenotypic variance explained by genome- or chromosome-wide SNPs for complex traits.

GenomeSyn, implemented in Perl, is a bioinformatics tool for visualizing genome synteny and structural variations. It provides genome synteny visualization for two or three genomic sequences. Additionally, genome annotation can be uploaded to visualize various SVs among different genomes.

GTDB-Tk is a software toolkit for assigning objective taxonomic classifications to bacterial and archaeal genomes based on the Genome Database Taxonomy

Gubbins is an algorithm that iteratively identifies loci containing elevated densities of base substitions while concurrently constructing a phylogeny based on the putative point mutations outside of these regions.

HyPhy (Hypothesis Testing using Phylogenies) is an open-source software package for the analysis of genetic sequences (in particular the inference of natural selection) using techniques in phylogenetics, molecular evolution, and machine learning.

Efficient phylogenomic software by maximum likelihood

KING is a toolset to explore genotype data from a genome-wide association study (GWAS) or a sequencing project. KING can be used to check family relationship and flag pedigree errors by estimating kinship coefficients and inferring IBD segments for all pairwise relationships.

The LTSOFT application implements a new approach to using information from known associated variants when conducting disease association studies. The approach is based in part on the classical technique of liability threshold modeling and performs estimation of model parameters for each known variant while accounting for the published disease prevalence from the epidemiological literature.

The Madeline 2.0 Pedigree Drawing Engine (PDE) is a pedigree drawing program for use in linkage and family-based association studies. The program is designed to handle large and complex pedigrees with an emphasis on readability and aesthetics.

MEGA, Molecular Evolutionary Genetics Analysis, is a software suite for analyzing DNA and protein sequence data from species and populations

Mega2 is a data-handling program for facilitating genetic linkage and association analyses.

MERLIN uses sparse trees to represent gene flow in pedigrees and is one of the fastest pedigree analysis packages around.

MetaWRAP is a modular pipeline for shotgun metagenomic data analysis. It deploys state-of-the-art software to handle metagenomic data processing starting from raw sequencing reads and ending in metagenomic bins and their analysis. It includes hybrid algorithms that leverage the strengths of a variety of software to extract and refine high-quality bins from metagenomic data through bin consolidation and reassembly.

mothur is a tool for analyzing 16S rRNA gene sequences generated on multiple platforms as part of microbial ecology projects.

MrBayes is a program for Bayesian inference and model choice across a wide range of phylogenetic and evolutionary models. MrBayes uses Markov chain Monte Carlo (MCMC) methods to estimate the posterior distribution of model parameters.

Ohana is a suite of software for analyzing population structure and admixture history using unsupervised learning methods.

OMA standalone is a standalone package that can infer orthologs using the OMA algorithm on custom genomes.

A package of programs for phylogenetic analyses of DNA and protein sequences using maximum likelihood.

infers haplotypes from population genotype data

Phylip is a package of programs for inferring phylogenies (evolutionary trees). Includes methods for parsimony, distance matrix and likelihood methods.

Pplacer places query sequences on a fixed reference phylogenetic tree to maximize phylogenetic likelihood or posterior probability according to a reference alignment. Pplacer is designed to be fast, to give useful information about uncertainty, and to offer advanced visualization and downstream analysis.

PRIMUS, Pedigree Reconstruction and Identification of a Maximum Unrelated Set, is used to read genome-wide IBD estimates and identify a maximum unrelated set and pedigree reconstruction.

PRSice is a Polygenic Risk Score software for calculating, applying, evaluating and plotting the results of polygenic risk scores (PRS) analyses.

QIIME is an open source software package for comparison and analysis of microbial communities, primarily based on high-throughput amplicon sequencing data (such as SSU rRNA) generated on a variety of platforms, but also supporting analysis of other types of data (such as shotgun metagenomic data).

RAxML (randomized axelerated maximum likelihood) is a sequential and parallel program for inference of large phylogenies with maximum likelihood (ML).

RAxML-NG is a phylogenetic tree inference tool which uses maximum-likelihood (ML) optimality criterion. Its search heuristic is based on iteratively performing a series of Subtree Pruning and Regrafting (SPR) moves, which allows to quickly navigate to the best-known ML tree. Successor to raxml.

Bayesian phylogenetic inference using probabilistic graphical models and an interpreted language

The scDRS application implements an approach that links scRNA-seq with polygenic disease risk at single-cell resolution, independent of annotated cell types. scDRS identifies cells exhibiting excess expression across disease-associated genes implicated by genome-wide association studies (GWASs). Genes whose expression was correlated with the scDRS score across cells (reflecting coexpression with GWAS disease-associated genes) were strongly enriched for gold-standard drug target and Mendelian disease genes.

SHAPEIT is a fast and accurate haplotype inference software

SMR integrates summary-level data from GWAS with data from expression quantitative trait locus (eQTL) studies to identify genes whose expression levels are associated with a complex trait because of pleiotropy. It implements methods to test for pleiotropic association between the expression level of a gene and a complex trait of interest using summary-level data from GWAS and expression quantitative trait loci (eQTL) studies (Zhu et al. 2016 Nat Genet).

SOLAR-Eclipse is an extensive, flexible software package for genetic variance components analysis, including linkage analysis, quantitative genetic analysis, SNP association analysis (QTN and QTLD), and covariate screening.

Summarize non-parameteric bootstrap or Bayesian posterior probability support for splits or clades on phylogenetic trees.

ensoorQTL leverages general-purpose libraries and graphics processing units (GPUs) to achieve high efficiency of computations at low costR. Using PyTorch or TensorFlow it allows > 200-fold decreases in runtime and ~ 5–10-fold reductions in cost when running on GPUs relative to CPUs.

TreeMix is a method for inferring the patterns of population splits and mixtures in the history of a set of populations.

TreeTime provides routines for ancestral sequence reconstruction and inference of molecular-clock phylogenies.

AlphaPept: a modern and open framework for MS-based proteomics

DIA-NN - a universal software for data-independent acquisition (DIA) proteomics data processing

diaTracer is a computational tool that enables spectrum-centric analysis of Bruker's diaPASEF data-independent acquisition proteomics data, facilitating direct (“spectral-library free”) peptide identification and quantification. Usually used via Fragpipe.

FragPipe is a Java Graphical User Interface (GUI) and CLI workflow tool for a suite of computational tools enabling comprehensive analysis of mass spectrometry-based proteomics data. It is powered by MSFragger.

IonQuant is a fast and comprehensive tool for MS1 precursor intensity-based quantification for timsTOF PASEF DDA and non-timsTOF (e.g., Orbitrap) data. It enables label-free quantification with false discovery (FDR) controlled match-between-runs (MBR). It can also be used for quantification in labelling-based experiments such as those involving SILAC, dimethyl, or similar labelling strategies. IonQuant is available as part of FragPipe (recommended option), but can also be run as a command-line tool.

The Mascot search engine uses mass spectrometry data to identify proteins from primary sequence databases. Mascot searches can be run directly on the NIH Mascot server at https://biospec.nih.gov, or by using the Mascot daemon on your own desktop PC.

MaxQuant is a quantitative proteomics software package designed for analyzing large mass-spectrometric data sets. It is specifically aimed at high-resolution MS data. Several labeling techniques as well as label-free quantification are supported.

An ultrafast database search tool for peptide identification in mass spectrometry-based proteomics.

A software for postprocessing of shotgun proteomics data

Philosopher is fast, easy-to-use, scalable, and versatile data analysis software for mass spectrometry-based proteomics. Philosopher is dependency-free and can analyze both traditional database searches and open searches for post-translational modification (PTM) discovery.

Advanced Normalization Tools (ANTs) extracts information from complex datasets that include imaging. Paired with ANTsR (answer), ANTs is useful for managing, interpreting and visualizing multidimensional data.

Command line interface to stan

DENTIST (Detecting Errors iN analyses of summary staTISTics) is a quality control (QC) tool for summary-level data from genome-wide association studies (GWASs). It leverages the difference between the observed GWAS test-statistic of a variant and its predicted value (using the neighbouring variants and linkage equilibrium (LD) data from a reference panel) to remove problematic variants.

FEBio software suite implement a nonlinear implicit finite element (FE) framework, designed specifically for analysis in computational solid biomechanics. FEBio offers modeling scenarios, constitutive models, and boundary conditions, which are relevant to numerous applications in biomechanics. The open-source FEBio software is written in C++, with particular attention to scalar and parallel performance on modern computer architectures.

The GAUSS Mathematical and Statistical System is an easy-to-use data analysis environment based on the fast and powerful GAUSS Matrix Programming Language designed for computationally intensive tasks.

Gurobi is a mathematical optimization solver. It is a commercial product developed by gurobi.com. On Biowulf, Gurobi is licensed for use by the members of the CDSL_Gurobi_users group only. It is installed in /data/CDSL_Gurobi_users and is not accessible by any other users. A token license server, running on Biowulf, manages the Gurobi license.

IDL and ENVI are a complete computing environment for the interactive analysis and visualization of data. IDL integrates an array-oriented language with mathematical analysis and graphical display techniques. ENVI is designed for extracting information from geospatial and medical imagery.

Lifelines is a complete survival analysis library, written in pure Python. It has the benefits of easy installation, internal plotting methods, simple and intuitive API handles right, left and interval censored data. It also contains the most popular parametric, semi-parametric and non-parametric models.

Mathematica is an interactive system for doing mathematical computation. It performs numerical, symbolic and graphical computations, and incorporates a high-level programming language.

MATLAB is an interactive software package for scientific and engineering numeric computation. MATLAB integrates numerical analysis, matrix computation, signal processing, and graphics in an environment where problems and solutions are expressed just as they are written mathematically.

MCL implements Markov cluster algorithm. Among its applications is the assignment of proteins into families based on precomputed sequence similarity information. This approach does not suffer from the problems that normally hinder other protein sequence clustering algorithms, such as the presence of multi-domain proteins, promiscuous domains and fragmented proteins.

Meep (or MEEP) is a free finite-difference time-domain (FDTD) simulation software package developed at MIT to model electromagnetic systems, along with the MPB eigenmode package.

GNU Octave is a high-level language, primarily intended for numerical computations. It provides a convenient command line interface for solving linear and nonlinear problems numerically, and for performing other numerical experiments using a language that is mostly compatible with Matlab.

PEER stands for probabilistic estimation of expression residuals. It is a collection of Bayesian approaches to infer hidden determinants and their effects from gene expression profiles using factor analysis methods.

R (the R Project) is a language and environment for statistical computing and graphics. R is similar to S, and provides a wide variety of statistical and graphical techniques (linear and nonlinear modelling, statistical tests, time series analysis, classification, clustering, ...).

RStudio is a set of integrated tools designed to help you be more productive with R. It includes a console, syntax-highlighting editor that supports direct code execution, as well as tools for plotting, history, debugging and workspace management.

Base SAS provides a scalable, integrated software environment specially designed for data access, transformation and reporting.

"AutoDock CrankPep or ADCP is an AutoDock docking engine specialized for docking peptides. It combines technology form the protein folding filed with an efficient representation of a rigid receptor as affinity grids to fold the peptide in the context of the energy landscape created by the receptor."

Blender is the free and open source 3D creation suite. Blender on Biowulf is meant for command-line rendering.

CHAP is a tool for the functional annotation of ion channel structures

Chimera is a highly extensible program for interactive visualization and analysis of molecular structures and related data, including density maps, supramolecular assemblies, sequence alignments, docking results, trajectories, and conformational ensembles.

UCSF ChimeraX (or simply ChimeraX) is the next-generation molecular visualization program from the Resource for Biocomputing, Visualization, and Informatics (RBVI), following UCSF Chimera.

Chromeister is an ultra fast, heuristic approach to detect conserved signals in extremely large pairwise genome comparisons.

Coot is for macromolecular model building, model completion and validation, particularly suitable for protein modelling using X-ray data.

Cytoscape is an open source software platform for visualizing molecular interaction networks and biological pathways and integrating these networks with annotations, gene expression profiles and other state data.

free open-source molecular editor and visualization package. It offers a range of features including a molecular editor, surface generation (orbitals and densities) and animations (vibrational modes and reaction pathways).

LAMMPS is a classical molecular dynamics code, and an acronym for Large-scale Atomic/Molecular Massively Parallel Simulator. It runs on a variety of different computer systems, including single processor systems, distributed-memory machines with MPI, and GPU and Xeon Phi systems. LAMMPS is open source software, released under the GNU General Public License.

A collection of tools for characterizing membrane structure and dynamics within simulated trajectories of molecular systems.

Open Babel is a chemical toolbox designed to speak the many languages of chemical data.

Automates many of the common tasks of preparing structures for continuum solvation calculations as well as many other types of biomolecular structure modeling

PoseFilter is a PyMOL plugin and assists in the analysis of docked ligands through identification of unique oligomeric poses by utilizing RMSD and interaction fingerprint analysis methods.

Psi4 is an ab-initio electronic structure code that supports various methods for calculating energies and gradients of molecular systems.

The Rosetta++ software suite can perform de novo protein structure predictions, identify low free energy sequences for target protein backbones, predict the structure of a protein-protein complex from the individual structures of the monomer components, incorporate NMR data into the basic Rosetta protocol to accelerate the process of NMR structure prediction, and more...

A limited number of Schrödinger applications are available on the Biowulf cluster through the Molecular Modeling Interest Group. Most are available through the Maestro GUI.

VMD is a molecular visualization program for displaying, animating, and analyzing large biomolecular systems using 3-D graphics and built-in scripting. To use, type vmd at the prompt.

a tool for genotyping Variable Number Tandem Repeats (VNTR) from sequence data

Another Gtf/Gff Analysis Toolkit

Calculates genotype frequencies of a SNPMatrix. This component tests each SNP for its Hardy-Weinberg equilibrium. If there are NA values, the frequencies of missing value per sample in the input file are calculated.

AMRFinderPlus - Identify AMR genes and point mutations, and virulence and stress resistance genes in assembled bacterial nucleotide and protein sequence.

Processing and integrating RNA-Seq data in order to generate high-resolution annotations is challenging, time consuming and requires numerous different steps. ANNOgesic is a powerful and modular pipeline that provides the required analyses and simplifies RNA-Seq-based bacterial and archaeal genome annotation. It predicts and annotates numerous features, including small non-coding RNAs, with high precision.

antiSMASH allows the rapid genome-wide identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genomes.

Anvi’o is an open-source, community-driven analysis and visualization platform for microbial ‘omics. It brings together many aspects of today’s cutting-edge strategies including genomics, metagenomics, metatranscriptomics, pangenomics, metapangenomics, phylogenomics, and microbial population genetics in an integrated and easy-to-use fashion through extensive interactive visualization capabilities.

high resolution HLA typing from RNA seq

Arriba identifies gene fusions in RNA-Seq data. It also can detect other structural variants in genomic data, such as intron duplications and gene truncations.

AUGUSTUS is a program that predicts genes in eukaryotic genomic sequences.

Rapid & standardized annotation of bacterial genomes, MAGs & plasmids

Bartender is an accurate clustering algorithm to detect barcodes and their abundances from raw next-generation sequencing data. In contrast with existing methods that cluster based on sequence similarity alone, Bartender uses a modified two-sample proportion test that also considers cluster size. This modification results in higher accuracy and lower rates of under- and over-clustering artifacts.

bioBakery is a meta’omic analysis environment and collection of individual software tools with the capacity to process raw shotgun sequencing data into actionable microbial community feature profiles, summary reports, and publication-ready figures. It includes a collection of preconfigured analysis modules also joined into workflows for reproducibility. Each individual module has been developed to perform a particular task, e.g. quantitative taxonomic profiling or statistical analysis.

Birdsuite is a four-stage analytical framework instantiated in software for deriving integrated and mutually consistent copy number and SNP genotypes. The method sequentially assigns copy number across regions of common copy number polymorphisms (CNPs), calls genotypes of SNPs, identifies rare CNVs via a hidden Markov model (HMM), and generates an integrated sequence and copy number genotype at every locus.

NCBI's well-known sequence database searching program which compares a nucleotide or protein query sequence against all sequences in a database.

BLAT is a DNA/Protein Sequence Analysis program that is designed to quickly find sequences of 95% and greater similarity of length 40 bases or more.

The BOLT-LMM algorithm computes statistics for testing association between phenotype and genotypes using a linear mixed model (LMM)

Bracken is a companion program to Kraken 1, KrakenUniq, or Kraken 2 While Kraken classifies reads to multiple levels in the taxonomic tree, Bracken allows estimation of abundance at a single level using those classifications (e.g. Bracken can estimate abundance of species within a sample).

The next version of the bwa-mem algorithm in bwa.

Cactus is a reference-free whole-genome multiple alignment program.

CADD (Combined Annotation Dependent Depletion) is a tool for scoring the deleteriousness of single nucleotide variants as well as insertion/deletions variants in the human genome. Currently, it supports the builds: GRCh37/hg19 and GRCh38/hg38.

CAVIAR (CAusal Variants Identication in Associated Regions) is a statistical framework that quantifies the probability of each variant to be causal while allowing with arbitrary number of causal variants

Miscellaneous collection of Python and R scripts for processing Iso-Seq data

Cell2location: Comprehensive mapping of tissue cell architecture via integrated single cell and spatial transcriptomics (cell2location model)

Efficient genotyping bi-allelic SNPs on single cells

BigWig manpulation tools using libBigWig and htslib

Rapid assessment of genome bin quality using machine learning.

The CHESS (Comparison of Hi-C Experiments using Structural Similarity) application implements an algorithm for the comparison of chromatin contact maps and automatic differential feature extraction.

AQUAS Transcription Factor and Histone ChIP-Seq processing pipeline. The AQUAS pipeline is based off the ENCODE (phase-3) transcription factor and histone ChIP-seq pipeline specifications (by Anshul Kundaje)

CIRCLE-seq is a tool ifor analysis of circularization for in vitro reporting of cleavage effects by sequencing. It is a highly sensitive, sequencing-efficient in vitro screening strategy that outperforms existing cell-based or biochemical approaches for identifying CRIScrisPR–Cas9 genome-wide off-target mutations. It can be practiced using widely accessible next-generation sequencing technology and does not require reference genome sequences. Importantly, CIRCLEcircle-seq can be used to identify off-target mutations associated with cell-type-specific single-nucleotide polymorphisms, demonstrating the feasibility and importance of generating personalized specificity profiles.

Circos is a program for the generation of publication-quality, circularly composited renditions of genomic data and related annotations. Circos is particularly suited for visualizing alignments, conservation and intra and inter-chromosomal relationships. Also, Circos is useful to visualize any type of information that benefits from a circular layout. Thus, although it has been designed for the field of genomics, it is sufficiently flexible to be used in other data domains.

Clair3 is a small variant caller for Illumina, PacBio and ONT long reads. Compare to PEPPER (r0.4), Clair3 (v0.1) shows a better SNP F1-score with ≤30-fold of ONT data (precisionFDA Truth Challenge V2), and a better Indel F1-score, while runs generally four times faster.

cLoops2 is an enhanced and elegant flexible peak/loop/domain -calling and analysis tool for 1D/3D genomic data.

Clustal-Omega is a general purpose multiple sequence alignment (MSA) program for proteins and DNA/RNA. It produces high quality MSAs and is capable of handling data-sets of hundreds of thousands of sequences in reasonable time.

ClustalW is a general-purpose multiple alignment program for DNA or protein sequences.

CoNIFER (copy number inference from exome reads) uses exome sequencing data to find copy number variants (CNVs) and genotype the copy-number of duplicated genes. It can reliably be used to discover disruptive genic CNVs missed by standard approaches and should have broad application in human genetic studies of disease.

Cooltools is a suite of computational tools that enables flexible, scalable, and reproducible analysis of high-resolution contact frequency data. Cooltools leverages the widely-adopted cooler format which handles storage and access for high-resolution datasets. Cooltools provides a paired command line interface and Python application programming interface, which respectively facilitate workflows on high-performance computing clusters and in interactive analysis environments.

CoverageMaster (CoM) is a copy number variation (CNV) calling algorithm i ased on depth-of-coverage maps designed to detect CNVs of any size i n exome [whole exome sequencing (WES)] and genome [whole genome sequencing (WGS)] data. The core of the algorithm is the compression of sequencing coverage data in a multiscale Wavelet space and the analysis through an iterative Hidden Markov Model.

CutRunTools is a flexible, general pipeline for facilitating the identification of chromatin-associated protein binding and genomic footprinting analysis from antibody-targeted CutRun primary cleavage data. CutRunTools extracts endonuclease cut site information from sequences of short read fragments and produces single-locus binding estimates, aggregate motif footprints, and informative visualizations to support the high-resolution mapping capability of CutRun.

CytoSig is a data-driven infrastructure hosted by the National Cancer Institute. CytoSig includes both a database of target genes modulated by cytokines and a predictive model of cytokine signaling activity and regulatory cascade from transcriptomic profiles.

CytoSPACE implements an optimization method for mapping individual cells from a single-cell RNA sequencing atlas to spatial expression profiles. Across diverse platforms and tissue types, it outperforms previous methods with respect to noise tolerance and accuracy, enabling tissue cartography at single-cell resolution.

CytoTRACE (Cellular (Cyto) Trajectory Reconstruction Analysis using gene Counts and Expression) is a computational method that predicts the differentiation state of cells from single-cell RNA-sequencing data. CytoTRACE leverages a simple, yet robust, determinant of developmental potential—the number of detectably expressed genes per cell, or gene counts. We have validated CytoTRACE on ~150K single-cell transcriptomes spanning 315 cell phenotypes, 52 lineages, 14 tissue types, 9 scRNA-seq platforms, and 5 species.

DeCiFer is an algorithm that simultaneously selects mutation multiplicities and clusters somatic single-nucleotide variants (SNVs) by their corresponding descendant cell fractions (DCF), a statistic that quantifies the proportion of cells which acquired the SNV or whose ancestors acquired the SNV. DCF is related to the commonly used cancer cell fraction (CCF) but further accounts for SNVs which are lost due to deleterious somatic copy-number aberrations (CNAs), identifying clusters of SNVs which occur in the same phylogenetic branch of tumour evolution.

DeepLoc 2.0 predicts the subcellular localization(s) of eukaryotic proteins. It is is able to predict one or more localizations for any given protein.

DeepPolisher is a transformer-based sequencing correction method similar to DeepConsensus. DeepPolisher is designed to identify errors in genome assemblies. DeepPolisher takes haplotype-specific reads aligned to phased assemblies and produces a VCF file containing potential errors in the assembly.

DeepSEA is a deep learning-based algorithmic framework for predicting the chromatin effects of sequence alterations with single nucleotide sensitivity. DeepSEA can accurately predict the epigenetic state of a sequence, including transcription factors binding, DNase I sensitivities and histone marks in multiple cell types, and further utilize this capability to predict the chromatin effects of sequence variants and prioritize regulatory variants.

DeepSomatic is a google developed analysis pipeline that uses a deep neural network to call somatic variants from tumor-normal sequencing data.

DeNovoGear is az software for analyzing de novo mutations from familial and somatic tissue sequencing data. It uses likelihood-based error modeling to reduce the false positive rate of mutation discovery in exome analysis and fragment information to identify the parental origin of germ-line mutations. DeNovoGear has been used on human whole-genome sequencing data to produce a set of predicted de novo insertion and/or deletion (indel) mutations.

DEPICT is an integrative tool that employs predicted gene functions to systematically prioritize the most likely causal genes at associated loci, highlight enriched pathways and identify tissues/cell types where genes from associated loci are highly expressed. DEPICT is not limited to genes with established functions and prioritizes relevant gene sets for many phenotypes.

DIAMOND is a new high-throughput program for aligning DNA reads or protein sequences against a protein reference database such as NR, at up to 20,000 times the speed of BLAST, with high sensitivity.

dRep is a python program for rapidly comparing large numbers of genomes. dRep can also "de-replicate" a genome set by identifying groups of highly similar genomes and choosing the best representative genome for each genome set.

Duphold: scalable, depth-based annotation and curation of high-confidence structural variant calls

dyno is a meta package that installs several other packages from the dynvers (https://github.com/dynverse). It comprises a set of R packages to construct and interpret single-cell trajectories.

Eagle performs a reference-based haplotype phasing. It attains high accuracy across a broad range of cohort sizes by efficiently leveraging information from large external reference panels (such as the Haplotype Reference onsortium; HRC) using a new data structure based on the positional Burrows-Wheeler transform.

ESPRESSO is a novel method for processing alignment of long read RNA-seq data, which can effectively improve splice junction accuracy and isoform quantification. ESPRESSO jointly considers alignments of all long reads aligned to a gene and uses error profiles of individual reads to improve the identification of splice junctions and the discovery of their corresponding transcript isoforms.

Exonerate is a generic tool for pairwise sequence comparison. It allows you to align sequences using a many alignment models, either exhaustive dynamic programming or a variety of heuristics.

Expansion Hunter: a tool for estimating repeat sizes. There are a number of regions in the human genome consisting of repetitions of short unit sequence (commonly a trimer). Such repeat regions can expand to a size much larger than the read length and thereby cause a disease. Expansion Hunter aims to estimate sizes of such repeats by performing a targeted search through a BAM/CRAM file for reads that span, flank, and are fully contained in each repeat.

ExpansionHunter Denovo (EHdn) is a suite of tools for detecting novel expansions of short tandem repeats (STRs). EHdn is intended for analysis of a collection of BAM/CRAM files containing alignments of short (100-200bp) reads.

A tool designed to provide fast all-in-one preprocessing for FastQ files. This tool is developed in C++ with multithreading supported to afford high performance.

Fastq-pair rewrites paired end fastq files to make sure that all reads have a mate and to separate out singletons.

fastq_demux is a simple program to demultiplex a FASTQ file or a pair of FASTQ files based on the barcodes present in the FASTQ headers.

FitHiChIP is a computational method for identifying chromatin contacts among regulatory regions such as enhancers and promoters from HiChIP/PLAC-seq data. FitHiChIP jointly models the non-uniform coverage and genomic distance scaling of HiChIP data, captures previously validated enhancer interactions for several genes including MYC and TP53, and recovers contacts genome-wide that are supported by ChIA-PET, promoter capture Hi-C and Hi-C data.

FLAIR (Full-Length Alternative Isoform analysis of RNA) is a workflow leveraging the full-length transcript sequencing data that nanopore affords. It uses multiple alignment steps and splice site filters to increase confidence in the set of isoforms defined from noisy data.

FOCUS (Fine-mapping Of CaUsal gene Sets) is software to fine-map transcriptome-wide association study statistics at genomic risk regions. The software takes as input summary GWAS data along with eQTL weights and outputs a credible set of genes to explain observed genomic risk.

FreePSI is a new method for genome-wide percent spliced in (PSI) estimation that requires neither a reference transcriptome (hence, transcriptome-free) nor the mapping of RNA-seq reads (hence, alignment-free). The first attribute allows FreePSI to work effectively when a high quality reference transcriptome is unavailable and the second not only helps make FreePSI more efficient, it also eliminates the necessity of dealing with multi-reads.

Funannotate is a genome prediction, annotation, and comparison software package. It was originally written to annotate fungal genomes (small eukaryotes ~ 30 Mb genomes), but has evolved over time to accomodate larger genomes.

FuSeq is a software for discovering fusion genes from paired-end RNA sequencing data. It implements a fast and accurate method to discover fusion genes based on quasi-mapping to quickly map the reads, extract initial candidates from split reads and fusion equivalence classes of mapped reads, and finally apply multiple filters and statistical tests to get the final candidates.

In silico Validation of Fusion Transcript Predictions

GEM (Gene-Environment interaction analysis for Millions of samples) is a software program for large-scale gene-environment interaction testing in samples from unrelated individuals. It enables genome-wide association studies in up to millions of samples while allowing for multiple exposures, control for genotype-covariate interactions, and robust inference.

collection of bioinformatic tools

GenoTools is a suite of automated genotype data processing steps written in Python. The core pipeline was built for Quality Control and Ancestry estimation of data in the Global Parkinson's Genetics Program (GP2).

gfatools is a set of tools for manipulating sequence graphs in the GFA or the rGFA format. It has implemented parsing, subgraph and conversion to FASTA/BED.

gffcompare can be used to compare and evaluate the accuracy of RNA-Seq transcript assemblers (Cufflinks, Stringtie). It can collapse (merge) duplicate transcripts from multiple GTF/GFF3 files (e.g. resulted from assembly of different samples) and classify transcripts from one or multiple GTF/GFF3 files as they relate to reference transcripts provided in a annotation file (also in GTF/GFF3 format).

GimmeMotifs is a pipeline for transcription factor motif analysis written in Python. It incorporates an ensemble of computational tools to predict motifs de novo from ChIP-sequencing data. Similar redundant motifs are compared using the weighted information content similarity score and clustered using an iterative procedure. A comprehensive output report is generated with several different evaluation metrics to compare and evaluate the results.

GLIMPSE is a phasing and imputation method for large-scale low-coverage sequencing studies. It perform accurate imputed genotype calls and outperforms SNP arrays.

Joint variant calling for large cohort sequencing

G-Nomix does high resolution ancestry deconvolution for high-throughput sequencing data.

GraphMap is a highly sensitive and accurate mapper for long, error-prone reads. It offers a number of valuable features, such as mapping position agnostic to alignment parameters, high sensitivity and precision, handling circular genomes, meaningful mapping quality, various alignment strategies, and more.

GTDB-Tk is a software toolkit for assigning objective taxonomic classifications to bacterial and archaeal genomes based on the Genome Database Taxonomy

gutSMASH is a tool that has been developed to systematically evaluate the metabolic potential of anaerobic bacteria in the gut by predicting both known and novel anaerobic metabolic gene clusters from the gut microbiome.

The hap-ibd program detects identity-by-descent (IBD) segments and homozygosity-by-descent (HBD) segments in phased genotype data. The hap-ibd program can analyze data sets with hundreds of thousands of samples.

HemTools is a collection of NGS pipelines and bioinformatic analysis tools. It includes tools for data visualization, motif analysis, integrative analysis, bioinformatica analysis, differential analysis, CRISPR analysis, and more.

The HH-suite is an open-source software package for sensitive protein sequence searching based on the pairwise alignment of hidden Markov models (HMMs).

hicmaptools is a command line tool to access HiC maps. The complete program provides multi-query modes and analysis tools.

Pacific Biosciences workflow for analyzing human PacBio whole genome sequencing data.

Tool for genotyping short tandom repeats from Illumina sequencing data

HMMER is used for searching sequence databases for homologs of protein sequences, and for making protein sequence alignments. It implements methods using probabilistic models called "profile hidden Markov models" (profile HMMs). Compared to BLAST, FASTA, and other sequence alignment and database search tools based on older scoring methodology, HMMER aims to be significantly more accurate and more able to detect remote homologs because of the strength of its underlying mathematical models.

A pipeline for comprehensive analysis of HTGTS-Rep-seq.

iDiffIR is a tool for identifying differential IR from RNA-seq data. It accepts any sorted, indexed BAM file for single- or paired-end reads.

IgBlast is a sequence analysis tool for immunoglobulin variable domains.

The Immcantation framework provide a start-to-finish analytical ecosystem for high-throughput AIRR-seq datasets. Immcantation includes Alakazam, Change-O, Dowser, pRESTO, SCOPer, SHazaM, and TIgGER.

ImReP is a novel computational method for rapid and accurate profiling of the adaptive immune repertoire from regular RNA-Seq data. It is able to efficiently extract TCR- and BCR-derived reads from RNA-Seq data. ImReP can also accurately assemble the complementary determining regions 3 (CDR3s), the most variable regions of B and T cell receptors, and determine their antigen specificity.

Package for searching DNA sequence databases for RNA structure and sequence similarities

InterProScan is the software package that allows sequences (protein and nucleic) to be scanned against InterPro's signatures. Signatures are predictive models, provided by several different databases, that make up the InterPro consortium.

In 2015, the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) published updated standards and guidelines for the clinical interpretation of sequence variants with respect to human diseases on the basis of 28 criteria. However, variability between individual interpreters can be extensive because of reasons such as the different understandings of these guidelines and the lack of standard algorithms for implementing them, yet computational tools for semi-automated variant interpretation are not available. To address these problems, InterVar implements these criteria to help human reviewers interpret the clinical significance of variants. InterVar can take a pre-annotated or VCF file as input and generate automated interpretation on 18 criteria.

Identifies cancer genes and pinpoints their putative mechanism of action across tumor types

iVar is a computational package that contains functions broadly useful for viral amplicon-based sequencing. Additional tools for metagenomic sequencing are actively being incorporated into iVar. While each of these functions can be accomplished using existing tools, iVar contains an intersection of functionality from multiple tools that are required to call iSNVs and consensus sequences from viral sequencing data across multiple replicates.

JAX-CNV implemens an algorithm for copy number variant (CNV) calling from the whole-genome sequencing (WGS) data. On testing data, it demonstrated ~7-fold increase in the number of detected CNVs as compared to the chromosomal microarray assay (CMA) for clinical diagnosis.

Jellyfish is a tool for fast, memory-efficient counting of k-mers in DNA.

Kalign is a fast and accurate multiple sequence alignment algorithm designed to align large numbers of protein sequences.

kallisto is a program for quantifying abundances of transcripts from RNA-Seq data, or more generally of target sequences using high-throughput sequencing reads. It is based on the novel idea of pseudoalignment for rapidly determining the compatibility of reads with targets, without the need for alignment.

LAST is designed for moderately large data (e.g. genomes, DNA reads, proteomes). It's especially geared toward:

• Finding rearrangements and recombinations (last-split)
• Finding DNA-versus-protein related regions, especially protein fossils.
• Unusual data, e.g. AT-rich DNA, because it can fit parameters to the data and calculate significance.
• Sensitive DNA-DNA search, due to fitting, sensitive seeding, and calculating significance.

LASTZ is a tool for (1) aligning two DNA sequences, and (2) inferring appropriate scoring parameters automatically. LASTZ is a drop-in replacement for BLASTZ, and is backward compatible with BLASTZ's command-line syntax.

ldsc is a command line tool for estimating heritability and genetic correlation from GWAS summary statistics. ldsc also computes LD Scores.

Liger2LiGer is a Nanopore chimera splitting/detection tool. Automated end-to-end chimera detection and dataset evaluation starting from fastq.

The La Jolla Assembler (LJA) is a tool for assemblies of long and accurate reads. It reduces the error rate in these reads by three orders of magnitude (making them nearly error-free) and constructs the de Bruijn graph for large genomes and large k-mer sizes. Since the de Bruijn graph constructed for a fixed k-mer size is typically either too tangled or too fragmented, LJA uses a new concept of a multiplex de Bruijn graph with varying k-mer sizes.

LoMA is localized assembly method that constructs highly accurate consensus sequences (CSs) from long reads. LoMA's algorithm employs minimap2 and MAFFT and classifies diploid haplotypes based on structural variants and CSs. It allows analyzing human samples sequenced with the Oxford Nanopore sequencer.

A selection of software developed at the Australian Centre for Ecogenomics to aid in the analysis of metagenomic datasets: unitem, refinem, checkm, graftm, groopm, bamm, finishm, singlem, orfm, and coverm

FOCUS (Fine-mapping Of CaUsal gene Sets) is software to fine-map transcriptome-wide association study statistics at genomic risk regions. The software takes as input summary GWAS data along with eQTL weights and outputs a credible set of genes to explain observed genomic risk.

MAGMA is a tool for gene analysis and generalized gene-set analysis of GWAS data. It can be used to analyse both raw genotype data as well as summary SNP p-values from a previous GWAS or meta-analysis.

Microsatellite Analysis for Normal-Tumor InStability is a program developed for detecting microsatellite instability from paired-end BAM files.

a set of multiple scripts used to analyze PLAC-Seq and HiChIP data.

MashMap is an approximate algorithm for computing local alignment boundaries between long DNA sequences. Given a minimum alignment length and an identity threshold, it computes the desired alignment boundaries and identity estimates using kmer-based statistics, and maintains sufficient probabilistic guarantees on the output sensitivity.

MELT is an application for identifying mobile elements in genomic data

MEME is used to discover motifs in groups of DNA/protein sequences or databases.

tool to combine genotype data imputed against multiple reference panels.

MetaWRAP is a modular pipeline for shotgun metagenomic data analysis. It deploys state-of-the-art software to handle metagenomic data processing starting from raw sequencing reads and ending in metagenomic bins and their analysis. It includes hybrid algorithms that leverage the strengths of a variety of software to extract and refine high-quality bins from metagenomic data through bin consolidation and reassembly.

MINTIE is a tool for identifying novel, rare transcriptional variants in cancer RNA-seq data. MINTIE detects gene fusions, transcribed structural variants, novel splice variants and complex variants, and annotates all novel transcriptional variants.

MMseqs2 (Many-against-Many sequence searching) is a software suite to search and cluster huge protein and nucleotide sequence sets.

MoChA is a bcftools extension to call mosaic chromosomal alterations starting from phased VCF files with either B Allele Frequency (BAF) and Log R Ratio (LRR) or allelic depth (AD)

evaluates Microsatellite Instability (MSI) for cancer patients with next generation sequencing data.

Mummer is a system for aligning entire genomes extremely rapidly.

Fast Multiple Sequence Alignment program.

mutsig2cv analyzes somatic point mutations discovered in DNA sequencing, identifying genes mutated more often than expected by chance.

Napu (Nanopore Analysis Pipeline) is a collection of WDL workflows for variant calling and de novo assembly of ONT data,

The NCBI C++ Toolkit is a set of executables and libraries for a multitude of sequence analysis functions.

The SRA Toolkit and SDK from NCBI is a collection of tools and libraries for using data in the INSDC Sequence Read Archives.

tool for identifying tumorigenic mechanisms and providing diagnostic and therapeutic targets.

NeST-VNN is an interpretable neural network-based model that predicts cell response to a drug. This framework integrates information across multiple levels of cancer cell biology to understand drug response, and can serve to identify and explain biomarkers for clinical application.

nested (now also called TE-greedy) is software to analyze nested LTR transposable elements in DNA sequences, such as reference genomes.

netMHCpan predicts binding of peptides to any MHC molecule of known sequence using artificial neural networks

NetOglyc produces neural network predictions of mucin type GalNAc O-glycosylation sites in mammalian proteins.

Nirvana provides clinical-grade annotation of genomic variants (SNVs, MNVs, insertions, deletions, indels, and SVs (including CNVs). It can be run as a stand-alone package or integrated into larger software tools that require variant annotation.

NucleoATAC is a software for nucleosome calling using ATAC-seq. It can identify the rotational and translational positions of nucleosomes with up to base-pair resolution and provide quantitative measures of nucleosome occupancy.

Octopus is a mapping-based variant caller that implements several calling models within a unified haplotype-aware framework. Octopus takes inspiration from particle filtering by constructing a tree of haplotypes and dynamically pruning and extending the tree based on haplotype posterior probabilities in a sequential manner. This allows octopus to implicitly consider all possible haplotypes at a given loci in reasonable time.

OncodriveCLUSTL is a sequence-based clustering algorithm to detect significant clustering signals across genomic regions. It is based based on a local background model derived from the simulation of mutations accounting for the composition of trior penta-nucleotide context substitutions observed in the cohort under study.

OncodriveFML is a method designed to analyze the pattern of somatic mutations across tumors in both coding and non-coding genomic regions to identify signals of positive selection, and therefore, their involvement in tumorigenesis.

OpenCRAVAT is a new open source, scalable decision support system for variant and gene prioritization. It includses a modular resource catalog to maximize community and developer involvement, and as a result the catalog is being actively developed and growing every month. Resources made available via the store are well-suited for analysis of cancer, as well as Mendelian and complex diseases.

OptiType is a HLA genotyping algorithm based on integer linear programming, capable of producing accurate 4-digit HLA genotyping predictions from NGS data by simultaneously selecting all major and minor HLA Class I alleles.

ORF finder searches for open reading frames (ORFs) in the DNA sequence you enter. The program returns the range of each ORF, along with its protein translation. Use ORF finder to search newly sequenced DNA for potential protein encoding segments, verify predicted protein using newly developed SMART BLAST or regular BLASTP.

OrthoFinder is an accurate and comprehensive platform for comparative genomics. It finds orthogroups and orthologs, infers rooted gene trees for all orthogroups and identifies all of the gene duplication events in those gene trees.

PAINTOR (Probabilistic Annotation INtegraTOR) is a probabilistic framework that integrates association strength with genomic functional annotation data to improve accuracy in selecting plausible causal variants for functional validation.

Pangene constructs a pangenome gene graph in the GFA format from the miniprot alignment between a protein set and multiple genomes..

nf-core/pangenome is a bioinformatics best-practise analysis pipeline for the rendering of a collection of sequences into a pangenome graph.

Phylogenetic Assignment of Named Global Outbreak LINeages. PANGOLIN is a system for identifying phylogenetic COVID lineages that contribute most to active spread.

Parsnp is a command-line-tool for efficient microbial core genome alignment and SNP detection. Parsnp was designed to work in tandem with Gingr, a flexible platform for visualizing genome alignments and phylogenetic trees

PATRIC is an integration of different types of data and software tools that support research on bacterial pathogens. For users that wish command-line access to PATRIC there are the P3-scripts. They are intended to run on locally, going over the network to access the services provided by PATRIC.

PEPPER-Margin-DeepVariant is a haplotype-aware variant calling pipeline for processing third-generation nanopore sequence data. This pipeline is also applicable to PacBio HiFi data, providing an efficient solution with superior performance over the current WhatsHap-DeepVariant standard.

pangenome graph builder.

phasebook is a novel approach for reconstructing the haplotypes of diploid genomes from long reads de novo, that is without the need for a reference genome.

phasing and Allele Specific Expression from RNA-seq. Performs haplotype phasing using read alignments in BAM format from both DNA and RNA based assays, and provides measures of haplotypic expression for RNA based assays.

This Python/C++ code is the accompanying software for the paper PhyloWGS: Reconstructing subclonal composition and evolution from whole-genome sequencing of tumors, with authors Amit G. Deshwar, Shankar Vembu, Christina K. Yung, Gun Ho Jang, Lincoln Stein, and Quaid Morris.

Plotsr generates high-quality visualisation of synteny and structural rearrangements between multiple genomes. For this, it uses the genomic structural annotations between multiple chromosome-level assemblies.

PB-metagenomics-tools is a suite of tools for performing metagenomic analyses using HiFi sequencing data. It comprises four pipelines: HiFi-MAG-Pipeline, Taxonomic-Profiling-Diamond-Megan, Taxonomic-Profiling-Minimap-Megan and Taxonomic-Profiling-Sourmash.

A suite of population scale analysis tools for single-cell genomics data including implementation of Demuxlet/Freemuxlet methods and auxiliary tools

PRANK is a probabilistic multiple alignment program for DNA, codon and amino-acid sequences. PRANK is based on a novel algorithm that treats insertions correctly and avoids over-estimation of the number of deletion events.

predicting library complexity and genome coverage in high-throughput sequencing

Primer3 is a program for designing PCR primers, hybridization probes, and sequencing primers.

Prokka is a software tool for the rapid annotation of prokaryotic genomes.

PSIPRED is a simple and accurate secondary structure prediction method, incorporating two feed-forward neural networks which perform an analysis on output obtained from PSI-BLAST (Position Specific Iterated - BLAST).

Pysamstats is a fast Python and command-line utility for extracting simple statistics against genome positions based on sequence alignments from a SAM or BAM file.

pySCENIC is a lightning-fast python implementation of the SCENIC pipeline (Single-Cell rEgulatory Network Inference and Clustering). It enables biologists to infer transcription factors, gene regulatory networks and cell types from single-cell RNA-seq data.

A tool set for molecular QTL discovery and analysis. It allows to go from the raw sequence data to collection of molecular Quantitative Trait Loci (QTLs) in few easy-to-perform steps.

RandFold computes the probability that, for a given RNA sequence, the Minimum Free Energy (MFE) of the secondary structure is different from a distribution of MFE computed with random sequences.

The BNS (RNA Bind-N-Seq Analysis) pipeline is a set of bioinformatics tools to analyze data from high-throughput sequencing experiments of protein-bound RNAs. The current version includes read splitting, calculation of kmer frequencies and enrichments, QC metrics, production of motif sequence logos, and RNA secondary structure analysis.

regenie is a C++ program for whole genome regression modelling of large genome-wide association studies. It is developed and supported by a team of scientists at the Regeneron Genetics Center. regenie employs the BGEN library.

Modified base model training and application.

RGI (Resistance Gene Identifier) is a robust antimicrobial resistance (AMR) gene predicting tool. It is based on newly curated Comprehensive Antibiotic Research Database (CARD) and allows detection detect AMR genes from thirteen genomes of Pseudomonas strains.

MATS is a computational tool to detect differential alternative splicing events from RNA-Seq data.

RMBlast is a RepeatMasker-compatible version of the standard NCBI blastn program. RMBlast supports RepeatMasker searches by adding a few necessary features to the stock NCBI blastn program.

RNAmmer predicts ribosomal RNA genes in full genome sequences by utilising two levels of Hidden Markov Models: An initial spotter model searches both strands. The spotter model is constructed from highly conserved loci within a structural alignment of known rRNA sequences. Once the spotter model detects an approximate position of a gene, flanking regions are extracted and parsed to the full model which matches the entire gene.

Roary is a high speed stand alone pan genome pipeline, which takes annotated assemblies in GFF3 format (produced by Prokka) and calculates the pan genome.

ROSE (Rank Ordering of Super-Enhancers) is tool for (1) creating stitched enhancers, and (2) separating super-enhancers from typical enhancers. given sequencing data (.bam) and a file of previously identified constituent enhancers (.gff)

R package for large-scale genetic association studies.

Scanpy is a scalable toolkit for analyzing single-cell gene expression data. It includes preprocessing, visualization, clustering, pseudotime and trajectory inference and differential expression testing. The Python-based implementation efficiently deals with datasets of more than one million cells.

scChromHMM provides a suite of tools for rapid processing of single-cell histone modification data to perform chromatin states analysis of the genome within each single-cell. It is an extention of bulk ChromHMM framework, which consumes the HMM model learned from ChromHMM and perform chromatin state analysis by running forward-backward algorithm for each single-cell.

The scDRS application implements an approach that links scRNA-seq with polygenic disease risk at single-cell resolution, independent of annotated cell types. scDRS identifies cells exhibiting excess expression across disease-associated genes implicated by genome-wide association studies (GWASs). Genes whose expression was correlated with the scDRS score across cells (reflecting coexpression with GWAS disease-associated genes) were strongly enriched for gold-standard drug target and Mendelian disease genes.

SCENIC+ is a new method for the inference of enhancer-driven gene regulatory networks (eGRN). It predicts genomic enhancers along with candidate upstream transcription factors (TF) and links these enhancers to candidate target genes. Specific TFs for each cell type or cell state are predicted based on the concordance of TF binding site accessibility, TF expression, and target gene expression.

scVelo is a method to describe the rate of gene expression change for an individual gene at a given time point based on the ratio of its spliced and unspliced messenger RNA (mRNA). It avoids errors in the velocity estimates by solving the full transcriptional dynamics of splicing kinetics using a likelihood-based dynamical model. This generalizes RNA velocity to systems with transient cell states, which are common in development and in response to perturbations.

SEACR is intended to call peaks and enriched regions from sparse Cleavage Under Targets and Release Using Nuclease (CUT&RUN) or chromatin profiling data in which background is dominated by "zeroes" (i.e. regions with no read coverage).

selscan is a tool for haplotype-based scans to detect natural selection, which are useful to identify recent or ongoing positive selection in genomes. It is an efficient multithreaded application that implements Extended Haplotype Homozygosity (EHH), Integrated Haplotype Score (iHS), and Cross-population EHH (XPEHH). selscan accepts phased genotypes in multiple formats, including TPED.

seq2HLA computationally determines human leukocyte antigen (HLA) genotypes of a sample using RNA-Seq sequencing reads.

A JavaScript module for the visualization of genomic sequence graphs. It automatically generates a "tube map"-like visualization of sequence graphs which have been created with vg.

Shimmer is a software package for the characterization of genetic differences between two very similar samples, e.g., a tumor sample and its matched normal tissue sample.

SignalP predicts the presence and location of signal peptide cleavage sites in amino acid sequences from different organisms: Gram-positive bacteria, Gram-negative bacteria, and eukaryotes. The method incorporates a prediction of cleavage sites and a signal peptide/non-signal peptide prediction based on a combination of several artificial neural networks.

SignatureAnalyzer is a tool for the identification of somatic mutational signatures. It employs Bayesian non-negative matrix factorization (NMF).

slivar is a set of command-line tools that enables rapid querying and filtering of VCF files. It facilitates operations on trios and groups and allows arbitrary expressions using simple javascript.

SMC++ is a program for estimating the size history of populations from whole genome sequence data.

smoove simplifies and speeds calling and genotyping SVs for short reads. It also improves specificity by removing many spurious alignment signals that are indicative of low-level noise and often contribute to spurious calls.

SnapATAC is a software package for analyzing scATAC-seq datasets. Conventional assays to map regulatory elements via open chromatin analysis of primary tissues is hindered by sample heterogeneity. SnapATAC dissects cellular heterogeneity in an unbiased manner and map the trajectories of cellular states.

Splice junction recognition model based on a deep residual convolutional neural network to assess splice junctions.

StrainGE is a set of tools to analyse the within-species strain diversity in bacterial populations. It consists of two main components: 1) StrainGST: Strain Genome Search tool, a tool to find close reference genomes for strains present in a sample and 2) StrainGR: Strain Genome Recovery, a tool to perform strain-aware variant calling at low coverages.

Straw is library which allows rapid streaming of contact data from .hic files.

Stripenn is a command line interface python package developed for detection of atchitectural stripes from chromatin conformation capture (3C) data. It implements an algorithm rooted in computer vision for demarcation and quantification of the architectural stripes. Stripenn was demonstrated to outperform existing methods, be applicable in the context of analysis of B and T lymphocytes, and to allow examination of the role of sequence variation on the architectural stripes by studying the conservation of these features in inbred strains of mice.

Comprehensive utilities to explore structural variations in genomes

SynthDNM is a random-forest based classifier that can be readily adapted to new sequencing or variant-calling pipelines by applying a flexible approach to constructing simulated training examples from real data. The optimized SynthDNM classifiers predict de novo SNPs and indels with robust accuracy across multiple methods of variant calling.

tandem-genotypes finds changes in length of tandem repeats, from "long" DNA reads aligned to a genome.

A tool to mask low complexity and short period tandem repeats

Tapestry is a tool to validate and edit small eukaryotic genome assemblies using long sequence reads. It is designed to help identify complete chromosomes, symbionts, haplotypes, complex features and errors in close-to-complete genome assemblies.

TElocal: a tool that utilizes both uniquely and ambiguously mapped reads to quantify transposable element expression at the locus level.

TelomereHunter is a software for the detailed characterization of telomere maintenance mechanism footprints in the genome. The tool is implemented for the analysis of large cancer genome cohorts and provides a variety of diagnostic diagrams as well as machine-readable output for subsequent analysis.

TMHMM predicts transmembrane helices in proteins.

TOGA is a new method that integrates gene annotation, inferring orthologs and classifying genes as intact or lost. TOGA implements a novel machine learning based paradigm to infer orthologous genes between related species and to accurately distinguish orthologs from paralogs or processed pseudogenes. This tutorial explains how to get started using TOGA. It shows how to install and execute TOGA, and how to handle possible issues that may occur.

a suite of tools primarily for the identification of modified nucleotides from nanopore sequencing data.

A tandem repeat in DNA is two or more adjacent, approximate copies of a pattern of nucleotides. Tandem Repeats Finder is a program to locate and display tandem repeats in DNA sequences.

TRGT is a tool for targeted genotyping of tandem repeats from PacBio HiFi data. In addition to the basic size genotyping, TRGT profiles sequence composition, mosaicism, and CpG methylation of each analyzed repeat. TRGT can also create a visualization of reads overlapping the repeats.

trimAl is a tool for the automated removal of spurious sequences or poorly aligned regions from a multiple sequence alignment. It can consider several parameters, alone or in multiple combinations, in order to select the most-reliable positions in the alignment. These include the proportion of sequences with a gap, the level of residue similarity and, if several alignments for the same set of sequences are provided, the consistency level of columns among alignments. Moreover, trimAl is able to manually select a set of columns to be removed from the alignment.

triodenovo is an application for calling de novo mutations in trios for NGS data.

Tcr Receptor Utilities for Solid Tissue (TRUST) is a computational tool to analyze TCR and BCR sequences using unselected RNA sequencing data, profiled from solid tissues, including tumors. TRUST4 performs de novo assembly on V, J, C genes including the hypervariable complementarity-determining region 3 (CDR3) and reports consensus of BCR/TCR sequences. TRUST4 then realigns the contigs to IMGT reference gene sequences to report the corresponding information. TRUST4 supports both single-end and paired-end sequencing data with any read length.

uLTRA implements an alignment method for long RNA sequencing reads based on a novel two-pass collinear chaining algorithm. uLTRA is guided by a database of exon annotations, but it can also be used as a wrapper around minimap2 to align reads outside annotated regions.

Unicycler is an assembly pipeline for bacterial genomes. It can assemble Illumina-only read sets where it functions as a SPAdes-optimiser. It can also assembly long-read-only sets (PacBio or Nanopore) where it runs a miniasm+Racon pipeline.

UROPA is a command line based tool for genomic region annotation

USEARCH is a unique sequence analysis tool with thousands of users world-wide. USEARCH offers search and clustering algorithms that are often orders of magnitude faster than BLAST.

VADR stands for Viral Annotation DefineR. It is a suite of tools for classifying and analyzing sequences homologous to a set of reference models of viral genomes or gene families. It has been mainly tested for analysis of Norovirus, Dengue, and SARS-CoV-2 virus sequences in preparation for submission to the GenBank database.

A platform-independent, technology-independent software tool for identifying SNPs and indels in massively parallel sequencing of individual and pooled samples.

VCF-kit is a collection of utility tools for processing and analyzing the VCF (variant call format) files, including primer generation for variant validation, dendrogram production,genotype imputation from sequence data in linkage studies, and additional tools to be used by statistical and population geneticists.

vConTACT2 is a tool to perform guilt-by-contig-association classification of viral genomic sequence data. It's designed to cluster and provide taxonomic context of metagenomic sequencing data.

Velocyto is a library for the analysis of RNA velocity. It includes a command line tool and an analysis pipeline.

Velvet is a de novo genomic assembler specially designed for short read sequencing technologies, such as Solexa or 454

VICTOR (Variant Interpretation for Clinical Testing Or Research) is a pipeline that can be used for disease gene discovery research or clinical genetic testing

RNA Secondary Structure Prediction and Comparison

VIPER combines the use of several dozen RNA-seq tools, suites, and packages to create a complete pipeline that takes RNA-seq analysis from raw sequencing data all the way through alignment, quality control, unsupervised analyses, differential expression, and downstream pathway analysis

VirSorter2 is a DNA and RNA virus identification tool. It leverages genome-informed database advances across a collection of customized automatic classifiers to improve the accuracy and range of virus sequence detection.

contains seqlogo utility to create sequence logo summarizing sequence alignments

WhatsHap is a software for phasing genomic variants using DNA sequencing reads, also called read-based phasing or haplotype assembly. It is especially suitable for long reads, but works also well with short reads.

A fast xenograft read sorter based on space-efficient k-mer hashing.

xenome is a tool for classifying reads from xenograft source.

The HLA gene complex on human chromosome 6 is one of the most polymorphic regions in the human genome and contributes in large part to the diversity of the immune system. Accurate typing of HLA genes with short-read sequencing data has historically been difficult due to the sequence similarity between the polymorphic alleles. xHLA iteratively refines the mapping results at the amino acid level to achieve high typing accuracy for both class I and II HLA genes.

Xtail is an analysis pipeline tailored for ribosome profiling data that comprehensively and accurately identifies differentially translated genes in pairwise comparisons. Applied on simulated and real datasets, Xtail exhibits high sensitivity with minimal false-positive rates, outperforming existing methods in the accuracy of quantifying differential translations.

xTea (x-Transposable element analyzer), is a tool for identifying TE insertions in whole-genome sequencing data. Whereas existing methods are mostly designed for shortread data, xTea can be applied to both short-read and long-read data. xTea outperforms other short read-based methods for both germline and somatic TE insertion discovery.

This package provides an implementation of the protein structure inference pipeline of AlphaFold v2.0.

This package provides an implementation of the inference pipeline of AlphaFold 3

AlphaLink predicts protein structures using deep learning given a sequence and a set of experimental contacts. It extends OpenFold with crosslinking MS data or other experimental distance restraint by explicitly incorporating them in the OpenFold architecture.

AlphaPulldown is a Python package that streamlines protein-protein interaction screens and high-throughput modelling of higher-order oligomers using AlphaFold-Multimer. It provides a convenient command-line interface, a variety of confidence scores and a graphical analysis tool.

Autodock-GPU performs docking calculations, and processes ligand-receptor poses in parallel over multiple compute units on GPUs.

A pipeline for fully automated prediction of protein coding gene structures with GeneMark-ES/ET and AUGUSTUS in novel eukaryotic genomes.

CCP4 is a suite of programs for protein crystallography and structural biology.

ColabFold batch scripts

Automated particle picker for cryo-EM

The Cambridge Structural Database is the world repository of small molecule crystal structures.

The three-dimensional co-ordinates of each protein are used to calculate residue - residue distance matrices.

This project aims to offer easy access to Deep Learning for segmentation of structures of interest in biomedical 3D scans. It is a system that allows the easy creation of a 3D Convolutional Neural Network, which can be trained to detect and segment structures if corresponding ground truth labels are provided for training. The system processes NIFTI images, making its use straightforward for many biomedical tasks.

DeepMM implements fully automated de novo structure modeling method, MAINMAST, which builds three-dimensional models of a protein from a near-atomic resolution EM map. The method directly traces the protein’s main-chain and identifies Cα positions as tree-graph structures in the EM map.

This application is intended to evaluate coevolutionary structural predictions of fold-switching proteins. State-of-the art algorithms predict that these fold-switching proteins assume only one stable structure. We hypothesize that coevolutionary signatures are being missed. Fold-switching proteins have the ability to transition between two sets of stable secondary and tertiary structure. The approach successfully revealed coevolution of amino acid pairs uniquely corresponding to both conformations of 56 fold-switching proteins from distinct families.

The DSSP program was designed by Wolfgang Kabsch and Chris Sander to standardize secondary structure assignment. DSSP is a database of secondary structure assignments (and much more) for all protein entries in the Protein Data Bank (PDB). DSSP is also the program that calculates DSSP entries from PDB entries. DSSP does not predict secondary structure.

Code and pre-trained weights for Transformer protein language models from the Meta Fundamental AI Research Protein Team (FAIR)

Fast structural similarity search

RoseTTAFold All-Atom (RFAA), a deep network capable of modeling full biological assemblies containing proteins, nucleic acids, small molecules, metals, and covalent modifications given the sequences of the polymers and the atomic bonded geometry of the small molecules and covalent modifications. The heme_binder_diffusion pipeline employs RoseTTAFold All-Atom to perform de novo heme binding protein design.

LAMMPS is a classical molecular dynamics code, and an acronym for Large-scale Atomic/Molecular Massively Parallel Simulator. It runs on a variety of different computer systems, including single processor systems, distributed-memory machines with MPI, and GPU and Xeon Phi systems. LAMMPS is open source software, released under the GNU General Public License.

MDTraj is a python library that allows users to manipulate molecular dynamics (MD) trajectories and perform a variety of analyses, including fast RMSD, solvent accessible surface area, hydrogen bonding, etc.

Mustache is a tool for identifying chromatin loops from HiC and MicroC contact maps

The naccess program calculates the atomic accessible surface defined by rolling a probe of given size around a van der Waals surface.

OmegaFold is the first computational method to successfully predict high-resolution protein structure from a single primary sequence alone. Using a new combination of a protein language model that allows us to make predictions from single sequences and a geometry-inspired transformer model trained on protein structures, OmegaFold outperforms RoseTTAFold and achieves similar prediction accuracy to AlphaFold2 on recently released structures.

PHENIX is a software suite for the automated determination of macromolecular structures using X-ray crystallography and other methods.

A comprehensive molecular visualization product for rendering and animating 3D molecular structures.

PyRosetta is an interactive Python-based interface to the powerful Rosetta molecular modeling suite. It enables users to design their own custom molecular modeling algorithms using Rosetta sampling methods and energy functions.

ReLeaSE (Reinforcement Learning for Structural Evolution) is an application for de-novo Drug Design based on Reinforcement Learning. It integrates two deep neural networks: generative and predictive, that are trained separately but are used jointly to generate novel targeted chemical libraries. ReLeaSE uses simple representation of molecules by their simplified molecular input line entry specification (SMILES) strings only.

Rosetta Fold (RF) dissusion is an open source method for structure generation, with or without conditional information (a motif, target etc). It can perform motif scaffolding, unconditional protein generation, and other tasks.

Accurate prediction of protein structures and interactions using a 3-track network, , in which information at the 1D sequence level, the 2D distance map level, and the 3D coordinate level is successively transformed and integrated.

Accurate prediction of nucleic acid and protein-nucleic acid complexes using RoseTTAFoldNA

A limited number of Schrödinger applications are available on the Biowulf cluster through the Molecular Modeling Interest Group. Most are available through the Maestro GUI.

Sidesplitter reduces over-fitting in both idealised and experimental settings, while maintaining independence between the two sides of a split refinement. It can improve the final resolution in refinements of structures prone to severe over-fitting, such as membrane proteins in detergent micelles.

Subtom is a pipeline for subvolume alignment and averaging of electron cryo-tomography data.

VMD is a molecular visualization program for displaying, animating, and analyzing large biomolecular systems using 3-D graphics and built-in scripting. To use, type vmd at the prompt.

Xplor-NIH is a structure determination program which builds on the X-PLOR v3.851 program, including additional tools developed at the NIH.

ZDOCK predicts protein-docking models, and uses a fast Fourier transform to search all possible binding modes for proteins, evaluating based on shape complementarity, desolvation energy, and electrostatics.

Rapid & standardized annotation of bacterial genomes, MAGs & plasmids

BioNetGen is software for the specification and simulation of rule-based models of biochemical systems, including signal transduction, metabolic, and genetic regulatory networks. BioNetGen is presently a mixture of Perl, C++, and Python.

A publicly available repository of curated receptors, ligands and their interactions. Subunit architecture is included for both ligands and receptors, representing heteromeric complexes accurately.

CytoSig is a data-driven infrastructure hosted by the National Cancer Institute. CytoSig includes both a database of target genes modulated by cytokines and a predictive model of cytokine signaling activity and regulatory cascade from transcriptomic profiles.

dyno is a meta package that installs several other packages from the dynvers (https://github.com/dynverse). It comprises a set of R packages to construct and interpret single-cell trajectories.

eggNOGmapper is a tool for functional annotation of large sets of sequences based on fast orthology assignments using precomputed clusters and phylogenies from the eggNOG database. Orthology assignment is ideally suited for functional inference. However, predicting orthology is computationally intensive at large scale, and most other pipelines are relatively inaccessible (e.g., new assignments only available through database updates), so less precise homology-based functional transfer was previously the default for (meta-)genome annotation.

The GCN_Cancer application employs graph convolutional network (GCN) models to classify the gene expression data samples from The Cancer Genonme Atlas (TCAG) as 33 designated tumor types or as normal. It has been trained on 10,340 cancer samples and 731 normal tissue samples from TCGA dataset.

Gene Set Enrichment Analysis (GSEA) is a computational method that determines whether an a priori defined set of genes shows statistically significant, concordant differences between two biological states (e.g. phenotypes).

MAKER is a portable and easily configurable genome annotation pipeline. Its purpose is to allow smaller eukaryotic and prokaryotic genome projects to independently annotate their genomes and to create genome databases. MAKER identifies repeats, aligns ESTs and proteins to a genome, produces ab-initio gene predictions and automatically synthesizes these data into gene annotations having evidence-based quality values.

NEURON is a simulation environment for modeling individual neurons and networks of neurons. It provides tools for conveniently building, managing, and using models in a way that is numerically sound and computationally efficient. It is particularly well-suited to problems that are closely linked to experimental data, especially those that involve cells with complex anatomical and biophysical properties.

pySCENIC is a lightning-fast python implementation of the SCENIC pipeline (Single-Cell rEgulatory Network Inference and Clustering). It enables biologists to infer transcription factors, gene regulatory networks and cell types from single-cell RNA-seq data.

ROBOT is a command-line tool and library for automating ontology development tasks, with a focus on Open Biological and Biomedical Ontologies (OBO). It can be used as a command-line tool or as a library for any language on the Java Virtual Machine.

SCENIC+ is a new method for the inference of enhancer-driven gene regulatory networks (eGRN). It predicts genomic enhancers along with candidate upstream transcription factors (TF) and links these enhancers to candidate target genes. Specific TFs for each cell type or cell state are predicted based on the concordance of TF binding site accessibility, TF expression, and target gene expression.

A graph neural network model to estimate cell-wise metabolic using single cell RNA-seq data

scVelo is a method to describe the rate of gene expression change for an individual gene at a given time point based on the ratio of its spliced and unspliced messenger RNA (mRNA). It avoids errors in the velocity estimates by solving the full transcriptional dynamics of splicing kinetics using a likelihood-based dynamical model. This generalizes RNA velocity to systems with transient cell states, which are common in development and in response to perturbations.

A comprehensive analysis framework for T-cell and B-cell repertoire sequencing data.

Apptainer allows you to build and run Linux containers with emphasis on use in HPC. Apptainer is the Linux Foundation variant of and successor to the widely popular Singularity.

multiprotocol download utility
Type 'module load aria2'then 'aria2c --help' for more info.

asciinema [as-kee-nuh-muh] is a free and open source solution for recording terminal sessions and sharing them.
Type 'module load asciinema' then 'asciinema' to run.

High-speed fasp-powered file transfers. Mostly used to download data from NCBI, which has an Aspera server. See the data transfer page for details.

Autoconf is an extensible package of M4 macros that produce shell scripts to automatically configure software source code packages.

Command-line tools for Amazon Web Services. Use 'module load python; aws -help' to see the command-line help, or http://aws.amazon.com/cli/.

a command-line utility to copy blobs or files to or from Azure storage

Rapid & standardized annotation of bacterial genomes, MAGs & plasmids

Secure and fast copy utility

The BGEN library contains a reference implementation of the BGEN format, written in C++. The library can be used as the basis for BGEN support in other software.

Circos is a program for the generation of publication-quality, circularly composited renditions of genomic data and related annotations. Circos is particularly suited for visualizing alignments, conservation and intra and inter-chromosomal relationships. Also, Circos is useful to visualize any type of information that benefits from a circular layout. Thus, although it has been designed for the field of genomics, it is sufficiently flexible to be used in other data domains.

The GNU Core Utilities are the basic file, shell and text manipulation utilities of the GNU operating system.

Coherent PDF tools

Rust implementation for the Crypt4GH encryption format.

command line tool and library for transferring data with URLs

Datalad is a tool for uploading and downloading public up-t-to-date neuroimaging datasets.

GNU datamash is a command-line program which performs basic numeric,textual and statistical operations on input textual data files.

DNAnexus is a cloud-based commercial solution for next-generation sequence analysis and visualization. It has a command-line interface (CLI) which can be used to log in to the DNAnexus platform, upload and navigate data, and launch analyses.

Like du but more intuitive. dust provides an easy-to-use visualization tool for disk space usage allow you to easily find the subfolders eating up all of your quota.

CLI tool to manage the download of large quantities of files from DNAnexus

Entrez Direct (EDirect) is an advanced method for accessing the NCBI's set of interconnected databases (publication, sequence, structure, gene, variation, expression, etc.) from a UNIX terminal window.

Eye of Gnome image viewer

A modern, maintained replacement for ls. A replacement for the now-unmaintained exa.

fd is a simple, fast and user-friendly alternative to find. While it does not aim to support all of find's powerful functionality, it provides sensible (opinionated) defaults for a majority of use cases.

🌸 A command-line fuzzy finder

The GDC Data Transfer Tool provides an optimized method of transferring data to and from the GDC, and enables resumption of interrupted transfers.

gdrive is a command line utility for interacting with Google Drive.

Ghostscript is an interpreter for the PostScript language and for PDF.

Command-line client for setting up and monitoring recurring or scheduled Globus transfers

The GLPK (GNU Linear Programming Kit) package is intended for solving large-scale linear programming (LP), mixed integer programming (MIP), and other related problems. It is a set of routines written in ANSI C and organized in the form of a callable library.

Gnuplot is a portable command-line driven graphing utility to visualize mathematical functions and data interactively, and can support many non-interactive uses such as web scripting.
Type 'gnuplot' to run, or 'module avail gnuplot' to see other available versions.

Google Cloud SDK is a set of tools that you can use to manage resources and applications hosted on Google Cloud Platform. These include the gcloud, gsutil, and bq command line tools. See docs at https://cloud.google.com/sdk/docs/how-to.
Type 'module load google-cloud-sdk' to use on Biowulf.

Grace is a WYSIWYG 2D plotting tool for the X-Window system. It is a successor to Xmgr.
Type 'module load grace', then 'xmgrace' or 'gracebat' to run.

Graphviz is open source graph visualization software. Graph visualization is a way of representing structural information as diagrams of abstract graphs and networks. It has important applications in networking, bioinformatics, software engineering, database and web design, machine learning, and in visual interfaces for other technical domains.

Groff (GNU troff) is a typesetting system that reads plain text mixed with formatting commands and produces formatted output. Output may be PostScript or PDF, html, or ASCII/UTF8 for display at the terminal. Formatting commands may be either low-level typesetting requests (“primitives”) or macros from a supplied set. Users may also write their own macros. All three may be combined.

gThumb is an image viewer and browser

h5utils is a set of utilities for visualization and conversion of scientific data in the free, portable HDF5 format. Type 'module load h5utils' to access the executables (e.g. h5topng)

A modern alternative to the watch command, records the differences in execution results and can check these differences later.

HyperQueue (HQ) lets you build a computation plan consisting of a large amount of tasks and then execute it transparently over a system like SLURM/PBS. It dynamically groups tasks into SLURM/PBS jobs and distributes them to fully utilize allocated nodes.

ImageMagick is a software suite to create, edit, compose, or convert bitmap images. It can read and write images in a variety of formats.
The ImageMagick tools are available by default (e.g. type 'convert') or use 'module load ImageMagick' to load the latest versions.

Inkscape is a Free and open source vector graphics editor

JAX is a python library that brings Autograd and XLA (Accelerated Linear Algebra) together for high-performance machine learning research. JAX uses XLA to compile and run your NumPy programs on GPUs. Compilation happens under the hood by default, with library calls getting just-in-time compiled and executed. But JAX also lets you just-in-time compile your own Python functions into XLA-optimized kernels using a one-function API, jit

A small utility to create JSON objects from command line arguments.

Command line json processor

Jupyter is an open-source web application that allows you to create and share documents that contain live code, equations, visualizations and explanatory text. Uses include: data cleaning and transformation, numerical simulation, statistical modeling, machine learning and much more.

javascript shell based on v8, but with fasta/fastq parser and the ability to open files on the filesystem

Krona allows hierarchical data to be explored with zooming, multi-layered pie charts. Krona charts can be created using an Excel template or KronaTools, which includes support for several bioinformatics tools and raw data formats. The interactive charts are self-contained and can be viewed with any modern web browser.

lbzip2 is a multi-threaded compression utility with support for bzip2 compressed file format.

MariaDB Server is one of the most popular database servers in the world. It’s made by the original developers of MySQL and guaranteed to stay open source.

GNU Midnight Commander is a visual file manager, with a feature rich full-screen text mode application that allows you to copy, move and delete files and whole directory trees, search for files and run commands in the subshell. Type module load mc and then the command mc to get started.

Minimal Eclipe GIT. A graphical front end for git based off the eclipse IDE.

MySQL is an open-source relational database management system.

ncdu is an interactive disk usage analyzer with a text UI designed to find space-hogging directories.

In order to submit data to the National Institute of Mental Health Data Archives (NDA), users must validate their data to ensure it complies with the required format. This is done using the NDA validation tool, vtcmd. Additionally, users can package and download data from NDA as well, using the downloadcmd tool.

Netpbm is a toolkit for manipulation of graphic images, including conversion of images between a variety of different formats. There are over 300 separate tools in the package including converters for about 100 graphics formats. Examples of the sort of image manipulation we're talking about are: Shrinking an image by 10%; Cutting the top half off of an image; Making a mirror image; Creating a sequence of images that fade from one image to another.

nvchecker (short for new version checker) is for checking if a new version of some software has been released.

NVTOP stands for Neat Videocard TOP, a (h)top like task monitor for GPUs and accelerators. It can handle multiple GPUs and print information about them in a htop-familiar way.

OMERO.py provides Python bindings to the OMERO.blitz server as well as a pluggable command-line interface. ezomero is module with convenience functions for writing Python code that interacts with OMERO.

This command-line tool allows you to upload and download OpenNeuro.org datasets without a browser.

GNU parallel is a shell tool for executing jobs in parallel using one or more computers.

ParaView is an open-source, multi-platform data analysis and visualization application.

patchelf is a small utility to modify the dynamic linker and RPATH of ELF executables.

PATRIC is an integration of different types of data and software tools that support research on bacterial pathogens. For users that wish command-line access to PATRIC there are the P3-scripts. They are intended to run on locally, going over the network to access the services provided by PATRIC.

A simple PDF to SVG converter using the Poppler and Cairo libraries.

pigz (parallel implementation of gzip) is a fully functional replacement for gzip that exploits multiple processors and multiple cores to the hilt when compressing data.

Pilon is a software tool which can be used to: 1) Automatically improve draft assemblies 2) Find variation among strains, including large event detection

Portal_client is Python-based client for downloading data made available through portals powered by the GDC-based portal system. There are several portals running on the internet to support various research efforts. Notably, the Neuroscience Multi-omic Archive (NeMO, https://nemoarchive.org/) and the Human Microbiome Project Data Analysis and Coordination Center (hmpdacc.org) use the portal to enable data exploration and download. The client accepts a manifest file as an input. This file contains URLs to the files to be downloaded. Manifest files can be generated using the shopping cart functionality of the portal's query interface.

PostgreSQL is a powerful, open source object-relational database system with over 35 years of active development that has earned it a strong reputation for reliability, feature robustness, and performance.

POVRAY (Persistence of Vision RAYtracer) is a high-quality tool for creating three-dimensional graphics. Raytraced images are publication-quality and 'photo-realistic', but are computationally expensive so that large images can take many hours to create.

A framework for managing and maintaining multi-language pre-commit hooks

A download client for the European Genome-phenome Archive (EGA). The EGA is designed to be a repository for all types of sequence and genotype experiments, including case-control, population, and family studies.

QPDF is a command-line tool and C++ library that performs content-preserving transformations on PDF files

Quarto is an open-source scientific and technical publishing system built on Pandoc

Rclone is a utility for synchronizing directories on a file-based storage system (e.g. /home or /data) with an object store such as Amazon S3. It uses the S3 protocol, and it can be used with the HPC object storage system.

rdfind is a program that finds duplicate files. It is useful for compressing backup directories or just finding duplicate files. It compares files based on their content, NOT on their file names. After typing module load rdfind, type man rdfind for more information.

Local Ancestry and Admixture Inference

ripgrep, a modern line-oriented search tool providing rg

RStudio Server is a web-based R IDE similar to RStudio Desktop.

Use the Seven Bridges Command Line Interface (SB CLI) to programmatically access and automate your interaction with the Platform via the API. The CLI is called by a simple command: sb.

Screen is a full-screen window manager that multiplexes a physical terminal between several processes, typically interactive shells.

Singularity is a container platform focused on supporting ``Mobility of Compute``. It allows users to emulate, and share custom Linux environments allowing for the creation of self-contained development stacks.

Rapidly extracts SNPs from a multi-FASTA alignment

Apache Spark is a fast and general engine for large-scale data processing. It is commonly used as an in-memory alternative to Hadoop MapReduce.

Spielbash helps you automate asciicasts with asciinema.

SQLite is a software library that implements a self-contained, serverless, zero-configuration, transactional SQL database engine.

GNU Stow is a symlink farm manager which takes distinct packages of software and/or data located in separate directories on the filesystem, and makes them appear to be installed in the same place.

Swarm is a script designed to simplify submitting a group of commands to the Biowulf cluster.

SWIG is a software development tool that connects programs written in C and C++ with a variety of high-level programming languages.

The synapseclient package provides an interface to Synapse, a collaborative workspace for reproducible, data intensive research projects

tmux is a terminal multiplexer.
Type 'module load tmux' to load the module, then 'tmux --help'

UK Biobank tools for downloading and processing UKBB data. Type 'module load ukbb' to access them.

An extremely fast Python package and project manager, written in Rust.

VarTrix is a software tool for extracting single cell variant information from 10x Genomics single cell data.

vcf2db creates a gemini-compatible database from a VCF.

VisiData is an interactive multitool for tabular data

Interactive cli tool for HTTP inspection

eXtensible ARchiver

xclip is a command line utility that is designed to run on any system with an X11 implementation. It provides an interface to X selections ("the clipboard") from the command line. It can read data from standard in or a file and place it in an X selection for pasting into other X applications. xclip can also print an X selection to standard out, which can then be redirected to a file or another program.

Xpdf is a free PDF viewer and toolkit, including a text extractor, image converter, HTML converter, and more. Most of the tools are available as open source.

Zstandard, or zstd as short version, is a fast lossless compression algorithm, targeting real-time compression scenarios at zlib-level and better compression ratios.

A Python-based wrapper for the Cromwell pipeline system used by ENCODE pipelines. Includes the croo and qc2tsv commands for parsing Cromwell and quality control output.

A Workflow Management System geared towards scientific workflows.

A Workflow Description Language local runner. An alternative to Cromwell for running WDL pipelines

Data-driven computational pipelines

Snakemake aims to reduce the complexity of creating workflows by providing a fast and comfortable execution environment, together with a clean and modern domain specific specification language (DSL) in python style. It is well suited for bioinformatic workflows.