Address of the bookmark: http://pevsnerlab.kennedykrieger.org/php/?q=software

1. Strand NGS
   • offers many different tools including alignment, RNA-Seq, DNA-Seq, ChIP-Seq, Small RNA-Seq, Genome Browser, visualizations, Biological Interpretation, etc. Supports workflows “one can import the sample data in FASTA, FASTQ or tag-count format. In addition, prealigned data in SAM, BAM or Illumina-specific ELAND format can be directly imported for analysis.”
   • Alignment feature: Supports alignment from Illumina, Ion Torrent, 454 (Roche), and Pac Bio
   • DNA-Seq Feature, can annotate with dbSNP
2. CLC Genomics Workbench
   
   • (QIAGEN). Features include: resequencing, workflow, read mapping, de novo assembly, variant detection, RNA-Seq, ChIP-Seq, Genome Browser, etc (entire list on website); Main Workbench offers database search (Genbank, Blast, Pubmed); 2000 organizations have invested in CLC
   • Accepts VCF files from 1000 Genomes Project
   • Accepts downloaded tracks from dbSNP
   • Also accepts: FASTA, GFF/GTF/GVF, BED, Wiggle, Cosmic, UCSC variant database, complete genomics master var file
   • Read mapping: “In addition to Sanger sequence data, reads from these high-throughput sequencing machines are supported: The 454 FLX System and the 454 GS Junior System from Roche, Illumina Genome Analyzer, Illumina HiSeq, Illumina HiScan, and Illumina MiSeq sequencing systems, SOLiD system from Life Technologies, Ion Torrent system from Life Technologies, Helicos from Helicos BioSciences”
   • De novo assembly: “In addition to Sanger sequence data, reads from these high-throughput sequencing machines are supported The 454 FLX System and the 454 GS Junior System from Roche, Illumina Genome Analyzer, Illumina HiSeq, Illumina HiScan, and Illumina MiSeq sequencing systems, SOLiD system from Life Technologies, Ion Torrent system from Life Technologies”
   • Annotation tracks from Ensembl
3. DNAnexus
   • Private cloud repository -- formerly a redistributor of SRA and other NCBI resources; command-line or via web, can fetch data from a URL, build custom pipeline/ workflow has sra.dnanexus.com site: data downloads come directly from NCBI
4. Ingenuity Variant Analysis
   • (QIAGEN) allows for variant identification and analysis, uses NCI-60 data set for cancer, Supported third part informatin: Entrez Gene, RefSeq, ClinVar; gives contextual details of results instead of just A to B relationship
   • Has own database-- “knowledge base” based on COSMIC, OMIM, and TCGA databases
5. Lasergene Genomics Suite
   • Comprehensive NGS software pipeline for assembly, alignment, variant calling and analysis of NGS data
   • Supported workflows include: reference-guided and de novo genome and transcriptome assembly and analysis, metagenomics sample assembly, targeted resequencing, exome alignment, gene panels with validation control, variant analysis, and RNA-Seq, ChIP-Seq and miRNA alignment and analysis.
   • #1 in accuracy: fewer false negatives and better sensitivity compared to results obtained from other aligners
   • Aligns exome data and performs variant calling an average of 3 times faster than alternative pipelines
   • Annotates genomic data with allele and genotype frequency, functional impact predictions, evolutionary conservation scores and pathogenicity
   • Supports all major NGS technologies (Illumina, Ion Torrent, Pac Bio and Roche 454) and project types
   • Available on Windows, Mac OS X, Linux, and the Amazon Cloud
6. NextGENe
   • “perfect analytical partner for the analysis of desktop sequencing data produced by the ION PGM™, Roche Junior, Illumina MiSeq as well as high throughput systems as the Ion Torrent Proton, Roche FLX, Applied BioSystems SOLiD™ and Illumina® platforms.” runs on Windows, free-standing multi-application package-- SNP/Indel analysis, CNV prediction and disease discovery, whole genome alignment, etc.
   • Data can be imported from Clinvar, dbSNP, Genbank:http://www.softgenetics.com/PDF/NextGene_UsersManual_web.pdf
7. Partek Genomics Suite
   • Cited in over 3,500 peer-reviewed scientific publications
   • Workflows for microarray and PCR data include: Gene expression including alternative splicing, miRNA expression, Genome Wide Association Studies, Mother-Father-Child Trio analysis, DNA Copy number including allele specific copy number and Loss of Heterozygosity (LOH), and ChIP, and methylation. Next Generation Sequencing (NGS) workflows include: RNA-Seq, miRNA-Seq, ChIP-Seq, DNA-Seq, and Methylation
   • Powerful statistics and interactive, publication ready visualizations
   • Supports all commercial next generation sequencing and microarray file format as well as text files
   • Can input GEO SOFT files
8. Partek Flow
   • Installation can be cloud-based or on a local cluster or Linux server
   • Easy to use point-and-click interface
   • Takes NGS data (.fastq, BAM, SAM), microarrays (Affymetrix, Illumina) and text files
   • Supports custom genome builds and annotation databases
   • Performs base trimming, alignment, quantification, quality analysis, statistics, and visualization
   • Includes ten fully customizable aligners (Bowtie, Bowtie 2, BWA, GSNAP, Isaac 2, SHRiMP 2, STAR, TMAP, TopHat and TopHat 2)
   • Applications for RNA-Seq, Small RNA-Seq, WGS/WES, Pathway enrichment, Fusion detection and Variant calling
   • Allows users to create, save, share, or download analysis pipelines for automated and repeatable analysis
   • Collaborate with others without transferring data
   • Integrates microarray and next generation sequencing data
9. Golden Helix: SNP and Variation Suite
   • used for managing, analyzing and visualizing genotypic and phenotypic data; Features: Genome-wide association studies, genomic prediction, copy number analysis, small sample DNA-Seq workflows, large sample DNA-seq analysis, RNA-seq analysis. Supported files: .txt, excel XLS & XLSX, CEL, CHP, CNT, Illumina, Plink PED, TPED, BED, Agilent files, NimbleGen data summary files, VCF files, Impute2 GWAS files, HapMap format, MACH output, + 50 other formats consumes NCBI data directly
10. Genomatix
    • Applications: ChIP-Seq, DNA-Seq, RNA-Seq, DNA methylation; enable personalized medicine,
    • Mining Stations: Supports all established NGS sequencing platforms- SOLiD, 454 Life Sciences, Genome Analyzer, HiSeq, MiSeq, IonTorrent
    • Software Suite: can upload sequence of BED files
    • Genome browser: BED and BAM files, Public data- 1500 BED files available for every user
11. Biodatomics
    • Open source platform (SaaS), analysis and genome sequencing tools, integrates over 400 genomic analysis open source tools and pipelines, have a private and public cloud version. Features: genomic data visualization, drag and drop interface, accelerated analysis, real-time collaboration
    • They have a couple modules to do so, and have enabled parts of the sra toolkit
12. SolveBio
    • Software product, for clinical genomics professionals, manage, curate, report genomic variation
    • Has own data library -- data from NCBI
13. Basepair
    • Offers high quality workflows for all common NGS applications (RNA-Seq, ChIP-Seq, DNA-Seq, etc.)
    • Very fast - get all results in a 1-2 hours. Cloud-based, no storage or computing limits.
    • Easy to use - less than a minute to run an analysis
    • REST and Python API to mange large projects.
     

Variant Identification

Germline Callers

1. IMPUTE2
   • Description: phasing observed genotypes and imputing missing genotypes uses reference panels to provide all available halotypes, does not use population labels or genome-wide measures; designed to represent variation in one population; Fairly popular
   • Input:
   • Reference Haplotypes: Links to 1000 Genomes and HapMap downloads
   • Output:
2. FreeBayes
   • Description: finds SNPs, Indels, MNPs; reports variants based on alignment; haplotype based
   • Input: BAM- uses BAMtools API to parse
   • Reference genome: FASTA
   • Output: VCF
3. SOAPindel
   • Description: detects indels from NGS paired-end sequencing
   • Input: files with read alignment can be SOAP or SAM formats, users must also give raw reads in Fasta or Fastq
   • Reference Sequence used to align reads: FASTA
   • Output:
4. 2Kplus2
   • Description: algorithm searches graphs produced by de novo assembler Cortex; c++ source code for SNP detection “2kplus2.cpp is a c++ source code for the detection and the classification of single nucleotide polymorphisms in transformed De Bruijn graphs using Cortex assembler.”
   • Input:
   • Output:
5. Atlas 2
   • Description: specializes in separation of true SNPs and indels from sequencing and mapping errors, last update January 2013
   • Input: takes BAM file,
   • Reference Genome: FASTA
   • Output: produces VCF
6. CRISP
   • Description: identifies SNPs and INDELs from pooled high-throughput NGS, not used for analysis of single samples; implemented in C and uses SAMtools API; latest version should work with diploid genomes
   • Input: requires BAM files (aligned with GATK)
   • Reference Genome: indexed FASTA file
   • Output: VCF files
7. Dindel
   • Description: (Wellcome Trust Sanger) calls small indels from short-read sequences, only can handle Illumina data; cannot test candidate indels; written in C++, used on Linux based and Mac computers (not tested in windows)
   • Input: BAM files
   • Output: VCF
8. discoSnp++
   • Description: detects homozygous and heterozygous SNPs and Indels; software composed of 2 modules (kissnp2 and kissreads)
   • Input: raw NGS datasets; fasta, fastq, gzipped or not;
   • no reference genome required; read pairs can be given
   • Output: FASTA
9. FamSeq
   • Description: family-based sequencing studies- provides probability of an individual carrying variant based on family’s raw measurements; accommodates de novo mutations, can perform variant calling at chrX;
   • Input: VCF
   • Output: VCF
10. GeneticThesaurus
    • Description: “Annotation of genetic variants in repetitive regions”
    • Input: Initial variant calling from bam → vcf output
    • Reference Genome: need to provide own fasta file for hg19 genome,
    • Output: vcf.gz, vtf.gz, and baf.tsv.gz output
11. glfMultiples
    • Description: command-line, variant caller
    • Input: GLF
    • Output: VCF
12. glfSingle
    • Description: uses likelihood-based model for variant calling, starts from genotype likelihoods that have been computed from other tools (ex. Samtools BAQ), the likelihoods combine with individual-based prior p(genotype) to generate posterior probabilities
    • Input: GLF
    • Output: VCF
13. Halvade
    • Description: command-line; written in Java, “to run halvade a reference is needed for both GATK and BWA and a SNP (dbSNP!) database is required
    • Input: FASTQ
    • Output: VCF
14. indelMINER
    • Description: identifies indels from paired-end reads
    • Input: BAM (aligned in SAMtools API)
    • Output: VCF
15. Indelocator
    • Description: (Broad Institute): does not perform realignment, relies on alignments in BAM files (BAM files need aligned before put into indelocator); recommended to use GATK prior;
    • Input: 2 BAM files(tumor & normal), annotated as germline or somatic; also has single sample mode
    • Output: “Output of Indelocator is a high-sensitivity list of putative indel events containing large numbers of false positives. The statistics reported for each event have to be used to custom-filter the list in order to lower false positive rate”
