BOL: Related items

Ancestral sequence reconstruction (ASR) or ancestral gene/sequence reconstruction/resurrection tools to study molecular evolution

Jit — Tue, 30 May 2017 04:20:05 -0500

Ancestral sequence reconstruction (ASR) – also known as ancestral gene/sequence reconstruction/resurrection – is a technique used in the study of molecular evolution. The method consists of the synthesis of an ancestral gene and expression of the corresponding ancestral protein. The idea of protein 'resurrection' was suggested in 1963 by Pauling and Zuckerkandl. Some early efforts were made in the eighties-nineties, led by the laboratory of Steven A. Benner, showing the potential of this technique – one that only started to be fulfilled in the post-genomic era. Thanks to the improvement of algorithms and of better sequencing and synthesis techniques, the method was developed further in the early 2000s to allow the resurrection of a greater variety of and much more ancient genes. Over the last decade, ancestral protein resurrection has developed as a strategy to reveal the mechanisms and dynamics of protein evolution.

Following are the list of Ancestral /sequence/ reconstruction (ASR) tools:

inferCars

Reconstructs contiguous regions of an ancestral genome. Given information about adjacencies between conserved segments in each modern species, our goal is to infer segment order in the ancestral genome. To get a clean and precise statement of the problem, we formalize it using graph theory. We develop an algorithm that identifies a most parsimonious scenario for the history of each individual adjacency, although the whole-genome prediction is not guaranteed to optimize traditional measures like the number of breakpoints. We introduce weights to the graph edges to model the reliability of each adjacency.

ANGES:reconstructing ANcestral GEnomeS maps

A suite of Python programs that allows reconstructing ancestral genome maps from the comparison of the organization of extant-related genomes. ANGES can reconstruct ancestral genome maps for multichromosomal linear genomes and unichromosomal circular genomes. It implements methods inspired from techniques developed to compute physical maps of extant genomes.

Cocos

Constructs phylogenies of multi-domain proteins. With a given species tree and domain phylogenies, the procedure infers the composition of ancestral multi-domain proteins. Cocos implements and extend a suggested algorithmic approach by Behzadi and Vingron in an easy-to-use program. Such method could be applied to reconstruction of partial homologous units such as bacterial operons or protein complexes.

MySSP

Constructs an initial DNA sequence at the root of the tree and simulates evolution across the tree using a variety of common models of DNA evolution. MySSP is a program for the simulation of DNA sequence evolution across a phylogenetic tree. It is designed for large-scale studies, including simulation of multiple replicates and outputs sequences into NEXUS, MEGA, or FASTA formats. MySSP has a fairly simple graphical user interface (GUI) for basic use, but also has a specialized batch script interpreter to allow for more complicated or large-scale simulations.

PARANA: Parsimonious Ancestral Reconstruction And Network Analysis

Performs parsimony based inference of ancestral biological networks. Given multiple extant networks and phylogenetic information relating extant nodes, PARANA finds a parsimonious set of ancestral interaction events (edge gains and losses) which explain the extant networks. The framework adopted by PARANA is able to represent network evolution under models that support gene duplication and loss and independent interaction gain and loss. The method works on both directed and undirected networks and can incorporate asymmetric interaction gain and loss costs. In contrast to previous approaches, PARANA does not require knowing the relative ordering of unrelated duplication events and thus, works on phylogenetic trees even where branch lengths are not provided.

GapAdj: Gapped Adjacencies

A synteny-based method that is flexible enough to handle a model of evolution involving whole genome duplication events, in addition to rearrangements, gene insertions, and losses. Ancestral relationships between markers are defined in term of Gapped Adjacencies, i.e. pairs of markers separated by up to a given number of markers. It improves on a previous restricted to direct adjacencies, which revealed a high accuracy for adjacency prediction, but with the drawback of being overly conservative, i.e. of generating a large number of contiguous ancestral regions (CARs).

ANCESTOR

A web server allowing one to easily and quickly perform the last three steps of the ancestral genome reconstruction procedure. Ancestors implements several alignment algorithms, an indel maximum likelihood solver and a context-dependent maximum likelihood substitution inference algorithm. The results presented by the server include the posterior probabilities for the last two steps of the ancestral genome reconstruction and the expected error rate of each ancestral base prediction.

ProCARs

