Find open-source science resources

A small <720Kb C++ windows utility. That allows you to load Ancestry, 23andMe, FTDNA, or Genes for Good RAW DNA files search them, merge them. covert them to Ancestry format. But also create files from peer reviewed publications to compare with you loaded data to give your genetic disposition for the condition you have entered the data for an statistical risk if OR values are included. Included with the program are example files for Type 2 Diabetes risk factors. (As I have type 2 Diabetes so I could test the results).

The complexity of high-throughput quantitative omics experiments often leads to low replicates numbers and many missing values. We implemented a new test to simultaneously consider missing values and quantitative changes, which we combined with well-performing statistical tests for high confidence detection of differentially regulated features. The package contains functions to run the test and to visualize the results.

The application of 2nd and 3rd generation High Throughput Sequencing (HTS) technologies has deeply reshaped experimental method to investigate microbial communities and obtain a taxonomic and functional profile of the invetigated community. Shotgun Metagenomics allow to quickly obtain a representation of microorganisms genomes characterizing a particular environment. In order to obtain a fast e reliable taxonomic classification of microorganisms genomes we present kMetaShot, an alignment-free taxonomic classifier based on k-mer/minimizer counting.

Trackastra is a transformer-based cell-tracking tool for live-cell microscopy. It links already segmented cell instances across time by predicting associations between detections. It supports greedy tracking with or without cell divisions, optional ILP-based linking, pretrained tracking models, and export of tracked masks and lineage information in Cell Tracking Challenge format.

toscca is an R package to perform Thresholded Ordered Sparse Canonical Correlation Analysis (TOSCCA).

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

Generate HTML or PDF reports to explore a set of regions such as the results from annotation-agnostic expression analysis of RNA-seq data at base-pair resolution performed by derfinder. You can also create reports for DESeq2 or edgeR results.

This package provides functions for annotation-agnostic differential expression analysis of RNA-seq data. Two implementations of the DER Finder approach are included in this package: (1) single base-level F-statistics and (2) DER identification at the expressed regions-level. The DER Finder approach can also be used to identify differentially bounded ChIP-seq peaks.

Helper package for speeding up the derfinder package when using multiple cores. This package is particularly useful when using BiocParallel and it helps reduce the time spent loading the full derfinder package when running the F-statistics calculation in parallel.

maeparser is a parser for Schrodinger Maestro files.

This tutorial aims to illustrate the process of analyzing a membrane molecular dynamics (MD) simulation, step by step, using the BioExcel Building Blocks (biobb)

This tutorial aims to illustrate the process of checking a molecular structure before using it as an input for a Molecular Dynamics simulation, step by step, using the BioExcel Building Blocks (biobb).

This BioExcel Building Blocks library (BioBB) workflow provides a pipeline to setup DNA structures for the Ascona B-DNA Consortium (ABC) members. It follows the work started with the NAFlex tool to offer a single, reproducible pipeline for structure preparation, ensuring reproducibility and coherence between all the members of the consortium.

This tutorial aims to illustrate the process of generating protein conformational ensembles from 3D structures using Coarse-Grained tools from the FlexServ server and analysing its molecular flexibility

This tutorial aims to illustrate the process of computing a conformational transition between two known structural conformations of a protein, step by step, using the BioExcel Building Blocks (biobb).

This tutorial aims to illustrate the process of protein-ligand docking, step by step, using the BioExcel Building Blocks library (biobb).

This tutorial aims to illustrate the process of computing classical molecular interaction potentials from protein structures step by step, using the BioExcel Building Blocks library (biobb)

This tutorial aims to illustrate how to compute a fast-growth mutation free energy calculation, step by step, using the BioExcel Building Blocks library (biobb). The particular example used is the Staphylococcal nuclease protein (PDB code 1STN), a small, minimal protein, appropriate for a short tutorial.

RFantibody is a pipeline for structure-based de novo antibody and nanobody design, integrating backbone design with RFdiffusion, sequence design with ProteinMPNN, and structure prediction with RoseTTAFold2. It provides a comprehensive toolset for generating and filtering high-quality antibody designs.

Logistic Factor Analysis is a method for a PCA analogue on Binomial data via estimation of latent structure in the natural parameter. The main method estimates genetic population structure from genotype data. There are also methods for estimating individual-specific allele frequencies using the population structure. Lastly, a structured Hardy-Weinberg equilibrium (HWE) test is developed, which quantifies the goodness of fit of the genotype data to the estimated population structure, via the estimated individual-specific allele frequencies (all of which generalizes traditional HWE tests).

GCAT is an association test for genome wide association studies that controls for population structure under a general class of trait models. This test conditions on the trait, which makes it immune to confounding by unmodeled environmental factors. Population structure is modeled via logistic factors, which are estimated using the `lfa` package.

SCENIC+ is a python package to build gene regulatory networks (GRNs) using combined or separate single-cell gene expression (scRNA-seq) and single-cell chromatin accessibility (scATAC-seq) data.