16. Isaac Variant Caller
    • Description: detects SNPs and small indels from diploid sample; designed to run on “nux-like platforms”
    • Input: BAM
    • Output: VCF
17. KvarQ
    • Description: in silico genotyping for selected loci in bacterial genome, written in Python and C
    • Input: FASTQ
    • reference genome or de novo assembly not needed
    • Output:
18. LoFreq
    • Description: SNV caller, Python language, standalone program, uncovers cell-population heterogeneity from high-throughput sequencing datasets; calls variants found in <.05% of the population
    • Input: BAM file input→ suggest running through GATK
    • Output:
19. Manta
    • Description: Calls indels and SVs from paired end reads; standalone, command line program; Written in C++ and Python
    • Input: BAM (can tolerate non-paired-end reads); a matched tumor sample may be provided as well
    • Output: VCF
20. MarginAlign
    • Description: SNV caller, specifically tailored to Oxford Nanopore Reads, written in Python; Package comes with 3 programs, marginAlign, marginCaller (calls SNVs), marginStats (computes qc stats on sam files)
    • Input: SAM
    • Output: SAM
21. MendelScan
    • Description: Last release March 2014; for analyzing sequencing data in family studies of inherited diseases; variant calls for a family in VCF file; still in alpha-testing on github, example data uses 1000 genomes dataset
    • Input:
    • Output:
22. nanopore
    • Description: UCSC Nanopore group (group at UCSC studying using ion channels for analysis of single RNA/DNA structures) software pipeline; tailored to Oxford Nanopore Reads; command line program
    • Input: FASTQ
    • Reference files: FASTA
    • Output: “For each possible pair of read file, reference genome and mapping algorithm an experiment directory will be created in the nanopore/output directory.”