Reconstructs ancestral gene orders as contiguous ancestral regions (CARs) with a progressive homology-based method. ProCARs runs from a phylogeny tree (without branch lengths needed) with a marked ancestor and a block file. This homology-based method is based on iteratively detecting and assembling ancestral adjacencies, while allowing some micro-rearrangements of synteny blocks at the extremities of the progressively assembled CARs. The method starts with a set of blocks as the initial set of CARs, and detects iteratively the potential ancestral adjacencies between extremities of CARs, while building up the CARs progressively by adding, at each step, new non-conflicting adjacencies that induce the less homoplasy phenomenon. The species tree is used, in some additional internal steps, to compute a score for the remaining conflicting adjacencies, and to detect other reliable adjacencies, in order to reach completely assembled ancestral genomes.

FastML

A user-friendly tool for the reconstruction of ancestral sequences. FastML implements various novel features that differentiate it from existing tools: (i) FastML uses an indel-coding method, in which each gap, possibly spanning multiples sites, is coded as binary data. FastML then reconstructs ancestral indel states assuming a continuous time Markov process. FastML provides the most likely ancestral sequences, integrating both indels and characters; (ii) FastML accounts for uncertainty in ancestral states: it provides not only the posterior probabilities for each character and indel at each sequence position, but also a sample of ancestral sequences from this posterior distribution, and a list of the k-most likely ancestral sequences; (iii) FastML implements a large array of evolutionary models, which makes it generic and applicable for nucleotide, protein and codon sequences; and (iv) a graphical representation of the results is provided, including, for example, a graphical logo of the inferred ancestral sequences.

maxAlike

Reconstructs a genomic sequence for a specific taxon based on sequence homologs in other species. The input is a multiple sequence alignment and a phylogenetic tree that also contains the target species. For this target species, the algorithm computes nucleotide probabilities at each sequence position. Consensus sequences are then reconstructed based on a certain confidence level.

MLGO: Maximum Likelihood for Gene Order Analysis

A web tool for the reconstruction of phylogeny and/or ancestral genomes from gene-order data. MLGO was designed for analysis of large-scale genomic changes including not only rearrangements but also gene insertions, deletions and duplications. MLGO can be used to infer a phylogeny from genome rearrangement and gene order data, and can also obtain an estimation of ancestral genomes, given an input tree. MLGO takes the advantage of binary encoding on gene-order data, supports a fairly general model of genomic evolution (rearrangements plus duplications, insertions, and losses of genomic regions), and successfully accommodates itself into the framework of maximized likelihood.

Image Reference : Wiki

Computer simulation of genetic mechanism !!

Jit — Sun, 13 Mar 2016 09:29:56 -0500

Computer simulation is the discipline of designing a model of an actual or theoretical physical/biological system, executing the model on a digital computer, and analyzing the execution output. Simulation embodies the principle of ``learning by doing'' --- to learn about the system we must first build a model of some sort and then operate the model. The use of simulation is an activity that is as natural as a child who role plays. Children understand the world around them by simulating (with toys and figurines) most of their interactions with other people, animals and objects. As adults, we lose some of this childlike behavior but recapture it later on through computer simulation. To understand reality and all of its complexity, we must build artificial objects and dynamically act out roles with them. Computer simulation is the electronic equivalent of this type of role playing and it serves to drive synthetic environments and virtual worlds. Within the overall task of simulation, there are three primary sub-fields: model design, model execution and model analysis

Simulation models have become important tools in Bioinformatics studies. There are many reasons for this, but we emphasize three of the more important:

(1) they enable exploration of hypotheses, and as such, have become invaluable means to guide research;

(2) they are unique approaches to integrate (in the literal term of the word) biological knowledge, in the form of experimental results; and

(3) they enable connecting biology with other fields of study ranging from physiology to genomics;

This blog, and this software list, is intended to guide the potential user of simulation models.
It is not, in any way, meant to be comprehensive on the very diverse simulation tools that already exist, but focuses on mechanistic, dynamic models. Similarly, it is not meant to provide any coverage of the breadth of applications; however, for interested readers, we provide references to use as a possible starting point.