CebraEM is a bioinformatics tool for analyzing and processing large-scale imaging data, providing a pipeline for segmentation, annotation, and analysis with support for both Linux and Windows environments. It includes modules for core functionality, annotation, and network analysis, requiring specific dependencies and a conda environment for execution.

Pairwise SNP distance matrix from a FASTA sequence alignment

This package provides functionality for interactive visualization of RNA-seq datasets based on Principal Components Analysis. The methods provided allow for quick information extraction and effective data exploration. A Shiny application encapsulates the whole analysis.

This package provides functions for an Interactive Differential Expression AnaLysis of RNA-sequencing datasets, to extract quickly and effectively information downstream the step of differential expression. A Shiny application encapsulates the whole package. Support for reproducibility of the whole analysis is provided by means of a template report which gets automatically compiled and can be stored/shared.

Use this package to interface with the WikiPathways API. It provides programmatic access to WikiPathways content in multiple data and image formats, including official monthly release files and convenient GMT read/write functions.

S3segmenter is a Matlab-based set of functions that generates single cell (nuclei and cytoplasm) label masks.

The aim of TCGAbiolinks is : i) facilitate the GDC open-access data retrieval, ii) prepare the data using the appropriate pre-processing strategies, iii) provide the means to carry out different standard analyses and iv) to easily reproduce earlier research results. In more detail, the package provides multiple methods for analysis (e.g., differential expression analysis, identifying differentially methylated regions) and methods for visualization (e.g., survival plots, volcano plots, starburst plots) in order to easily develop complete analysis pipelines.

Tandem repeat genotyping with long reads, being a modified version of HipSTR.

Permutation analysis, based on Monte Carlo sampling, for testing the hypothesis that the number of conserved differentially methylated elements, between several generations, is associated to an effect inherited from a treatment and that stochastic effect can be dismissed.

MitoFinder: efficient automated large-scale extraction of mitogenomic data from high throughput sequencing data

A rule-based checker for the bio.tools database.

A simple toolbox of two ImageJ plugins for quantifying adipocyte morphology and function in tissues and in vitro.

Minigraph is a sequence-to-graph mapper and graph constructor. For graph generation, it aligns a query sequence against a sequence graph and incrementally augments an existing graph with long query subsequences diverged from the graph.

ValidSense is a toolbox for assessing agreement between two quantitative methods or devices measuring the same quantity using the Limits of Agreement (LoA) analysis, also known as the Bland-Altman analysis.

This tutorial aims to illustrate the process of setting up a simulation system containing a protein, step by step, using the BioExcel Building Blocks (biobb) REST API. The particular example used is the Lysozyme protein (PDB code 1AKI).

A database system designed to store, organize, and manage large-scale nucleotide sequencing read data (like PacBio reads) for the Dazzler genome assembler

Git repo for Bio::DB::HTS module on CPAN, providing Perl links into HTSlib

Cis-mQTL mapping protocol for placental methylome using TensorQTL, providing a step-by-step guide for genotype QC, imputation, and analysis with detailed scripts and commands.

Utility that performs integrated analyses of 'gene' data (a set of genes or other genomic features) with 'peak' data (a set of regions, for example ChIP peaks) to identify the genes nearest to each peak, and vice versa.

Pangolin is a deep-learning based method for predicting splice site strengths (for details, see Zeng and Li, Genome Biology 2022). It is available as a command-line tool that can be run on a VCF or CSV file containing variants of interest; Pangolin will predict changes in splice site strength due to each variant, and return a file of the same format. Pangolin's models can also be used with custom sequences.

SMILES-Mapper is a small web app that allows students to easily visualize how SMILES & InChI strings are created together with other mol-block file formats such as .mol and .sdf files.

Short Python script (using Biopython library functions) to extract sequences from a FASTA, QUAL, FASTQ, or SFF file based on the list of IDs given by a column of a tabular file. The output order follows that of the tabular file, and if there are duplicates in the tabular file, there will be duplicates in the output sequence file.

In silico derivatization for GC. The GC-derivatization tool converts carbonyl groups to C═N-OCH3 (MeOX) and transforms acidic protons into -Si(CH3)3 (TMS). Key functionalities include checking for specific groups, removing derivatization groups, and adding derivatization groups to molecules.

adapter trimmer for Oxford Nanopore reads

Scalable gVCF merging and joint variant calling for population sequencing projects

Circlator is a tool to circularize genome assemblies. It will attempt to identify each circular sequence and output a linearised version of it. It does this by assembling all reads that map to contig ends and comparing the resulting contigs with the input assembly.

Content-Aware Image Restoration for Cryo-Transmission Electron Microscopy Data

AlphaPickle is a Python tool that converts AlphaFold and ColabFold output files into user-friendly CSV files and plots, enabling easy analysis and visualization of protein prediction data without requiring programming expertise. It processes .pkl, .json, and PDB files to extract and visualize metrics like pLDDT and PAE.