23. Platypus
    • Description: Package program, written in C, Python, Cython; Can identify SNPs, MNPs, short indels, and larger variants; has been tested on very large datasets (1000 genomes)
    • Input: BAM
    • Reference Genome: FASTA (files must be indexed using Samtools or similar program
    • Output: VCF
24. QualitySNPng
    • Description: detection of SNPs; “can be used as a standalone application with graphical user interface as part of pipeline system”; does not require fully sequenced reference genome; haplotype strategy
    • Input:SAM, ACE
    • Output: GUI
25. ReviSTER
    • Description: command line program; automated pipeline; utilizes BWA, BLAT, and SAMTools; utilizes BWA mapping program;
    • Input: FASTQ,
    • Reference sequence file and list file containing STR locations as inputs
    • Output: SAM
26. RVD
    • Description: command-line program, detection of rare SNVs, relies upon Samtools, can be run in MATLAB
    • Input: BAM
    • Reference Genome: FASTA
    • Output: “The algorithm output is a call table -- a comma-separated file with one line for each base position and each line in the following format:
    • AlginmentReferencePosition, AlignmentBase, Call ,SecondBase, CenteredErrorPrc, ReferenceErrorPrc, SecondBasePrc”
27. SNVer
    • Description: calls common and rare variants in pool or individual NGS data, reports overall p-value, operating system independent statistical tool, identifies SNPs and INDELs, written in Java, no dependencies, straightforward command-line
    • (SNVerGUI=GUI version) --SNVerGUI: desktop tool for variant detection
    • Input: chrX annotation, sam.zip, bam.zip
    • reference file must be aligned to the data file
    • Output:
28. SNVMix
    • Description: detects SNVs from NGS, post-alignment tool
    • Input: pileupformat (Maq or Samtools)
    • Output:
29. SV-M
    • Description: Structural Variant Machine - predicts indels, uses split read alignment profiles, validated by Sanger Sequencng
    • Input:paired-end Illumina reads from 1001 genomes project (uses ref plant- 1001genomes.org)
    • Ouptut:
30. SNPest
    • Description: Standalone program, language C++, Perl
    • Input: mpileup (SAMtools)
    • Output: VCF
31. TrioCaller
    • Description:Command line program, relies on BWA and samtools; genotype calling for unrelated individuals and parent-offspring trios
    • Input: BAM (that has been aligned in BWA and Samtools
    • Output: BCF that can be formatted to VCF using bcftools
32. Snippy
    • Description: finds indels between haploid reference genome and NGS sequence reads
    • Input:read files- FASTQ or FASTA (can be .gz compressed), output- .aln, .tab, .txt
    • Reference genome in FASTA or GENBANK
    • Output:
33. VntrSeek
    • Description: pipeline for discovering microsatellite tandem repeats with high-throughput sequencing data
    • Input: gzip-compressed FASTA or FASTQ
    • Output: VCF files; one for TRs and observed alleles, another file contains link to viewer

Somatic Callers

1. Cake
   • Description: standalone program, “pipeline for the integrated analysis of somatic variants in cancer genomes”; integrates four algorithms; written in Perl; required tools: samtools, tabix, vcftools, VarScan2, bambino, cmake, somaticsniper (User guide; workflow page)
   • Input: tumor and normal reads in BAM files, run through variant calling programs to generate intermediate VCF
   • Output: VCF
2. MuTect
   • Description: Broad Institute, identification of somatic point mutations in cancer genomes; requires preprocessing of reads (GATK)
   • Input: same as GATK (FASTA reference genome, SAM read files)
   • Output: call-stats, VCF, wiggle files
3. Polymutt
   • Description: calls SNVs and detects de novo point mutations in families
   • Input: GLF or BAM or VCF (must have identical chromosome orders)
   • Output: VCF
4. Bassovac
   • Description: Improved Bayesian inversion somatic caller; unlike other software packages, treats effects fully probabilisticallys instead of using ad-hoc modeling; effects are integrated at the atomic level and standard probability theory integrates read tallies to the sample level and to the tumor-normal pair level; "pending public release"
   • Input:
   • Output:
5. CLImAT
   • Description: standalone program; “accurate detection of copy number alteration and loss of heterozygosity in impure and aneuploid tumor samples using whole genome sequencing data”
   • Input: depth file generated by DFExtract and a config file
   • Output: .results file, .Gtype, LOG.txt, also generates visualization
6. DeNovoGear
   • Description: de-novo variant calling and interpretation; standalone program; dependencies C++ compiler, CMake, HTSlib, Eigen, Boost
   • Input: PED and BCF
   • Output: “The output format is a single row for each putative de novo mutation (DNM), with the following fields”
7. EBCall
   • Description: Empirical Baysian Mutation Calling; standalone program; uses tumor/normal paired reads and non-paired normal reference samples; dependent on samtools, R and VGAM pack for R
   • Input: BAM
   • Output: not sure what exact type of file- “The format of the result is suitable for adding annotation by annovar.”
8. HapMuc
   • Description: standalone program; “utilizes the information of heterozygous germline variants near candidate mutations”; Dependent upon- Boost, SAMtools, BEDtools; 3 step workflow
   • Input: BAM
   • Output: BED
9. MultiGeMS
   • Description: Multi-sample Genotype Model Selection
   • Input: .txt, pileup (SAM/BAM converted to pileup format)
   • Output: VCF
10. MultiSNV
    • Description: command-line program; calls SNVs from NGS data from multiple samples from the same patient; dependent on R, Git, cmake, Boost and compile libraries
    • Input: BAM or pileup
    • Output: VCF
11. MutationSeq
    • Description: standalone program, somatic SNV detection in tumor/normal samples; dependent on python, bamtools, boost, and LAPACK
    • Input: BAM
    • Output: VCF4.1 consisting of two parts (meta information & data lines)
12. qSNP
    • Description: standalone program; SNV caller for somatic variants in “low cellularity cancer samples”
    • Input: BAM, dbSNP data, Illumina data, chrConv
    • Output: “qSNP output files are named using a 4-element pattern: ...”
13. RADIA
    • Description: RNA and DNA Integrated Analysis for Somatic Mutation Detection; DNA only Method(tumor/normal pair, ignores RNA) or Triple BAM Method (uses all three datasets from same patient); dependent upon python, samtoools, pysam API, BLAT, SnpEff
    • Input: BAM
    • Reference Genome: FASTA indexed with SAMtools faidx
    • Output: VCF
14. RVD2
    • Description: sensitive, variant detection for low-depth targeted NGS data; python module or command- line program;
    • Input: tab- deliminted depth chart format (converted from pileup files)
    • Output: three hdf5 files and a vcf file
15. Shimmer
    • Description: standalone program; detects somatic SNVs with multiple testing correction, uses Fisher’s exact test; dependent on git, samtools, R, R statmod package; for tumor/normal matched samples
    • Input: BAM
    • Output: VCF
16. SNV-PPILP
    • Description: Refines GATK’s Unified Genotyper SNV calls for “multiple samples assumed to form a phylogeny”
    • Input:
    • Output:
17. SomaticSniper
    • Description: command-line application to identify SNPs between tumor/normal pairs- predicts probability of difference between two
    • Input: BAM
    • Reference Genome in FASTA
    • Output: VCF
18. Strelka
    • Description: somatic variant calling workflow for matched tumor-normal samples; detects indels; runs on *nux-like platform
    • Input: BAM (must be sorted and indexed)- Strelka does own realignment around indels-- don’t need to do this type of pre-processing
    • Output: pair of VCF files
19. Triodenovo
    • Description: Bayesian framework for calling de novo mutations in trios
    • Input: VCF file with PL or GL fields (recommend using GATK or samtools to generate)
    • Output: out_vcf
20. UNCeqr
    • Description: finds somatic mutations using integration of DNA and RNA seq data-- boosts sensitivity for low purity tumors and rare mutations;
    • Input:”can accept a variety of sequencing inputs and configurations”
    • Output: “table of somatically mutated sites and associated information. These somatic mutations can be annotated with predicted transcript and protein effects using third party tools, such as Annovar”
21. Virmid
    • Description: Virtual Microdissection for SNP calling; Java based; for disease-control matched samples; uncovers SNPs with low allele frequency by considering alpha contamination
    • Input: BAM (must be sorted and indexed- samtools sort)
    • Output: VCF and report file

Germline + Somatic Callers

1. VarScan 2
   • Description: identify germline variants, private and shared variants, somatic mutations, and somatic CNVs; detects indels
   • Input: SAMtools pileup
   • Output: VCF
2. BAYSIC
   • Description: Bayesian method; combines variant calls from different methods (GATK, FreeBayes, Atlas, Samtools, etc)
   • Input: VCF format from one or more variant calling programs
   • Output: VCF file containing integrated set of variant calls
3. MSIsensor
   • Description: Microsatellite instability detection; C++ program, detects somatic and germline variants in tumor-normal paired data
   • Input: BAM index files (normal and tumor)
   • Output:
4. Beagle version 4
   • Description: software package: genotype calling, phasing, imputation of ungenotyped markers, and identity-by-descent segment detection:unsure if this one is in the right category; genotype calling, phasing, imputation of ungenotyped markers, and identity-by-descent segment detection;
   • Input: VCF
   • Output: VCF
5. QuadGT
   • Description: software package, SNV calling from normal-tumor pair and two parent genomes; quantifies descent-by-modification relationships; Written in Java
   • Input: BAM files (parsed by Picard/Samtools API)
   • Reference Genome; FASTA
   • Output: VCF
6. RAREVATOR
   • Description: RAre REference VAriant annotaTOR; command line; “identification and annotation of germline and somatic variants in rare reference allele loci from second generation sequencing data”; Bayesian genotype likelihood model
   • Input: BED or VCF files from GATK
   • Output: two VCF files (one for SNVs, one for Indels)
7. Scalpel
   • Description: Used for detecting indels in a reference genome; performs localized micro-assembly of specific regions of interest; can do single, de novo, somatic reads; requires that raw reads are aligned with BWA
   • Input: BAM
   • Output: either VCF or ANNOVAR
8. SOAPsnp
   • Description: based on Baye’s theorem; calls consensus genotype
   • Input:SOAP short read alignment results
   • Output: GLF, option of flat tabular format
9. VariantMaster
   • Description: “extract causative variants for monogenic and sporadic genetic diseases”; uses ANNOVAR;
   • Input: BAM or VCF files (from SAMtools, GATK)
   • Output:

Downstream Analysis of Variants

1. PrediXcan
   • Description: command-line, standalone package program; available in Perl, Python, and R versions; predicts liklihood of a gene being related to a certain phenotype- “that directly tests the molecular mechanisms through which genetic variation affects phenotype.”; no actual expression data used, only in silico expression; “PrediXcan can detect known and novel genes associated with disease traits and provide insights into the mechanism of these associations.”
   • Input: genotype and phenotype file (doesn’t specify file type)
   • Output:default values: genelist, dosages (file format: snpid rsid) , dosage_prefix, weights, output
2. ATHENA
   • Description: Analysis Tool for Heritable and Environmental Network Associations; software package, combines machine learning model with biology and statistics to predict non-linear interactions
   • Input: Configuration file, Data file, Map file (includes rsID)
   • Output: Summary file, Best model file, dot file, individual score file, cross-validation file
3. CCRaVAT and QuTie
   • Description: (Wellcome Trust Sanger) Case-Control Rare Variant Analysis Tool and Quantitative Trait; software packages for large-scale analysis of rare variants
   • Input: PED file and MAP file
   • Output: Five tab-delimited txt files
4. GCTA
   • Description: Genome Wide Complex Trait Analysis; package program, command line interface; estimates variance by all SNPs; 5 main functions: “data management, estimation of the genetic relationships from SNPs, mixed linear model analysis of variance explained by the SNPs, estimation of the linkage disequilibrium structure, and GWAS simulation”
   • Input: PLINK binary PED files, MACH output format
   • Output:
5. GenomeComb
   • Description: package for analysis of complete genome data; annotation using public data or custom tracks, automated primer desing for Sanger or Sequenom validation; “The cg process_illumina command can be used to generate annotated multisample data starting from fastq files, using tools such as bwa for alignment and GATK and samtools for variant calling. Sequencing data can also be imported from Complete Genomics (cg_process_sample command), Real Time Genomics (cg_process_rtgsample command) and VariantCallFormat (VCF) variant files (vcf2sft command).”
   • Input: Sequencing data from Complete Genomics, Illumina, SOLiD and VCF;
   • Output: standard file format used is a simple tab delimited file (.sft, .tsv)
6. Genome Track Analyzer
   • Description: compares genome tracks; allows user to compare DNA expression/binding;
   • Input: multiple: SGR/TXT, BED, BED6, GFF; if using prealigned sequence data- use MACS peak caller: BAM, BED, SAM, ELAND
   • Output:
7. GVCBLUP
   • Description: animal gene mapping; “genomic prediction and variance component estimation of additive and dominance effects”; standalone program, command line interface, writting in C++ and Java
   • Input:
   • Output:
8. HOMOG
   • Description: Analyzes heterogeneity with respect to single marker loci or known maps of markers; Carries out homogeneity test for alternative hypothesis “Two family types, one with linkage betweeen a trait to a marker or map of markers, the other without linkage”
   • Input: HOMOG.DAT - described on website
   • Output: HOMOG.OUT
9. INTERSNP
   • Description: GWIA for case-control SNP and quantitative traits; selected for joint analysis using priori information; Provides linear regression framework, Pathway Association Analysis, Genome-wide Haplotype Analysis,
   • Input: PLINK input formats (ped/map, tped/tfam, bed/bim/fam) Compatible with SetID files
   • Gene reference file: Ensembl Release 75
   • Output: covariance matrix for regression models
10. mtSet
    • Description: Currently only the standalone version available, but moving to LIMIX software suite; offers set tests- allows for testing between variants and traits; accounts for confounding factors ex. relatedness
    • Input: sample-to-sample genetic covariance matrix needs to be computed; multiple types of input; simulator requires input genotype and relatedness component;
    • Output: resdir (result file of analysis), outfile (test statistics and p-values), manhattan_plot (flag)
11. MultiBLUP
    • Description: Package program, command line interface; constructs linear prediction models; Best Linear Unbiased Prediction; improves upon BLUP involving kinship matrices; options: pre-specified kinships, regional kinships, adaptive multiblups, LD weightings
    • Input: PLINK format
    • Output:.reml, .indi.blp

Variant Annotation

1. ANNOVAR
   • Description: command-line tool, supports SNPs, INDELs, CNVs and block substitutions, provides wide variety of annotation techniques, depends upon multiple databases (each needing to be downloaded); annotates genetic variants; utilizes RefSeq, UCSC Genes, and the Ensembl gene annotation systems; can compare mutations detected in dpSNP or 1000 Genomes Project; Very popular *“The final command run TABLE_ANNOVAR, using dbSNP version 138, 1000 Genomes Project 2014 Oct version, NIH-NHLBI 6500 exome database version 2 (referred to as esp6400siv2), dbNFSP version 2.6 (referred to as ljb26), dbSNP version 138 (referred to as snp138) databases and remove all temporary files, and generates the output file called myanno.hg19_multianno.txt”
   • Input: VCF, ANNOVAR input format (simple text-based format); can convert other formats into ANNOVAR input format
   • Output: VCF (if input VCF), output file with multiple columns, tab-delimited output file
2. wANNOVAR
   • provides web-based access to ANNOVAR software
3. PolyPhen-2
   • Description: Very popular; Polymorphism Phenotyping; Web application; predicts impact of amino acid substitution on protein; Calculates Bayes posterior probability (Last update July 2015)
   • Input: FASTA
   • Output:
4. SIFT
   • Description: predicts how an amino acid substitution will affect protein function; Based on degree of conservation of amino acid residues- collected though PSI-BLAST; can be applied to nonsynonymous polymorphisms or laboratory-induced missense mutations; links to dbSNP 132, GRCh37; Standalone or web app program; Very popular
   • Input: Uniprot ID or Accession, Go term ID, Function name, Species Name or ID, etc
   • Output:
5. snpEff
   • Description: Genetic variant annotation and effect prediction toolbox; integrated with Galaxy, GATK, and GNKO; can annotate SNPs, INDELs, and multiple-nucleotide polymorphisms; categorizes effects into classes by functionality; Very popular; Standalone or Web app; Claims to calculate all SNPs in 1000 genomes (EMBI) in less than 15 minutes; can annotate SNPs, MNPs, and insertions and deletions; Provides assessment of impact of the variant ( low, medium or high)
   • Input: VCF, BED
   • Output: VCF (with new ANN field, also used in ANNOVAR and VEP), HTML summary files
6. SnpSIFT
   • Description: Filter and manipulate annotated files; Part of SnpEff main distribution; one variants have been annotated, this can be used to filter your data to find relevant variants
   • Input:
   • Output:
7. VAAST 2
   • Description: Variant Annotation, Analysis, and Search Tool; probabilistic search tool for identifying damage genes and the disease causing variants; can score both coding and non-coding variants; Four tools: VAT (Variant annotation tool), VST (Variant Selection Tool), VAAST, pVAAST (for pedigree data); updated April 2015
   • Input: FASTA, GFF3, GVF
   • Output: CDR (condenser file), VAAST file (both unique to VAAST)
8. VEP
   • Description: (Ensembl) Variant Effect Predictor; determines effect of variants on genes, transcripts, and protein sequence; uses SIFT and PolyPhen
   • Input: Coordinates of variants and nucleotide changes; whitespace- separated format, VCF, pileup, HGVS
   • Output: VCF, JSON, Statistics
9. ABSOLUTE
   • Description: (Broad Institute); can estimate purity and ploidy to compute absolute copy number and mutation multiplicitie; reextracts data from the mixed DNA population
   • Input: HAPSEQ segdat or segmentation file
   • Output: per-sample output directory and subdirectory providing per-sample text files containing standard out being emitted from R
10. Alamut Batch
    • Description: high-throughput annotation software for NGS analysis; for “intensive variant analysis workflows”; “enriches raw NGS variants with dozens of attributes”; based on clinically oriented Alamut database; Supports human genes; easy to integrate into pipeline (Latest Release- July 2015)
    • Input:VCF, tab-delimted file
    • Output: tab-separated file of annotations
11. AVIA
    • Description: Annotation, Visualization, and Impact Analysis; “The tool is based on coupling a comprehensive annotation pipeline with a flexible visualization method. We leveraged the ANNOVAR (Wang et. al, 2010) framework for assigning functional impact to genomic variations by extending its list of reference annotation databases (RefSeq, UCSC, SIFT, Polyphen etc.) with additional in-house developed sources (Non-B DB, PolyBrowse).”
    • Input: BED
    • Output: Table of annotations with gene annotation features
12. BioR
    • Description: (Mayo Clinic) (Page last updated June 2015) Biological Reference Repository; “data integration tool that enables coordinate based searches and joins based on strings”; “BioR consists of two parts 1) the BioR toolkit which depends on Java…. 2) the BioR catalogs which are the data files used by the system”
    • Input: VCF
    • BioR-Supported Catalogs (tar-gzip files): dbSNP, 1000 genomes, HapMap, OMIM, NCBIGene
    • Output: VCF + JSON
13. CADD
    • Description: Combined Annotation Dependent Depletion; tool for scoring SNV deletions/insertions; “integrates multiple annotations into one metric”; Score strongly correlates with allelic diversity and pathogenicity; links to 1000 Genome variants; uses Ensembl Variant Effect Predictor
    • Input: VCF
    • Output: CADD score
14. CandiSNPer
    • Description: web application, characterizes SNPs located in vicinity of SNP of interest;
    • Input: enter SNP ID (rsID), choose population, region, measure for LD, threshold plot format, color of SNPs, and chose to show genes
    • Output: Imagefile
15. CanvasDB
    • Description: “local database infrastructure for analysis of targeted- and whole genome re-sequencing projects”; dependent on MySQL, R, and ANNOVAR
    • Input:
    • Output:
16. CAROL
    • Description: (Wellcome Trust Sanger); Combined Annotation scoRing toOL; Combined functional annotation score of nonsynonymous coding variants; Combines information from PolyPhen-2 and SIFT
    • Input: tab-delimited with columns obtained from PolyPhen-2 and SIFT output
    • Output: tab-delimited file