Simulation models are meant to answer questions which scientists have in a dynamic, quantitative, and often, a pictorial way. Much of the bioinformatics research and its applications, in particular, involve a large number of components, actors, and factors. Assembling these in a coherent framework may seem a daunting task, especially for beginners, and can lead to confusion, even for experienced scientists, especially if the objectives of such an exercise are not well defined. Followings are the list of tools bioinformatician may use to analyze and provide answers to complex biological mechanisms and related problems.

Software Resource	Brief Description and Homepage
Aladyn	Tools to investigate how demographic parameters, populations genetics and abiotic conditions affect the rate of adaptation http://www.katja-schiffers.eu/research.html
ALF	A Simulation Framework for Genome Evolution http://www.cbrg.ethz.ch/alf
ART	ART is a set of simulation tools to generate synthetic next-generation sequencing data by mimicking real sequencing process with empirical error models or quality profiles. http://www.niehs.nih.gov/research/resources/software/biostatistics/art/
BAMSurgeon	Methods for realistic simulation of mutations in real data. https://github.com/adamewing/bamsurgeon
Bayesian Serial SimCoal	Bayesian Serial SimCoal, (BayeSSC) is a modification of SIMCOAL 1.0, a program written by Laurent Excoffier, John Novembre, and Stefan Schneider. http://www.stanford.edu/group/hadlylab/ssc/index.html
BaySICS	An integral platform with a graphical interface for statistical inference based on approximate Bayesian computation. https://sites.google.com/site/baysicsabc/
BEERS	BEERS was designed to benchmark RNA-Seq alignment algorithms and also algorithms that aim to reconstruct different isoforms and alternate splicing from RNA-Seq data http://cbil.upenn.edu/beers/
BOTTLENECK	Bottleneck is a program for detecting recent effective population size reductions from allele data frequencies http://www.ensam.inra.fr/urlb/bottleneck/bottleneck.html
BottleSim	BottleSim is a computer simulation program for simulating the process of population bottlenecks http://chkuo.name/software/bottlesim.html
CASS	Protein Sequence Simulation https://liberles.cst.temple.edu/software/cass/index.html
CDPOP	CDPOP is a landscape genetics tool for simulating the emergence of spatial genetic structure in populations resulting from specified landscape processes governing organism movement behavior. http://cel.dbs.umt.edu/cdpop
Classical Genetics Simulator	Web-based simulation software http://www.cgslab.com/
CoaSim	CoaSim is a tool for simulating the coalescent process with recombination and geneconversion under various demographic models. http://users-birc.au.dk/mailund/coasim/index.html
cosi	The cosi package is written in C and is available as a tar file. http://www.broadinstitute.org/~sfs/cosi/
CS-PSeq-Gen	A program to simulate the evolution of protein sequences under the constraints of the information of a particular reconstructed phylogeny http://bioserv.rpbs.univ-paris-diderot.fr/software/cs-pseq-gen/
DAWG	An application designed to simulate the evolution of recombinant DNA sequences in continuous time http://scit.us/projects/dawg
Easypop	EASYPOP is an individual based model intended to simulate datasets under a very broad range of conditions http://www.unil.ch/dee/en/home/menuinst/softwares--dataset/softwares/easypop.html
EggLib	EggLib is a C++/Python library and program package for evolutionary genetics and genomics. http://egglib.sourceforge.net/
EpiSIM	EpiSIM: simulation of multiple epistasis, linkage disequilibrium patterns and haplotype blocks for genome-wide interaction analysis https://sourceforge.net/projects/episimsimulator/files/
EvolSimulator	A simulation test bed for hypotheses of genome evolution http://acb.qfab.org/acb/evolsim/
EvolveAGene	A realistic coding sequence simulation program that separates mutation from selection and allows the user to set selection conditions http://bellinghamresearchinstitute.com/software/index.html
fastsimcoal	A continuous-‐time coalescent simulator of genomic diversity under arbitrarily complex evolutionary scenarios http://cmpg.unibe.ch/software/fastsimcoal/
FastSLINK	Simulation of Marker and Phenotype Data in Pedigrees https://watson.hgen.pitt.edu/
FFPopSim	C++/Python library for population genetics. http://webdav.tuebingen.mpg.de/ffpopsim/
FLUX SIMULATOR	The Flux Simulator aims at providing a deterministic in silico reproduction of the experimental pipelines for RNA-Seq, employing a minimal set of parameters. http://sammeth.net/confluence/display/sim/home
forqs	Forward-in-time simulation of Recombination, Quantitative Traits, and Selection https://bitbucket.org/dkessner/forqs
ForSim	ForSim: A Forward Evolutionary Computer Simulation http://anth.la.psu.edu/research/weiss-lab/research/research
ForwSim	The program given below is based on the algorithm described in Padhukasahasram et al. 2008 to simulate genetic drift in a standard Wright-Fisher process. http://badri-populationgeneticsimulators.blogspot.com/
FPG	Forward Population Genetic simulation https://bio.cst.temple.edu/~hey/software/software.htm#fpg
FREGENE	FREGENE is a C++ program that simulates sequence-like data over large genomic regions in large diploid populations. http://www.ebi.ac.uk/projects/bargen
FIGG	FIGG is a genome simulation tool that uses known or theorized variation frequency, per a given fragment size and grouped by GC content across a genome to model new genomes in FASTA format while tracking applied mutations for use in analysis http://insilicogenome.sourceforge.net/
fwdpp	A C++ template library for implementing efficient forward simulations. http://molpopgen.github.io/fwdpp/
GAMETES	Genetic Architecture Model Emulator for Testing and Evaluating Software: Simulates complex SNP models with pure, strict epistatic interactions with n-loci. http://sourceforge.net/projects/gametes/?source=navbar
GASP	Genometric Analysis Simulation Program. A software tool for testing and investigating methods in statistical genetics by generating samples of family data based on user specified models. http://research.nhgri.nih.gov/gasp/
GCTA	Genome-wide Complex Trait Analysis http://www.complextraitgenomics.com/software/gcta/download.html
GemSIM	Next generation sequencing read simulator http://sourceforge.net/projects/gemsim/
GeneArtisan	Simulation of Markers in Case-Control Study Designs http://www.rannala.org/?page_id=241
GENOME	A rapid coalescent-based whole genome simulator http://www.sph.umich.edu/csg/liang/genome/
GenomePop2	GenomePop2 is a specialization of the program GenomePop just to manage SNPs under more flexible and useful settings. If you need models with more than 2 alleles please use the GenomePop program version. https://ritchielab.psu.edu/research/research-areas/statistical-genetics-and-gen-epi/methods/genomesimla
GenomeSimla	GenomeSIMLA is currently under development- however, we have a beta release that we are asking to be tested http://chgr.mc.vanderbilt.edu/genomesimla/
GENS2	Simulates interactions among two genetic and one environmental factor and also allows for epistatic interactions. https://sourceforge.net/projects/gensim/
GWAsimulator	A rapid whole genome simulation program http://biostat.mc.vanderbilt.edu/wiki/main/gwasimulator
HAP-SAMPLE	An association simulator for candidate regions or genome scans http://www.hapsample.org/
HAPGEN	A simulator for the simulation of case control datasets at SNP markers https://mathgen.stats.ox.ac.uk/genetics_software/hapgen/hapgen2.html
HapSim	A simulation tool for generating haplotype data with pre-specified allele frequencies and LD coefficients http://cran.r-project.org/web/packages/hapsim/index.html
HAPSIMU	A program that simulates heterogeneous populations with various known and controllable structures under the continuous migration model or the discrete model http://l.web.umkc.edu/liujian/
IBDsim	IBDSim is a computer package for the simulation of genotypic data under general isolation by distance models. http://raphael.leblois.free.fr/
indel-Seq-Gen	A biological sequence simulation program that simulates highly divergent DNA sequences and protein superfamilies http://bioinfolab.unl.edu/~cstrope/isg/
Indelible	A powerful and flexible simulator of biological evolution http://abacus.gene.ucl.ac.uk/software/indelible/
invertFREGENE	InvertFREGENE is a forward-in-time simulator of inversions in population genetic data http://www.ebi.ac.uk/projects/bargen/
kernalPop	A spatially explicit population genetic simulation engine http://cran.r-project.org/src/contrib/archive/kernelpop/
MaCS	Markovian Coalescent Simulator http://www-hsc.usc.edu/~garykche/
Marlin	Marlin provides a user-friendly interface for performing forward-in-time population genetic simulations. http://www.patrickmeirmans.com/software/marlin.html
Mason	A package for the simulation of nucleotide data. http://www.seqan.de/projects/mason/
mbs	modifying Hudson's ms software to generate samples of DNA sequences with a biallelic site under selection http://www.sendou.soken.ac.jp/esb/innan/innanlab/software.html
Mendel's Accountant	Mendel's Accountant (MENDEL) is an advanced numerical simulation program for modeling genetic change over time and was developed collaboratively by Sanford, Baumgardner, Brewer, Gibson and ReMine http://mendelsaccount.sourceforge.net/
MetaPopGen	Simulates genetics in large size metapopulations https://sites.google.com/site/marcoandrello/metapopgen
MetaSim	A tool to generate collections of synthetic reads that reflect the diverse taxonomical composition of typical metagenome data sets http://ab.inf.uni-tuebingen.de/software/metasim/
mlcoalsim	Multilocus Coalescent Simulations http://code.google.com/p/mlcoalsim-v1/
ms	The purpose of this program is to allow one to investigate the statistical properties of such samples, to evaluate estimators or statistical tests, and generally to aid in the interpretation of polymorphism data sets. http://home.uchicago.edu/~rhudson1/source/mksamples.html
msHOT	The purpose of this program is to allow one to investigate the statistical properties of such samples, to evaluate estimators or statistical tests, and generally to aid in the interpretation of polymorphism data sets. http://home.uchicago.edu/~rhudson1/
msms	A coalescent Simlation tool with selection. http://www.mabs.at/ewing/msms/index.shtml
MySSP	A program for the simulation of DNA sequence evolution across a phylogenetic tree http://www.rosenberglab.net/software.html
Nemo	A forward-time, individual-based, genetically explicit, and stochastic simulation program designed to study the evolution of genetic markers, life history traits, and phenotypic traits in a flexible (meta-)population framework. http://nemo2.sourceforge.net/
NetRecodon	Coalescent simulation of coding DNA sequences with recombination (inter and intracodon), migration and demography http://code.google.com/p/netrecodon/
OncoSimulR	BioConductor package for Forward Genetic Simulation of Cancer Progresion with Epistasis https://github.com/rdiaz02/oncosimul
PEDAGOG	Software for simulating eco-evolutionary population dynamics https://bcrc.bio.umass.edu/pedigreesoftware/node/5
phenosim	A tool to add phenotypes to simulated genotypes http://evoplant.uni-hohenheim.de/doku.php?id=software:software
PhyloSim	An R package for the Monte Carlo simulation of sequence evolution http://www.ebi.ac.uk/goldman-srv/phylosim/
pIRS	Profile-based Illumina pair-end reads simulator https://code.google.com/p/pirs/
ProteinEvolver	Simulation of protein evolution along phylogenies under structure-based substitution models http://code.google.com/p/proteinevolver/
QMSim	QTL and Marker Simulator http://www.aps.uoguelph.ca/~msargol/qmsim/
quantiNEMO	An individual-based program for the analysis of quantitative traits with explicit genetic architecture potentially under selection in a structured population http://www2.unil.ch/popgen/softwares/quantinemo/
RECOAL	Simulates new haplotype data from a reference population of haplotypes. ftp://popgen.usc.edu/
Recodon	Coalescent simulation of coding DNA sequences with recombination, migration and demography http://code.google.com/p/recodon/
rlsim	A package for simulating RNA-seq library preparation with parameter estimation http://bit.ly/rlsim-git
Rmetasim	Rmetasim is a front-end for the metasim engine that is implemented as a package that runs in the statistical computing environment R http://cran.r-project.org/web/packages/rmetasim/index.html
RNA Seq Simulator	RSS takes SAM alignment files from RNA-Seq data and simulates over dispersed, multiple replica, differential, non-stranded RNA-Seq datasets. http://useq.sourceforge.net/cmdlnmenus.html#rnaseqsimulator
Rose	Random model of sequence evolution http://bibiserv.techfak.uni-bielefeld.de/rose/
scrm	A coalescent simulator optimized for long sequences and large samples. https://scrm.github.io/
SelSim	SelSim is a program for Monte Carlo simulation of DNA polymorphism data for a recom- bining region within which a single bi-allelic site has experienced natural selection http://www.well.ox.ac.uk/~spencer/selsim/
Seq-Gen	An application for the Monte Carlo simulation of molecular sequence evolution along phylogenetic trees. http://tree.bio.ed.ac.uk/software/seqgen/
SEQPower	Statistical power analysis for sequence-based association studies http://bioinformatics.org/spower/
SeqSIMLA	SeqSIMLA can simulate sequence data with user-specified disease and quantitative trait models. Family or unrelated case-control data can be simulated. http://seqsimla.sourceforge.net/
Serial NetEvolve	A flexible utility for generating serially-sampled sequences along a tree or recombinant network http://biorg.cis.fiu.edu/sne/
SFS_CODE	SFS_CODE can perform forward population genetic simulations under a general Wright-Fisher model with arbitrary migration, demographic, selective, and mutational effects. http://sfscode.sourceforge.net/sfs_code/index/index.html
SIBSIM	Quantitative phenotype simulation in extended pedigrees http://sourceforge.net/projects/sibsim/
SimAdapt	A spatially explicit, individual-based, forward-time, landscape-genetic simulation model combined with a landscape cellular automaton. https://www.openabm.org/model/3137
SIMCOAL2	A coalescent program for the simulation of complex recombination patterns over large genomic regions under various demographic models http://cmpg.unibe.ch/software/simcoal2/
SimCopy	An R package simulating the evolution of copy number profiles along a tree. http://bit.ly/simcopy
SIMLA	SIMLA is a SIMuLAtion program that generates data sets of families for use in Linkage and Association studies. http://dmpi.duke.edu/simla-simulation-software-version-32
SimPed	A Simulation Program to Generate Haplotype and Genotype Data for Pedigree Structures http://bioinformatics.org/simped/
Simprot	A program to simulate protein evolution by substitution, insertion and deletion http://www.uhnresearch.ca/labs/tillier/software.htm#3
SimRare	Rare variant simulation and analysis tool http://code.google.com/p/simrare/
simuGWAS	A forward-time simulator that simulates realistic samples for genome-wide association studies. http://simupop.sourceforge.net/cookbook/simugwas
simuPOP	simuPOP is a general-purpose individual-based forward-time population genetics simulation environment. http://simupop.sourceforge.net/
SISSI	A software tool to generate data of related sequences along a given phylogeny, taking into account user defined system of neighbourhoods and instantaneous rate matrices. http://www.cibiv.at/software/sissi/
SMARTPOP	Simulating Mating Alliance as a Reproductive Tactic for Populations http://smartpop.sourceforge.net/
SNPsim	Coalescent simulation of hotspot recombination http://code.google.com/p/phylosoftware/
SPIP	SPIP simulates the transmission of genes from parents to offspring in a population having demographic structure defined by the user http://swfsc.noaa.gov/textblock.aspx?division=fed&id=3434
Splatche	Spatial and Temporal Coalescences in Heterogeneous Environment http://www.splatche.com/
srv	Simulator of Rare Varaints (srv) is a simulator for the simulation of the introduction and evolution of (rare) genetic variants. http://simupop.sourceforge.net/cookbook/simurarevariants
SUP	SLINK/FastSLINK utility program http://mlemire.freeshell.org/software.html
TreesimJ	A flexible, forward-time population genetic simulator http://code.google.com/p/treesimj/
Variant Simulation Tools	A simulation tool for post-GWAS genetic epidemiological studies using whole-genome or whole-exome next-gen sequencing data, with an emphasis on user-friendliness and reproducibility. http://varianttools.sourceforge.net/simulation/homepage
Vortex	VORTEX is an individual-based simulation model for population viability analysis (PVA). http://www.vortex9.org/vortex.html
Wessim	Whole Exome Sequencing SIMulator http://sak042.github.io/wessim/

