• paired-end (PE) and/or mate-pair libraries (NGS-based mode)
• long reads (NGS-based mode)
• synteny to the genome of some related species (reference-based mode)

Scaffolding 

In reference-based mode, pyScaf uses synteny to the genome of closely related species in order to order contigs and estimate distances between adjacent contigs.

Contigs are aligned globally (end-to-end) onto reference chromosomes, ignoring:

• matches not satisfying cut-offs (--identity and --overlap)
• suboptimal matches (only best match of each query to reference is kept)
• and removing overlapping matches on reference.

In preliminary tests, pyScaf performed superbly on simulated heterozygous genomes based on C. parapsilosis (13 Mb; CANPA) and A. thaliana (119 Mb; ARATH) chromosomes, reconstructing correctly all chromosomes always for CANPA and nearly always for ARATH (Figures in dropbox, CANPA table, ARATH table).
Runs took ~0.5 min for CANPA on 4 CPUs and ~2 min for ARATH on 16 CPUs.

Important remarks:

• Reduce your assembly before (fasta2homozygous.py) as any redundancy will likely break the synteny.
• pyScaf works better with contigs than scaffolds, as scaffolds are often affected by mis-assemblies (no de novo assembler / scaffolder is perfect...), which breaks synteny.
• pyScaf works very well if divergence between reference genome and assembled contigs is below 20% at nucleotide level.
• pyScaf deals with large rearrangements ie. deletions, insertion, inversions, translocations. Note however, this is experimental implementation!
• Consider closing gaps after scaffolding.

Address of the bookmark: https://github.com/lpryszcz/pyScaf

Address of the bookmark: https://github.com/lh3/fermi

More at http://biorxiv.org/content/early/2016/09/21/076463 

Address of the bookmark: http://shinyheatmap.com/

Targeted bait capture is a technique for sequencing many loci simultaneously based on bait sequences. HybPiper pipeline starts with high-throughput sequencing reads (for example from Illumina MiSeq), and assigns them to target genes using BLASTx or BWA. The reads are distributed to separate directories, where they are assembled separately using SPAdes. The main output is a FASTA file of the (in frame) CDS portion of the sample for each target region, and a separate file with the translated protein sequence.

HybPiper also includes post-processing scripts, run after the main pipeline, to also extract the intronic regions flanking each exon, investigate putative paralogs, and calculate sequencing depth. For more information, please see our wiki.

HybPiper is run separately for each sample (single or paired-end sequence reads). When HybPiper generates sequence files from the reads, it does so in a standardized directory hierarchy. Many of the post-processing scripts rely on this directory hierarchy, so do not modify it after running the initial pipeline. It is a good idea to run the pipeline for each sample from the same directory. You will end up with one directory per run of HybPiper, and some of the later scripts take advantage of this predictable directory structure.

Address of the bookmark: https://github.com/mossmatters/HybPiper

More at

https://github.com/jts/sga

SGA dependencies:
-google sparse hash library (http://code.google.com/p/google-sparsehash/)
-the bamtools library (https://github.com/pezmaster31/bamtools)
-zlib (http://www.zlib.net/)
-(optional but suggested) the jemalloc memory allocator (http://www.canonware.com/jemalloc/download.html)

Address of the bookmark: https://github.com/jts/sga

Cutadapt helps with these trimming tasks by finding the adapter or primer sequences in an error-tolerant way. It can also modify and filter reads in various ways. Adapter sequences can contain IUPAC wildcard characters. Also, paired-end reads and even colorspace data is supported. If you want, you can also just demultiplex your input data, without removing adapter sequences at all.

Cutadapt comes with an extensive suite of automated tests and is available under the terms of the MIT license.

If you use cutadapt, please cite DOI:10.14806/ej.17.1.200 .

More at https://github.com/marcelm/cutadapt

Address of the bookmark: http://cutadapt.readthedocs.io/en/stable/guide.html

MeGAMerge (A tool to merge assembled contigs, long reads from metagenomic sequencing runs)

Description

MeGAMerge is a perl based wrapper/tool that can accept any number of sequence (FASTA) files containing assembled contigs of any length in Multi-FASTA format to produce an improved contig set based on OLC based assembly. All overlap parameters (Minimum Overlap Length, Identity, etc) are user-declarable at runtime. It is written to run on Linux.

Requirements:

You will need to have the following tools installed and in $PATH, or added to $binpath in the tool:

Newbler (specifically runAssembly)
Minimus2 (part of AMOS, also requires MUMmer)

Address of the bookmark: https://github.com/LANL-Bioinformatics/MeGAMerge

Bioinformatics supports to researchers
Customized training in Bioinformatics for researchers and faculty members
Support in Installing, implementing and maintaining software on computer.
Create awareness for taking preventive measure against data security
Organize workshops on thrust ares of Bioinformatics
Research Training to students of Biotechnology and Bioinformatics

More at http://bioinfo.bisr.res.in/index.php

Rajiv Gandhi Centre for Biotechnology (RGCB) invites applications for the following three faculty scientist positions:

Scientist E-II or F in Bioinformatics & Computational Biology

SCIENTIST E-II OR F IN COMPUTATIONAL BIOLOGY & BIOINFORMATICS

Highly motivated and innovative individual who will pursue basic research, solve biological problems with emphasis on computational and quantitative experimental methods and build active bridges to translational research. The scientist will also provide computational biology support to ongoing research programs in disease biology, provide assistance to analyze complex data sets generated by RGCB scientists and collaborators inclusive of including high dimensional “omics” data and next generation sequencing data, such as whole genome, exome, RNA-seq and ChIP-seq as well as provide leadership for high quality training for junior scientists and regular teaching programs of the institute. Areas of research of interest to RGCB include but are not limited to computational, systems, or quantitative biology with applications to cell biology, developmental biology, metabolism, genomics, proteomics, biophysics, biological information systems, network pharmacology, drug design and cancer research. The scientist’s responsibilities include efforts for the integration of DNA variant annotation with statistical genetic analysis methods including linkage, imputation and association methods, adopting novel and innovative methodologies to analyze, integrate and interpret high dimensional data sets, provision of annotation to robust genetics and genomics findings using several data sources and methods, data management of exploratory clinical and R&D studies in partnership with other lines of genetic data generated from internal and external studies, delivery and documentation of genomic information to support genetic studies, ensuring high-quality genetic and genomic data is incorporated into exploratory- clinical research programs, developing tools that make maximum use of multiple data sources to support annotation of DNA variation and contributes to systems biology initiatives within RGCB

More at http://rgcb.res.in/scientist-positions/

Application Form http://rgcb.res.in/wp-content/uploads/2015/11/APPLICATION-FORMAT-FOR-SCIENTISTS.docx