17. CHASM
    • Description: Cancer-specific High-throughput Annotation of Somatic Mutations; Last updated May 2014; uses Random Forest Method to “distinguish between driver and passenger somatic mutations”; Positive driver class curated from COSMIC database; packed together with SNVBox (database)
    • Input:Passenger mutation rates, Transcript and amino acid change, Genomic coordinates
    • Output: CHASM score, p-value, FDR
18. CRAVAT
    • Description: Cancer-Related Analysis of Variants Toolkit; Web application; Uses CHASM, VEST, SNVGet; “CRAVAT provides predictive scores for germline variants, somatic mutations and relative gene importance, as well as annotations from published literature and databases” Latest Release May 2015;
    • Input: VCF, CRAVAT format
    • Output: CRAVAT report- MS Excel spreadsheet or tab-separated file (emailed)
19. CUPSAT
    • Description: Cologne University Protein Stability Analysis Tool; “tool to predict changes in protein stability upon point mutations”; web service program; Can predict mutant stability from existing PDB structures or custom protein structures
    • Input:for PDB- provide PDB ID and Amino Acid Residue Number; for custom- PDB file format
    • Output:
20. DANN
    • Description: Deleterious Annotation of genetic variants; standalone program, uses “the same feature set and training data as CADD to train a deep neural network”; can catch nonlinear relationships; “There are four different datasets: training, validation, testing, and ClinVar_ESP...The ClinVar_ESP dataset is also a testing set containing a set of “gold standard” pathogenic and benign variants”
    • Input:
    • Output:
21. ESEfinder
    • Description: Exonic Splicing Enhancer; useful for interpretation of point mutations/polymorphisms that are disease-associated; GUI interface; web app program
    • Input: FASTA
    • Output: html or plain text format, graphical display of results
22. Exomiser
    • Description: Wellcome Trust Sanger; functionally annotates variants from whole-exome sequencing data; Based on Jannovar and uses UCSC KnownGene; Java program; web app program (Page last modified Feb 2015)
    • Input: VCF
    • Output: TSV, VCF
23. FamAn
    • Description: Automated variant annotation pipeline for family-based sequencing studies; Annotaties SNVs and INDELs; 4 models- autosomal dominant, autosomal recessive, de novo mutations and a general model; “A variety of annotations are provided for each segregating variant: number of family (and family ID) each variant hits, variant genomic location and coding effect (based on snpEff), loss-of-function mutation annotation, selected ENCODE annotation, allele frequency in the 1000 Genomes Project, allele frequency in the Exome Variant Server (ESP6500), segmental duplication annotation, SIFT, PolyPhen2, LRT, MutationTaster, GERP++, PhyloP, SiPhy, etc.” (Last updated May 2014)
    • Input: VCF
    • Output: two excel compatible outputs