Bpipe - a tool for running and managing bioinformatics pipelines

Radha Agarkar — Sat, 21 May 2016 22:42:16 -0500

Bpipe provides a platform for running big bioinformatics jobs that consist of a series of processing stages - known as 'pipelines'.

January 20th, 2016 - New! Bpipe 0.9.9 released!
Download latest, all
Documentation
Mailing List (Google Group)

Bpipe has been published in Bioinformatics! If you use Bpipe, please cite:

Sadedin S, Pope B & Oshlack A, Bpipe: A Tool for Running and Managing Bioinformatics Pipelines, Bioinformatics

Address of the bookmark: http://docs.bpipe.org/

methylKit

Jit — Fri, 03 Jun 2016 10:09:29 -0500

methylKit is an R package for DNA methylation analysis and annotation from high-throughput bisulfite sequencing. The package is designed to deal with sequencing data from RRBS and its variants, but also target-capture methods such as Agilent SureSelect methyl-seq. In addition, methylKit can deal with base-pair resolution data for 5hmC obtained from Tab-seq or oxBS-seq. It can also handle whole-genome bisulfite sequencing data if proper input format is provided.

Address of the bookmark: https://github.com/al2na/methylKit

Graph Genome Suite

Jit — Fri, 28 Oct 2016 07:59:54 -0500

Seven Bridges is the biomedical data analysis company accelerating breakthroughs in genomics research for cancer, drug development and precision medicine. We build self-improving systems to analyze millions of genomes, including the Graph Genome Suite — the most advanced population genomics tools in the world.

Address of the bookmark: https://www.sbgenomics.com/graph/

LINKS

Jit — Fri, 04 Nov 2016 06:19:01 -0500

LINKS is a genomics application for scaffolding genome assemblies with long reads, such as those produced by Oxford Nanopore Technologies Ltd. It can be used to scaffold high-quality draft genome assemblies with any long sequences (eg. ONT reads, PacBio reads, another draft genomes, etc)

Paper at https://gigascience.biomedcentral.com/articles/10.1186/s13742-015-0076-3

Address of the bookmark: https://github.com/warrenlr/LINKS/

Maq: Mapping and Assembly with Quality

Jit — Tue, 22 Nov 2016 04:51:39 -0600

Maq stands for Mapping and Assembly with Quality It builds assembly by mapping short reads to reference sequences. Maq is a project hosted by SourceForge.net. The project page is available athttp://sourceforge.net/projects/maq/. Maq is previously known as mapass2.

Run Maq Now

Follow these steps to try Maq. All you need is a reference sequence file in the FASTA format.

Prepare a reference sequence (ref.fasta). Better a bacterial genome.
Download maq, maq-data and maqview at the download page.
Copy maq, maq.pl and maq_eval.pl to the $PATH or to the same directory.
Simulate diploid reference and read sequences, map reads, call variants and evaluate the results in one go:
```
maq.pl demo ref.fasta calib-30.dat
```
where calib-30.dat is contained in maq-data.