24. GeneTalk
    • Description: Combines tool for filtering and data analysis with an online network for genetic professionals; Different degrees- basic license, premium license, in-house solution (the last ones are paid for- Commercial tool?)
    • Input: VCF
    • Output: GeneTalk Annotation- includes clinical data, medical relevance, scientific relevance (http://www.gene-talk.de/public/GeneTalk_Whitepaper_Annotations.pdf)
25. GeneVetter
    • Description: “GeneVetter is a tool designed for investigation of the background prevalence of exonic variation in the Phase 3 1000 Genomes data under user defined filtering criteria”; web app program; GeneVetter uses GRch37p4 (hs37d5.fa.gz), dbSNP build 138, 1000G Phase 3, clinvar_2014072
    • Input: VCF
    • Output: TIMS score, summary table, PCA plot
26. GSITIC
    • Description: (Broad Institute) Last update- July 2014; Identifies genomic regions that are significantly “amplified or deleted”; Each is given a G score; gives genomic locations and q-values from aberrant regions
    • Input: segmentation file -seg, markers file -mk (required); -array file list -alf, CNV file -cnv
    • Reference genome: -refgene (created in MATLAB, GISITIC provides four reference genomes: hg16.mat, hg17.mat, hg18.mat, hg19.mat
    • Output: All lesions file (text file), amplifications file (text file), deletion genes file (text file), Gistic Scores file, Segmented copy number (pdf file), amplification score GISTIC plot (pdf file), Deletion score/q-vale GISTIC plot (pdf file)
27. HOPE
    • Description: Have yOur Protein Explained; Web app program; Automatic mutant analysis server that provides structural effects of a mutation; Uses BLAST against UniProt and PDB along with homology modeling
    • Input: FASTA protein sequence, or accession code of protein of interest
    • Output: a report containing information from a “decision tree” and illustrated figures and animations
28. Human Splicing Finder
    • Description: Last update: May 2013; aimed to help study pre-mRNA splicing; combines 12 algorithms to identify mutations’ effect on splicing motifs; uses ensembl database 70
    • Input: Gene Name, Ensembl transcript ID, Ensembl Gene ID, Consensus CDS, RefSeq Peptide ID, or own sequence (looks like you can enter FASTA)
    • Output: Chart with columns for predicted signal, predicted algorithm, cDNA position and interpretation
29. LARVA
    • Description: Large-scale Analysis of Variants in noncoding Annotations; New version released July 2015; Command-line program; used for studying noncoding variants; integrates comprehensive set of noncoding elements, modeling their mutation count; Dependent on C++ and BEDtools
    • Input: multiple
    • Output:
30. LINKAGE
    • Description:three main programs: mlink (calculates lod scores at fixed values for the recombination fraction in one interval of a genetic map), linkmap (calculates location scores for positions of a disease locus along a marker), and ilink (estimates parameters including recombination fractions, allele frequencies, penetrances, etc)
    • Input: pedfile (processed by MAKEPED) and datafile (reflects loci for each individual; set in PREPLINK)
    • Output:
31. MAC
    • Description: MNV Annotation Corrector; Ad hoc software, fixes incorrect amino acid predictions that are caused by multiple nucleotide variations; Uses existing annotators ANNOVAR, SnpEff, VEP (last update April 2015) (only 1 download this week → not popular)
    • Input: List of called SNVs and corresponding BAM
    • Output: Report identifying block of mutation within codon (BMCs)
32. mit-o-matic
    • Description: focuses on mtDNA, provides clinically relevant information from different resources; two component pipeline: command link for alignment of NGS reads and online version that provides genetic report on mitocondrial variants
    • Input:FASTQ, pileup
    • Reference sequence: rCRSm
    • Output: Online version gives comprehensive genetic report
33. Mutadelic
    • Description: Web App program; “This application generates reports on inherited mutations in five genes (ANK1, SLC4A1, SPTA1, SPTB and EPB42) associated with the following rare Mendelian blood disorders: Hereditary Spherocytosis (HS), Hereditary Elliptocytosis (HE) and Hereditary Pyropoikilocytosis”; Newer program- recently validated on omictools
    • Input: Can upload coordinates of DNA variants or VEP
    • Output: Displayed on web or can be downloaded in Excel or RDF format
34. MutationTaster
    • Description: (Last post on site 2014) Web app program; Rapid evaluation of disease causing alterations; uses NCBI 37 and Ensembl 69
    • Input: HGNC symbol, NCBI GeneID, or Ensembl ID,
    • Output: Report containing prediction, summary, name of alteration, etc
35. MutPred
    • Description: web app tool; Classifies amino acids substituation as disease associated or neutral in humans; Last modified Feb. 2014; Based on SIFT, trained using Human Gene Mutation Database
    • Input:
    • Output: “The output of MutPred contains a general score (g), i.e., the probability that the amino acid substitution is deleterious/disease-associated, and top 5 property scores (p), where p is the P-value that certain structural and functional properties are impacted.”
36. MutSigCV
    • Description: (Broad Institute) Mutation Significance (CV= covariates); Analyzes mutations discovered in DNA sequencing to identify genes that were mutated more often than expected
    • Input: mutations.maf, coverage.txt, covariates.txt
    • Output: output.txt
37. NGS-SNP
    • Description: Collection of command-line scripts for providing rich SNP annotations; “NCBI, Ensembl, and Uniprot IDs are provided for genes, transcripts and proteins when applicable”;
    • Input: Samtools consensus pileup, Maq, diBayes, Genetic format, VCF
    • Output: File containing annotated SNPs is copied from SNP list and some classes are added
38. Oncotator
    • Description: (Broad Institute) “Tool for annotating human genomic point mutations and data relevant to cancer researchers”; Web app; Supports annotation of data from ClinVar, dbSNP, 1000 genomes (plus many other external sites); Only GRCh27 coordinates supported; Last update: April 2015
    • Input: tal-delimited file
    • Output: tab-delimited MAF
39. PANTHER
    • Description: Protein ANalysis THrough Evolutionary Relationships; Web app program, also has its own database; Classification system used to classify proteins and their genes; Also, “Estimates the likelihood of a particular nonsynonymous (amino-acid changing) coding SNP to cause a functional impact on the protein”; Updated in 2015
    • Input: Data from PANTHER, IDs from Ensembl, EntrezGene, NCBI GI numbers, NCBI UniGene IDs HUGO, UniProt; if ID type is not one of the above, can input txt file or excel format
    • Output: Analysis results displayed online
40. PESX
    • Description: Putative Exonic Splicing Enhancers/Silencers; (Can’t tell if this is outdated or not)
    • Input: FASTA or plain text
    • Output: Excel spread sheet
41. Phen-Gen
    • Description: Combines patient's’ disease symptoms with sequencing data; Standalone or Web app version; Only excepts 1 family per run, in order to evaluate unrelated individuals, each sample needs to be run individually
    • Input: Variant- VCF; Pheotype- HPO; Pedigree- PED
    • Output: Combined scores file, variants for top genes file
42. PMUT
    • Description: Aimed at annotation and prediction of pathological mutations; based on different kinds of sequence info and neural networks to process information
    • Input: FASTA
    • Output; Simple yes/no and reliability index
43. PROVEAN
    • Description: Protein Variation Effect Analyzer; predicts whether an amino acid substitution or indel has impact on biological function of the protein; “comparable to SIFT or Polyphen-2”; Standalone, Web app, Command line or GUI; Last update May 2014
    • Input: FASTA, list of variants;
    • Output: tab-separated columns including Variant, Provean Score and prediciton
44. Rescue-ESE
    • Description: “An online tool for identifying candidate ESEs in vertebrate exons”; Web application; For human, mouse, zebrafish, pufferfish
    • Input: multi-FASTA or plain text
    • Output:
45. SCAN
    • Description: Web application program, includes a database as well; Database contains physical-based SNP annotations and functional annotations; “Information on physical, functional, and LD annotation served on the SCAN database comes directly from public resources, including the HapMap (release 23a), NCBI (dbSNP 129), or is information created by us using data downloaded from these public resources”; “SCAN can be utilized in several ways including: (i) queries of the SNP and gene databases; (ii) analysis using the attached tools and algorithms; (iii) downloading files with SNP annotation for various GWA platforms”
    • Input:
    • Output: HTML, comma-delimited, tab-delimited
46. SeattleSeq Annotation
    • Description: “SeattleSeqAnnotation137 was most recently updated October 13, 2013. The current version is 8.08. The most recent site, based on dbSNP build 141, and hg38/NCBI 38”; Provides annotations for SNVs and Indels- includes dbSNP rsID, gene names and accession numbers, variation functions, protein positions and amino acid changes, conservation scores, HapMap frequencies, PolyPhen predictions and clinical association.
    • Input: Maq, gff, CASAVA, VCF, GATK bed, custom
    • Output: “default output file format is a header line (starting with "#") followed by tab-separated annotations”; VCF
47. seqminer 3.7
    • Description: “Efficiently Read Sequence Data (VCF Format, BCF Format and METAL Format) into R”; Command line package program; Published August 2015
    • Input: VCF, BCF
    • Output: VCF
48. SG Adviser
    • Description: Scripps Genome Annotation and Distributed Variant Interpretation Server, web developed applications for variant annotation, “Downstream applications of variant annotation include: Clinical sequencing applications including: carrier testing, or identification of causal variants in molecular diagnosis, tumor sequencing, or diagnostic odyssey. Prioritization of variants prior to statistical analysis of sequence based disease association studies, especially for automated set-generation and enrichment of likely functional variants within sets. Identification of causal variants in post-GWAS/linkage sequencing studies. Identification of causal variants in forward genetic screens (stay tuned for non-human annotation)”
    • Input: SNV- VCF, BED, and a few others; CNV- BED, CNVator, plus others
    • Output: tab-delimited file
49. SNAP-2
    • Descriptio

More at https://gatb.inria.fr/

Address of the bookmark: https://gatb.inria.fr/

LPI scores are inversely proportional to the phylogenetic distance between database match sequences and the query genome. These scores are useful not only for large-scalede novo predictions of horizontally transferred proteins, but can also serve as an independent quality control test for potential horizontal transfer candidates identified by alternative methods, especially those based on nucleic acid signatures. Candidates having high LPI scores are unlikely to have been horizontally transferred, since they are highly conserved among closely related organisms.

One unique and powerful feature of the DarkHorse HGT Candidate database is the opportunity to explore the phylogenetic background of potential HGT donors as well as recipients. The breadth of the database allows not only query sequences, but also their database match partners to be evaluated for sequence similarity or novelty compared to taxonomically related organisms.

DarkHorse is configurable for varying degrees of phylogenetic granularity and protein sequence conservation. Users should consult the references cited below for a complete explanation of parameter selection and result interpretation. A brief tutorial page is also available on-line.

Address of the bookmark: http://darkhorse.ucsd.edu/download.html

We therefore list many genomle assembly tools here. We mainly reported for the assembly of genomes while the others are designed aiming at handling complex genomes.

• TriMetAss 1.2 – The Trinity-based Iterative Metagenomics Assembler
  • TriMetAss is an extension to the Trinity software [1], which can assemble select regions surrounding interesting features in metagenomic data. The software is particularly useful for very common and well-conserved genes (and – in theory – non-coding regions) that can occur in multiple contexts in the microbial community under study. It uses Vmatch [2] to extend seed reads (or contigs generated by another assembler) into longer contigs, by iteratively calling Vmatch and Trinity, until some stop criteria are met. Currently, TriMetAss lacks a thorough documentation, but you can direct questions to me if the README.txt file and the “-h” option is not sufficient to understand the software.
    
    
• OMWare 1.0 – Efficient Assembly of Genome-wide Physical Maps

  • The purpose of this Python module is help scientists use optical map data.
    Once complete, it will encapsulate and abstractify optical maps and their most common manipulations as they exist in a variety of formats.

• LightAssembler – Lightweight Resources Assembly Algorithm

  • Lightweight resources assembly algorithm for high-throughput sequencing reads.
    System requirements
    64-bit machine with g++ compiler or gcc in general, pthreads,and zlib libraries.

• QUAST 4.1 – Quality Assessment Tool for Genome Assemblies

  • QUAST evaluates genome assemblies.
    QUAST works both with and without a reference genome. 
    The tool accepts multiple assemblies, thus is suitable for comparison.

• DNA Baser 4.36 – DNA Sequence Assembly & Analysis
  • DNA Sequence Assembler is revolutionary bioinformatics software for automatic DNA sequence assembly , DNA sequence analysis, contig editing, file format conversion and mutation detection.
    
    
• COCACOLA – Binning Metagenomic Contigs using Sequence COmposition, Read CoverAge, CO-alignment, and Paired-end Read LinkAge
  
  • COCACOLA: a general framework for binning contigs in metagenomic studies incorporating read COverage, CorrelAtion, sequence COmposition and paired-end read LinkAge
    
    
• MaxBin 2.2 – Binning Assembled Metagenomic Sequences
  • MaxBin is software for binning assembled metagenomic sequences based on an Expectation-Maximization algorithm. Users can understand the underlying bins (genomes) of the microbes in their metagenomes by simply providing assembled metagenomic sequences and the reads coverage information or sequencing reads. 
    
    
• GAML 0.1 – Genome Assembly by Maximum Likelihood
  
  • GAML is a prototype genome assembly tool based on maximizing likelihood of the assembly in a model encompaasing error rate, insert length and other features of indvidual sequencing technologies. It can combine datasets produced by different technologies (currently Illumina, 454 and Pacific Biosciences).
    
    
• NanoMark – DNA Assembly Benchmark for Nanopore long reads

  • DNA Assembly Benchmark for Nanopore long reads
    A system for benchmarking DNA assembly tools, based on 3rd generation sequencers.

• ARC 1.1.4-beta – Assembly by Reduced Complexity

  • ARC is a pipeline which facilitates iterative, reference guided de novo assemblies with the intent of: 
    1.Reducing time in analysis and increasing accuracy of results by only considering those reads which should assemble together.
    2.Reducing/removing reference bias as compared to mapping based approaches.

• TransPS 1.1.0 – Transcriptome Post Scaffolding
  • TransPS is a pipeline for post-processing of pre-assembled transcriptomes using reference based method. It applies an align-layout-consensus structure, consisting of three major stages. First, query sequences are aligned with a reference genome. Second, query sequences are ordered based on the alignment to the reference. Third, non-redundant sequences matched to the same gene of reference genome are scaffolded into one contig. 
    
    
• assemblyManager – Computing the Robotic Commands for 2ab Assembly
  • Clotho provides persistence to such objects through relational databases that at least partially correspond the Clotho data model. Beyond database access and data model API support, Clotho Apps provide more specific functionality to Clotho such as viewing and editing data, running simulations, and automating various tasks. When thinking about Clotho Apps, an appropriate analogy would be Apps running on the Android operating system rather than the add-ons that extend the functionality of Firefox
    
    
• BinPacker 1.1 – Packing-Based De Novo Transcriptome Assembly from RNA-seq Data
  • BinPacker is a novel de novo assembler by modeling the transcriptome assembly problem as tracking a set of trajectories of items with their sizes representing coverage of their corresponding isoforms by solving a series of bin-packing problems
    
    
• FermiKit 0.13 – De novo Assembly based Variant Calling pipeline for Illumina Short Reads
  • FermiKit is a de novo assembly based variant calling pipeline for deep Illumina resequencing data. It assembles reads into unitigs, maps them to the reference genome and then calls variants from the alignment to an accuracy comparable to conventional mapping based pipelines (see evaluation in the tex directory). The assembly does not only encode SNPs and short INDELs, but also retains long deletions, novel sequence insertions, translocations and copy numbers
    
    
• REPdenovo – A tool to Construct Repeats directly from Raw Reads

  • REPdenovo is designed for constructing repeats directly from sequence reads. It based on the idea of frequent k-mer assembly. REPdenovo provides many functionalities, and can generate much longer repeats than existing tools. The overall pipeline is shown in the mannual file. REPdenovo supports the following main functionalities.
    1.Assembly. This step performs k-mer counting. Then we find frequent k-mers whose frequencies are over certain threshold. We then assemble these frequent k-mers into consensus repeats (in the form of contigs). Then we merge the constructed contigs to more completeness ones.
    2.Scaffolding. We use paired-end reads to connect repeat contigs into scaffolds, also provide the average coverage (indicates the copy number) for each constructed repeats.

• Xander – Gene-targeted Metagenomic Assembler
  • Metagenomics can provide important insight into microbial communities. However, assembling metagenomic datasets has proven to be computationally challenging. We present a novel method for targeting assembly of specific protein-coding genes using a graph structure combining both de Bruijn graphs and protein HMMs. The inclusion of HMM information guides the assembly, with concomitant gene annotation. 
    
    
• SWAP-Assembler 2 – A scalable and fully parallelized Genome Assembler
  • There is a growing gap between the output of new generation massively parallel sequencing machines and the ability to process and analyze the sequencing data. We present SWAP-Assembler, a scalable and fully parallelized genome assembler designed for massive sequencing data. Intend of using traditional de Bruijn Graph, SWAP-Assembler adopts multi-step bi-directed graph (MSG). With MSG, the standard genome assembly (SGA) is equivalent to the edge merging operations in a semi-group. Then a computation model, SWAP, is designed to parallelize semi-group computation. Experimental results showed that SWAP-Assembler is the fastest and most efficient assemblers ever, it can generated contigs with highest accuracy over all five selected assemblers and longest contig N50 in all selected parallel assemblers. Specially, in the scalability test, SWAP-Assembler can scales up to 1024 cores when processing Fish and Yanhuang dataset, and finishes the assembly work in only 15 and 29 minutes respecitively
    
    
• TGNet – Visualization and Quality Assessment of de novo Genome Assemblies
  • TGNet is a Cytoscape-based tool for visualization and quality assessment of de novo genome assemblies. Specifically it facilitates rapid detection of inconsistencies between a genome assembly and an independently derived transcriptome assembly.
    
    
• Circlator 1.1.3 – A tool to Circularize Genome Assemblies
  • A tool to circularize genome assemblies. The algorithm and benchmarks are described in the Genome Biology manuscript. 
    
    
• misFinder v0.4.05.05 – Identify Mis-assemblies in an unbiased manner using Reference and Paired-end Reads
  • misFinder is a tool that aims to identify the assembly errors with high accuracy in an unbiased way and correct these errors at their mis-assembled positions to improve the assembly accuracy for downstream analysis. It combines the information of reference (or close related reference) genome and aligned paired-end reads to the assembled sequence. Structure variation and mis-assembly can be detected by comparing the reference genome and assembled sequence.
    
    
• Scaffold_builder v2.2 – Order Contigs generated by draft sequencing along a Reference Sequence
  • The abundance of repeat elements in genomes can impede the assembly of a single sequence. The tool Scaffold_builder was designed to generate scaffolds (super contigs of sequences joined by N-bases) using the homology provided by a closely related reference sequence. Scaffold_builder is an advanced wrapper for Nucmer, written in Python that resolves several situations that may arise when mapping contigs to the reference genome.
    
    
• Rnnotator 3.5.0 – de novo Transcriptome Assembly pipeline from stranded RNA-Seq reads
  • Comprehensive annotation and quantification of transcriptomes are outstanding problems in functional genomics. Rnnotator is an automated software pipeline that generates transcript models by de novo assembly of RNA-Seq data without the need for a reference genome. The contigs produced by Rnnotator are highly accurate and reconstruct full-length genes when transcripts are sequenced sufficiently deep, roughly 30X for a given transcript. Rnnotator was designed to assemble Illumina single or paired-end reads. Rnnotator is also able to incorporate strand-specific RNA-Seq reads into the assembly in order to further improve the assembly.
    
    
• SATRAP 0.2 – SOLiD Assembler TRAnslation Program

  • A color space assembly must be translated into bases before applying bioinformatics analyses. SATRAP is designed to accomplish this important task adopting a very efficient strategy. The package integrates the Oases pipeline and several optimizations specifically designed for color space management. All steps of the pipeline allow to produce a SOLiD de novo transcriptome assembly and the subsequent color space translation. Alternatively, SATRAP can be used as a stand alone program to perform color space translation for either RNA-seq or DNA-seq SOLiD assemblies.

• Bandage v0.7.1 – Navigating De novo Assembly Graphs Easily
  • Bandage is a program for visualising de novo assembly graphs. By displaying connections which are not present in the contigs file, Bandage opens up new possibilities for analysing de novo assemblies.
    
    
• HapCol 1.1.1 – Haplotype Assembly from Long Gapless Reads
  • A fast and memory-efficient method for haplotype assembly from long gapless reads, like those produced by SMRT sequencing technologies (PacBio RS II) and Oxford Nanopore flow cell technologies (MinION).
    
    
• REAGO 1.1 – REconstruct 16S ribosomal RNA Genes from MetagenOmic data
  
  • an assembly tool for 16S ribosomal RNA recovery from metagenomic data
    
    
• FGAP 1.8.1 – Automated Gap Closing tool
  • FGAP aims to improve genome sequences by merging alternative assemblies or incorporating alternative data, analyzing the gap region and indicating the best sequence to close the gap.
    
    
• DETONATE 1.10 – DE novo TranscriptOme rNa-seq Assembly with or without the Truth Evaluation
  • DETONATE consists of two component packages, RSEM-EVAL and REF-EVAL. Both packages are mainly intended to be used to evaluate de novo transcriptome assemblies, although REF-EVAL can be used to compare sets of any kinds of genomic sequences.
    
    
• Trinity 2.1.1 – RNA-Seq De novo Assembly
  
  • Trinity represents a novel method for the efficient and robust de novo reconstruction of transcriptomes from RNA-Seq data. Trinity combines three independent software modules: Inchworm, Chrysalis, and Butterfly, applied sequentially to process large volumes of RNA-Seq reads. Trinity partitions the sequence data into many individual de Bruijn graphs, each representing the transcriptional complexity at at a given gene or locus, and then processes each graph independently to extract full-length splicing isoforms and to tease apart transcripts derived from paralogous genes.
    
    
• IsoSCM 2.0.11 – Transcript Assembly tool using Multiple Change-point Inference to improve 3’UTR Annotation
  • IsoSCM (Isoform Structural Change Model) is a new method for transcript assembly  that incorporates change-point analysis to improve the 3′ UTR annotation process.
    
    
• IVA 1.0.3 – Iterative Virus Assembler
  • 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.
    
    
• SFA-SPA 0.2.1 – A Suffix Array based Short Peptide Assembler for Metagenomic Data
  • SFA-SPA is a suffix array based short peptide assembler for metagenomic data
    
    
• RAMPART 0.12.2 – A Workflow Management System for de novo Genome Assembly
  • RAMPART is a de novo assembly pipeline that makes use of third party-tools and High Performance Computing resources. It can be used as a single interface to several popular assemblers, and can perform automated comparison and analysis of any generated assemblies
    
    
• Celera Assembler 8.3 – Whole Genome Shotgun Assembler
  • Celera Assembler (wgs-assembler) is scientific software for DNA research. It can reconstruct long sequences of genomic DNA given the fragmentary data produced by whole-genome shotgun sequencing. The Celera Assembler has enabled discovery in microbial genomes, large eukaryotic genomes, diploid genomes, and genomes from environmental samples. Celera Assembler contributed the first diploid sequence of an individual human, and metagenomics assemblies of the Global Ocean Sampling
    
    
• A5-miseq 20150522 – de novo Assembly & Analysis of Illumina Sequence data
  • de novo assembly & analysis of Illumina sequence data, including the A5 pipeline, A5-miseq, tools to evaluate assembly quality, and scripts to facilitate data submission to NCBI and the RAST annotation system
    
    
• Trans-ABySS 1.5.3 – Analyze ABySS multi-k-assembled Shotgun Transcriptome Data.
  • Trans-ABySS is a software pipeline for analyzing ABySS-assembled contigs from shotgun transcriptome data. The pipeline accepts assemblies that were generated across a wide range of k values in order to address variable transcript expression levels. It first filters and merges the multi-k assemblies, generating a much smaller set of nonredundant contigs. It contains scripts that map assembled contigs to known transcripts, currently supporting Blat and Exonerate contig-to-genome aligners. It identifies novel splicing events like exon-skipping, novel exons, retained introns, novel introns, and alternative splicing. Its scripts can also estimate gene expression levels, identify candidate polyadenylation sites, and identify candidate gene-fusion events.
    
    
• SAT-Assembler 20160120 – Scalable and Accurate Targeted Gene Assembly Tool
  • SAT-Assembler can perform targeted gene assembly for both RNA-Seq and metagenomic data. It addresses the above challenges of de novo assembly of large scale NGS data by conducting family-specic gene assembly, homology-guided overlap graph construction, and careful graph traversal.
    
    
• Opera 2.0.2 – Sequence Assembly Program
  • Opera (Optimal Paired-End Read Assembler) is a sequence assembly program . It uses information from paired-end reads to optimally order and orient contigs assembled from shotgun-sequencing reads.
    
    
• Sequencher 5.4.1 – DNA Sequence Assembly and Analysis
  • Sequencher is the industry standard software for DNA sequence analysis. It works with all automated sequencers and is widely known for its lightning-fast contig assembly, short learning curve, user-friendly editing tools, and superb technical support. First released almost 15 years ago, Sequencher is currently used for sequence analysis tasks in every major genomic and pharmaceutical company as well as numerous academic and government labs in over 40 countries around the world. Life Science researchers use Sequencher for many diverse DNA sequence analysis applications including de novo gene sequencing, mutation detection, forensic human identification, systematics, and more.
    
    
• Minia 2.0.3 – Short-read Assembler based on a de Bruijn graph
  • Minia is a short-read assembler based on a de Bruijn graph, capable of assembling a human genome on a desktop computer in a day
    
    
• MaSuRCA 3.1.3 – Whole Genome Short Read Assembler
  • MaSuRCA is whole genome assembly software. It combines the efficiency of the de Bruijn graph and Overlap-Layout-Consensus (OLC) approaches. MaSuRCA can assemble data sets containing only short reads from Illumina sequencing or a mixture of short reads and long reads (Sanger, 454).
    
    
• KmerGenie 1.6982 – K-mer size Selection for Genome Assembly
  • KmerGenie estimates the best k-mer length for genome de novo assembly. Given a set of reads, KmerGenie first computes the k-mer abundance histogram for many values of k. Then, for each value of k, it predicts the number of distinct genomic k-mers in the dataset, and returns the k-mer length which maximizes this number. Experiments show that KmerGenie’s choices lead to assemblies that are close to the best possible over all k-mer lengths.
    
    
• pilon v1.16 – Automated Assembly Improvement
  • pilon uses read alignment analysis to diagnose, report, and automatically improve de novo genome assemblies.
    
    
• Phred/Phrap/Consed 29.0 – DNA Sequence Assembler & Finishing Tools
  
  • phrap is a program for assembling shotgun DNA sequence data. Among other features, it allows use of the entire read and not just the trimmed high quality part, it uses a combination of user-supplied and internally computed data quality information to improve assembly accuracy in the presence of repeats, it constructs the contig sequence as a mosaic of the highest quality read segments rather than a consensus, it provides extensive assembly information to assist in trouble-shooting assembly problems, and it handles large datasets.
    
    
• CLC Genomics Workbench 8.5.1 – Assembly & Analysis of Sequencing Data
  • CLC Genomics Workbench, for analyzing and visualizing Next Generation Sequencing data, incorporates cutting-edge technology and algorithms, while also supporting and integrating with the rest of your typical NGS workflow.
    
    
• Metassembler 1.5 – Combines multiple Whole Genome de novo Assemblies into a combined Consensus Assembly
  • Metassembler is a software package for reconciling assemblies produced by de novo short-read assemblers such as SOAPdenovo and ALLPATHS-LG. The goal of assembly reconciliation, or “metassembly,” is to combine multiple assemblies into a single genome that is superior to all of its constituents
    
    
• Tablet 1.15.09.01 – Next Generation Sequence Assembly Visualization
  • Tablet is a lightweight, high-performance graphical viewer for next generation sequence assemblies and alignments.Supporting a range of input assembly formats, Tablet provides high-quality visualizations showing data in packed or stacked views, allowing instant access and navigation to any region of interest, and whole contig overviews and data summaries. Tablet is both multi-core aware and memory efficient, allowing it to handle assemblies containing millions of reads, even on a 32-bit desktop machine.
    
    
• ABySS 1.9.0 – de novo, parallel, paired-end Sequence Assembler
  • ABySS (Assembly By Short Sequences) is a de novo, parallel, paired-end sequence assembler that is designed for short reads. The single-processor version is useful for assembling genomes up to 100 Mbases in size. The parallel version is implemented using MPI and is capable of assembling larger genomes.
    
    
• CLEAT 2.0 – Identifies 3′ UTR Ends of Transcripts in de novo RNA-Seq Assemblies
  • CLEAT is a post-processing tool for CLEavage site Analysis of Transcriptomes. CLEAT is designed to work on trans-ABySS output.
    
    
• StriDe – novel Assembler
  • The StriDe Assembler integrates string and de Bruijn graph by decomposing reads within error-prone regions, while extending paire-end read into long reads for assembly through repetitive regions.
    
    
• REAPR 1.0.18 – Genome Assembly Evaluation
  • REAPR (Recognising Errors in Assemblies using Paired Reads) is a tool that evaluates the accuracy of a genome assembly using mapped paired end reads, without the use of a reference genome for comparison. It can be used in any stage of an assembly pipeline to automatically break incorrect scaffolds and flag other errors in an assembly for manual inspection. It reports mis-assemblies and other warnings, and produces a new broken assembly based on the error calls.
    
    
• GapFiller 1.10 – Close Gaps within Pre-assembled Scaffolds
  • GapFiller is a stand-alone program for closing gaps within pre-assembled scaffolds. It is unique in offering the possibility to manually control the gapclosure process. By using the distance information of paired-read data, GapFiller seeks to close the gap from each edge in an iterative manner. From a good number of tests we see the program yields excellent results both on bacterial en eukaryotic  datasets. The command-line Perl script and additional files van be downloaded below. The input data is given by pre-assembled scaffold sequences (FASTA) and NGS paired-read data (FASTA or FASTQ).
    
    
• SSAKE 3.8.4 – Assembling Millions of short DNA Sequences
  • SSAKE is a genomics application for assembling millions of very short DNA sequences.SSAKE is designed to help leverage the information from short sequence reads by stringently assembling them into contiguous sequences that can be used to characterize novel sequencing targets.
    
    
• SGA 0.10.14 – String Graph Assembler
  • SGA is a de novo assembler designed to assemble large genomes from high coverage short read data. The major goal of SGA is to be very memory efficient, which is achieved by using a compressed representation of DNA sequence reads.
    
    
• r2cat – Synteny Plots & Comparative Assembly
  
  • r2cat (related reference based contig arrangement tool) can be used to order a set of contigs with respect to a single reference genome. This is done by mapping the contigs onto the reference using a q-gram filter. The mapping is visualized in a synteny plot.
    
    
• TASR 1.6 – Targeted Assembly of Sequence Reads
  • TASR (Targeted Assembly of Sequence Reads)  is a genomics application that allows hypothesis-based interrogation of genomic regions (sequence targets) of interest.
    
    
• Rainbow v2.0.4 – Clustering and Assembling Short Reads, especially for RAD
  • Rainbow package consists of several programs used for RAD-seq related clustering and de novo assembly.
    
    
• CAFTOOLS 2.0.2 – Tools for the Common Assembly Format (CAF)
  • CAFTOOLS comprises a set of libraries and programs for manipulating DNA sequence assemblies using CAF files, a comprehensive representation of a sequence assembly as a text file.
• Gap Resolution – Improving Newbler Genome Assemblies. Gap Resolution was developed by DOE Joint Genome Institute to improve Newbler genome assemblies by automating the closure of sequence gaps caused by repetitive regions in the DNA.
  
  
• Meraculous 2.0.5 – De novo Genome Assembler from Short Reads
  • Meraculous is a new algorithm for whole genome assembly of deep paired-end short reads, and apply it to the assembly of a dataset of paired 75-bp Illumina reads derived from the 15.4 megabase genome of the haploid yeast Pichia stipitis.
    
    
• COPE 1.2.5 – Pair-end Reads Connection tool to facilitate Genome Assembly
  • COPE (Connecting Overlapped Pair-End reads) is a method to align and connect the illumina sequenced Pair-End reads of which the insert size is smaller than the sum of the two read length.The connected reads can be used in genome assembly, resequencing and transcriptome research.
    
    
• PEAR 0.9.6 – Pair-End reads AssembleR
  • PEAR is an ultrafast, memory-efficient and highly accurate pair-end reads assembler. It is fully parallelized and can run with as low as just a few kilobytes of memory.
    
    
• EBARDenovo 2.0.1 – Highly-accurate de novo Assembler of Paired-end RNA-Seq
  • EBARDenovo is a highly-accurate search-based de novo assembler of paired-end RNA-Seq for advance transcriptomic study.
    
    
• EagleView 2.2 – Genome Assembler Viewer
  • EagleView is an information-rich genome assembler viewer with data integration capability. EagleView can display a dozen different types of information including base qualities, machine specific trace signals, and genome feature annotations. It provides an easy way for inspecting visually the quality of a genome assembly and validating polymorphism candidate sites (e.g., SNPs) reported by polymorphism discovery tools. It can also facilitate data interpretation and hypothesis generation.
    
    
• MAIA 0.5 – Integrating Genome Assemblies

  • MAIA (Multiple Assembly IntegrAtion) is an algorithm to integrate multiple genome assemblies. For example, assemblies originating from:
    – Different runs of a de novo assembler
    – Assemblies of different data types
    – Comparative assemblies

• InteMAP 1.0 – Integrated Metagenomic Assembly pipeline for NGS Short Reads
  
  • InteMAP is a pipeline which integrates individual assemblers for assembling metagenomic short sequencing reads.
    
    
• MAP 20121108 – A de novo Metagenomic Assembly program for Shotgun DNA reads
  • MAP (Metagenomic Assembly program) is a de novo assembly approach and its implementation based on an improved Overlap/Layout/Consensus (OLC) strategy incorporated with several special algorithms.MAP uses the mate pair information, resulting in being more applicable to shotgun DNA reads (recommended as > 200 bp) currently widely-used in metagenome projects. Results of extensive tests on simulated data show that MAP can be superior to both Celera and Phrap for typical longer reads by Sanger sequencing, as well as has an evident advantage over Celera, Newbler, and the newest Genovo, for typical shorter reads by 454 sequencing.
    
    
• Phusion 2.1c – Assembly Genome Sequences from Whole Genome Shotgun(WGS) Reads
  • Phusion is a software package for assembling genome sequences from whole genome shotgun(WGS) reads.
    
    
• CodonCode Aligner 6.0.2 – DNA Sequence Assembly & Alignment
  • CodonCode Aligner is a program for sequence assembly, contig editing, and mutation detection, available for Windows and Mac OS X. Aligner is compatible with Phred-Phrap and fully supports sequence quality scores, while offering a familiar, easy-to-learn user interface.
    
    
• Cerulean 0.1.1 – Hybrid Genome Assembler
  • Cerulean is a hybrid assembly using high throughput short and long reads
    
    
• Ragout 1.2 – Tool for Reference-assisted Assembly
  • Ragout (Reference-Assisted Genome Ordering UTility) is a tool for assisted assembly using multiple references. It takes a short read assembly (a set of contigs), a set of related references and a corresponding phylogenetic tree and then assembles the contigs into scaffolds.
    
    
• laSV 1.0.2 – Local Assembly based Structural Variation Discovery tool
  • laSV is a software that employs a local de novo assembly based approach to detect genomic structural variations from whole-genome high-throughput sequencing datasets.
    
    
• SPAdes 3.6.2 – Single-cell Genome Assembler
  • SPAdes (St. Petersburg genome assembler) is intended for both standard isolates and single-cell MDA bacteria assemblies.
    
    
• PERGA 0.5.03.02 – Paired End Reads Guided Assembler
  • PERGA is a novel sequence reads guided de novo assembly approach which adopts greedy-like prediction strategy for assembling reads to contigs and scaffolds.
    
    
• Telescoper 0.2 – De novo Assembly Algorithm
  • Telescoper is a local assembly algorithm designed for short-reads from NGS platforms such as Illumina. The reads must come from two libraries: one short insert, and one long insert.
    
    
• MetaCompass 1.0 – Comparative Assembly of Metagenomic Sequences
  • MetaCompass is a software package for comparative assembly of metagenomic reads. MetaCompass achieves comparable assembly performance to the state of the art de novo assemblers, but these two different approaches complement each other a lot. So combining contigs between MetaCompass and other independent de novo assemblers give us the best overall metagenomic assembly.
    
    
• SCARF – Scaffolded and Corrected Assembly of Roche 454
  • SCARF is a next-gen sequence assembly tool for evolutionary genomics. Designed especially for assembling 454 EST sequences against high quality reference sequences from related species.
    
    
• MetaCAA – Assembly of Metagenomic Datasets
  • MetaCAA is a sequence-assembly tool specifically intended for metagenomes.
    
    
• Contiguity 1.0.4 – Contig Adjacency Graph Construction and Visualisation
  • Contiguity is interactive software for the visualization and manipulation of de novo genome assemblies.
    
    
• ScaffoldScaffolder 0.1 – Solving Contig Orientation via Bidirected to Directed Graph Reduction
  • ScaffoldScaffolder is a stand-alone scaffolding algorithm which was designed specifically for scaffolding diploid genomes.
    
    
• HaploClique 0.1 – Viral Quasispecies Assembly from Paired-end data
  • HaploClique is a computational approach to reconstruct the structure of a viral quasispecies from next-generation sequencing data as obtained from bulk sequencing of mixed virus samples.
    
    
• TAG 0.91 – Transcript Assembly by Mapping Reads to Graphs
  • TAG is a tool for metatranscriptome assembly using de Bruijn graph of matched metagenome as the reference
    
    
• EPGA2 – De Novo Assembler
  • EPGA2 updates some modules in EPGA which can improve memory efficiency in genome asssembly.
    
    
• GMcloser 1.5.1 / GMvalue 1.3 – Closing the Gaps in Scaffolds with Preassembled Contigs
  • GMcloser fills and closes the gaps present in scaffold assemblies, especially those generated by the de novo assembly of whole genomes with next-generation sequencing (NGS) reads.
    
    
• SLICEMBLER – Meta-assembler Designed for Ultra-deep Sequencing data
  • SLICEMBLER is a meta-assembler designed for ultra-deep sequencing data
    
    
• SEQLandscape v1 – Generation and Visualization of Sequence Landscape

  • SEQLandscape is an application allowing the generation and visualization of a sequence landscape. HyDA-Vista: Towards Optimal Guided Selection of k-mer Size for Sequence Assembly.

• misSEQuel v1.0beta – Misassembly Detection in Draft Genomes
  • misSEQuel is a software that enhances the quality of draft genomes by identifying misassembly errors and their breakpoints using paired-end sequence reads and optical mapping data.
    
    
• Dawg 1.2 – Simulating Sequence Evolution
  • Dawg (DNA Assembly with Gaps) is an application designed to simulate the evolution of recombinant DNA sequences in continuous time based on the robust general time reversible model with gamma and invariant rate heterogeneity and a novel length-dependent model of gap formation.
    
    
• BUSCO v1.1b1 – Assessing Genome Assembly and Annotation Completeness with Single-copy Orthologs
  • BUSCO completeness assessment employs sets of Benchmarking Universal Single-Copy Orthologs from OrthoDB to provide quantitative measures of the completeness of genome assemblies, annotated gene sets, and transcriptomes in terms of expected gene content.
    
    
• FinisherSC 2.0 – A Repeat-aware tool for upgrading de-novo Assembly using Long Reads
  • FinisherSC is a repeat-aware and scalable tool for upgrading de-novo assembly using long reads.
    
    
• WhatsHap – Haplotype Assembly for Future-Generation Sequencing Reads
  • WhatsHap is a software for phasing genomic variants using DNA sequencing reads, also called haplotype assembly. It is especially suitable for long reads, but works also well with short reads.
    
    
• Compartmentalized Assembler – Assembly of Physical Maps
  • Compartmentalized assembler is a novel method for the assemlby of high quality physical maps from fingerprinted clones.
    
    
• Elviz – Exploration of Metagenomic Assemblies
  • Elviz (Environmental Laboratory Visualization) is an interactive web-based tool for the visual exploration of assembled metagenome data and their complex metadata.
    
    
• SSP – de novo Transcriptome Assembler
  • SSP is a de novo transcriptome assembler that assembles RNA-seq reads into transcripts. SSP aims to reconstructs all the alternatively spliced isoforms and estimates the expression level of them.
    
    
• VirAmp – Galaxy-based Viral Genome Assembly pipeline
  • VirAmp is a web-based semi-de novo fast virus genome assembly pipeline designed for extremely high coverage NGS data. VirAmp is a collection of existing tools, combined into a single Galaxy interface. Users without further computational knowledge can easily operate the pipeline.
    
    
• aTRAM 1.04 – automated Target Restricted Assembly Method
  • aTRAM performs targeted de novo assembly of loci from paired-end Illumina runs.
    
    
• Ray 2.3.1 – Parallel Genome Assemblies for Parallel DNA sequencing
  • Ray is a parallel software that computes de novo genome assemblies with next-generation sequencing data.
    
    
• CAR – Contig Assembly of Prokaryotic Draft Genomes Using Rearrangements
  • CAR is an efficient and more accurate tool for assembling contigs of a prokaryotic draft genome based on a reference genome.
    
    
• VTBuilder – Assembly of Multi Isoform Transcriptomes
  • VTBuilder is a tool for the inference of non-chimeric contigs from read data that has been sequenced from complex multi-isoformic transcriptomes, such as snake venom glands, or rapidly evolving viral populations, such as HIV-1.
    
    
• TruHmm – TRanscription Unit Assembly by a Hidden Markov model
  • TruHmm is a reference based transcriptome assembler for prokaryotes, and is suitable for assembling transcripts for directional RNA-seq library.
    
    
• Bridger 20141201 – RNA-Seq Assembly
  • Bridger is a new de novo transcriptome assembler which takes advantage of techniques employed in Cufflinks to overcome limitations of the existing de novo assemblers.
    
    
• GRASP 0.0.4 – Guided Reference-based Assembly of Short Peptides
  • GRASP is a gene annotation tool for metagenomic studies. GRASP assembles the fragmented short-peptides, which are called from the NGS reads, and aligns the assembled contigs to the query reference protein. GRASP achieves much higher sensitivity than BLASTP for gene annotation purpose.
    
    
• Cortex 1.05.21 – Genome Assembly and Variation Analysis
  • Cortex is an efficient and low-memory software framework for analysis of genomes using sequence data. There are two main executables, being developed in parallel streams: cortex_con (primary contact Mario Caccamo) is for consensus genome assembly, and cortex_var (primary contact Zamin Iqbal) is for variation and population assembly.
    
    
• MEGAHIT v0.1.4 – Large and Complex Metagenomics Assembly via Succinct de Bruijn graph
  • MEGAHIT is a single node assembler for large and complex metagenomics NGS reads, such as soil. It makes use of succinct de Bruijn graph to achieve low memory usage, whereas its goal is not to make memory usage as low as possible.
    
    
• CISA 20140304 – Contig Integrator for Sequence Assembly
  • CISA has been developed to integrate the assemblies into a hybrid set of contigs, resulting in assemblies of superior contiguity and accuracy, compared with the assemblies generated by the state-of-the-art assemblers and the hybrid assemblies merged by existing tools
    
    
• Cufflinks 2.2.1 – Transcript Assembler & Abundance Estimator for RNA-Seq
  • Cufflinks assembles transcripts, estimates their abundances, and tests for differential expression and regulation 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.
    
    
• mapsembler 2.2.4 – Targetted Assembly of Short Sequence Reads
  • Mapsembler is a targeted assembly software. It takes as input a set of NGS raw reads and a set of input sequences (starters). It first determines if each starter is read-coherent, e.g. whether reads confirm the presence of each starter in the original sequence. Then for each read-coherent starter, Mapsembler outputs its sequence neighborhood as a linear sequence or as a graph, depending on the user choice.
    
    
• Tedna 1.2.2 – Transposable Element De Novo Assembler
  • Tedna is a lightweight de novo transposable element assembler. It assembles the transposable elements directly from the raw reads.
    
    
• HyDA 1.3.1 / Squeezambler 2.0.3 – Hybrid De Novo Assembler
  • HyDA is a multipurpose assembler, particularly tested for single cell and normal multicell genome co-assembly
    
    
• PANDASEQ 2.8 / Pandaseq-sam 1.3 – PAired-eND Assembler for DNA sequences
  • PANDASEQ is a program to align Illumina reads, optionally with PCR primers embedded in the sequence, and reconstruct an overlapping sequence.
    
    
• ZORRO 2.2 – Hybrid Sequencing Technology Assembler
  • ZORRO is a hybrid sequencing technology assembler. It merges two sets of pre-assembled contigs into a more contiguous and consistent assembly.
    
    
• FLASH 1.2.11 – Fast Length Adjustment of SHort reads
  • FLASH (Fast Length Adjustment of SHort reads) is a very accurate fast tool to merge paired-end reads from fragments that are shorter than twice the length of reads. The extended length of reads has a significant positive impact on improvement of genome assemblies.
    
    
• ALLPATHS-LG 51750 – Whole Genome Shotgun Assembler
  • ALLPATHS-LG (Large Genome) is a whole genome shotgun assembler that can generate high quality assemblies from short reads. It works on both small and large (mammalian size) genomes. To use it, you should first generate ~100 base Illumina reads from two libraries: one from ~180 bp fragments, and one from ~3000 bp fragments, both at about 45x coverage. Sequence from longer fragments will enable longer-range continuity.
    
    
• More Tools at http://bioinformaticsonline.com/pages/view/30440/genome-assembly-tools-and-software-part2

With advances in Cancer Genomics, Mutation Annotation Format (MAF) is being widley accepted and used to store variants detected. The Cancer Genome Atlas Project has seqenced over 30 different cancers with sample size of each cancer type being over 200. The resulting data consisting of genetic variants is stored in the form of Mutation Annotation Format. This package attempts to summarize, analyze, annotate and visualize MAF files in an efficient manner either from TCGA sources or any in-house studies as long as the data is in MAF format. Maftools can also handle ICGC Simple Somatic Mutation format.

maftools is on bioRxiv

Please cite the below if you find this tool useful for you.

Mayakonda, A. and H.P. Koeffler, Maftools: Efficient analysis, visualization and summarization of MAF files from large-scale cohort based cancer studies. bioRxiv, 2016. doi: http://dx.doi.org/10.1101/052662

Address of the bookmark: https://github.com/PoisonAlien/maftools

Prokka – Standalone command line tool, takes just a few minutes per genome. This is the best way to get good quality annotation in a flash, which is particularly useful if you have loads of genomes or need to annotate a pangenome or metagenome. Note however that the quality of functional information is not as good as RAST, and you will need several extra steps if you want to do functional profiling and pathway analysis of your genome(s)… which is in-built in RAST.

NCBI Prokaryotic Genome Annotation Pipeline is designed to annotate bacterial and archaeal genomes (chromosomes and plasmids).

Genome annotation is a multi-level process that includes prediction of protein-coding genes, as well as other functional genome units such as structural RNAs, tRNAs, small RNAs, pseudogenes, control regions, direct and inverted repeats, insertion sequences, transposons and other mobile elements.

PGAP: NCBI has developed an automatic prokaryotic genome annotation pipeline that combines ab initio gene prediction algorithms with homology based methods. The first version of NCBI Prokaryotic Genome Automatic Annotation Pipeline (PGAAP; see Pubmed Article) developed in 2005 has been replaced with an upgraded version that is capable of processing a larger data volume.  NCBI's annotation pipeline depends on several internal databases and is not currently available for download or use outside of the NCBI environment.

BEACON (automated tool for Bacterial GEnome Annotation ComparisON), a fast tool for an automated and a systematic comparison of different annotations of single genomes. The extended annotation assigns putative functions to many genes with unknown functions. BEACON is available under GNU General Public License version 3.0 and is accessible at: http://www.cbrc.kaust.edu.sa/BEACON/.

BlastKOLA: Assigns K numbers to the user's sequence data by BLAST searches, respectively, against a nonredundant set of KEGG GENES. KOALA (KEGG Orthology And Links Annotation) is KEGG's internal annotation tool for K number assignment of KEGG GENES using SSEARCH computation. Annotate Sequence in KEGG Mapper and Pathogen Checker in KEGG Pathogen are special interfaces to this server and can be executed in an interactive mode. BlastKOALA is suitable for annotating fully sequenced genomes.

PAGIT: Provides a toolkit for improving the quality of genome assemblies created via an assembly software. PAGIT compiled four tools: (i) ABACAS which classifies and orientates contigs and estimates the sizes of gaps between them; (ii) IMAGE uses paired-end reads to extend contigs and close gaps within the scaffolds; (iii) ICORN for identifying and correcting small errors in consensus sequences and; (iv) RATT for help annotation. The software was mainly created to analyze parasite genomes of up to about 300 Mb.

MAKER: A portable and easily configurable genome annotation pipeline. MAKER allows smaller eukaryotic and prokaryotic genome projects to independently annotate their genomes and to create genome databases. It 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. MAKER's inputs are minimal and its ouputs can be directly loaded into a Generic Model Organism Database (GMOD). They can also be viewed in the Apollo genome browser; this feature of MAKER provides an easy means to annotate, view and edit individual contigs and BACs without the overhead of a database. MAKER is available for download and can be tested online via the MAKER Web Annotation Service (MWAS).

MyPro is a software pipeline for high-quality prokaryotic genome assembly and annotation. It was validated on 18 oral streptococcal strains to produce submission-ready, annotated draft genomes. MyPro installed as a virtual machine and supported by updated databases will enable biologists to perform quality prokaryotic genome assembly and annotation with ease.

Following are the list of tools commonly used to PPIs predictions:

Protein-Protein Interaction Sites

PPISP

A consensus neural network method for predicting protein-protein interaction sites

HOMCOS

A server to predict interacting protein pairs and interacting sites by homology modeling of complex structures

HotPOINT

Prediction of protein interfaces using an empirical model

ISIS

Prediction of interaction hotspots from sequence

KFC server

Automated decision-tree approach to predicting protein-protein interaction hot spots

meta-PPISP

A meta server for predicting protein-protein interaction sites. meta-PPISP is built on three individual web servers: cons-PPISP, PINUP, and Promate

ODA

Identification of optimal surface patches with the lowest docking desolvation energy values

PINUP

Protein binding site prediction with an empirical scoring function

Other Sites (DNA, RNA, Metals)

CHED 

Web server for predicting soft metal binding sites in proteins

DBD-Hunter

A knowledge-based method for the prediction of DNA-protein interactions

DISPLAR

Given the structure of a protein known to bind DNA, the method predicts residues that contact DNA using neural network method

iDBPs

Predicts DNA binding proteins for proteins with known 3D structure.

PFplus

Below you can find a few basic tutorials for how to run MITObim and I encorage you to give them a try with the testdata that comes with this Repo, just to make sure everything is running smoothly on your system. It'll only take a few minutes and will potentially safe you a lot of time down the line.

I provide further examples here as Jupyter notebooks. Get in touch if you feel like sharing your particular MITObim solution and I'd be happy to put it up here, too!

Address of the bookmark: https://github.com/chrishah/MITObim

https://www.nature.com/articles/nmeth.4432

Address of the bookmark: https://github.com/xiaochuanle/MECAT