View the alignment:

cd maqdemo/easyrun;
maqindex -i -c consensus.cns all.map;
maqview -c consensus.cns all.map

Even for advanced maq users, running `maq.pl demo' is recommended. You may find something helpful.

Address of the bookmark: http://maq.sourceforge.net

HGA

Jit — Tue, 29 Nov 2016 07:25:53 -0600

HGA tool version 1.0 This tool helps to apply the Hierarchical Genome Assembly (HGA) method. The tool will apply: 1. Partitioning a given reads dataset into a given number of partitions. 2. Assembling each partitions using a pre-specified assembler (Velvet or SPAdes in this version) and using a given kmer size. 3. Merging all the assemblies of the partition. 4. Combining all the assemblies of the partition (using velvet with kmer value of 31). 5. Finaly, re-assembling the whole dataset with the merged contigs or the combined contigs, using a given kmer size.

https://github.com/aalokaily/Hierarchical-Genome-Assembly-HGA

Address of the bookmark: https://github.com/aalokaily/Hierarchical-Genome-Assembly-HGA

Scripts

Jit — Wed, 30 Nov 2016 10:35:15 -0600

Useful script for NGS analysis.

Address of the bookmark: http://augustus.gobics.de/binaries/scripts/

e-RGA: enhanced Reference Guided Assembly of Complex Genomes

Jit — Mon, 19 Dec 2016 05:56:14 -0600

Next Generation Sequencing has totally changed genomics: we are able to produce huge amounts of data at an incredibly low cost compared to Sanger sequencing. Despite this, some old problems have become even more difficult, de novo assembly being on top of this list. Despite efforts to design tools able to assemble, de novo, an organism sequenced with short reads, the results are still far from those achievable with long reads. In this paper, we propose a novel method that aims to improve de novo assembly in the presence of a closely related reference. The idea is to combine de novo and reference-guided assembly in order to obtain enhanced results.

Address of the bookmark: http://journal.embnet.org/index.php/embnetjournal/article/view/208