Reference

http://www.sciencedaily.com/releases/2013/12/131204132018.htm

Following are the genome assembly tools based on string graph:

1.SGA (String Graph Assembler) https://github.com/jts/sga

Assembles large genomes from high coverage short read data. SGA is designed as a modular set of programs, which are used to form an assembly pipeline. SGA implements a set of assembly algorithms based on the FM-index. As the FM-index is a compressed data structure, the algorithms are very memory efficient. The SGA assembly has three distinct phases. The first phase corrects base calling errors in the reads. The second phase assembles contigs from the corrected reads. The third phase uses paired end and/or mate pair data to build scaffolds from the contigs. The output of this software is a PDF report that allows the properties of the genome and data quality to be visually explored. By providing more information to the user at the start of an assembly project, this software will help increase awareness of the factors that make a given assembly easy or difficult, assist in the selection of software and parameters and help to troubleshoot an assembly if it runs into problems.

2. SAGE: String-overlap Assembly of GEnomes https://github.com/lucian-ilie/SAGE2

SAGE, for de novo genome assembly. As opposed to most assemblers, which are de Bruijn graph based, SAGE uses the string-overlap graph. SAGE builds upon great existing work on string-overlap graph and maximum likelihood assembly, bringing an important number of new ideas, such as the efficient computation of the transitive reduction of the string overlap graph, the use of (generalized) edge multiplicity statistics for more accurate estimation of read copy counts, and the improved use of mate pairs and min-cost flow for supporting edge merging. The assemblies produced by SAGE for several short and medium-size genomes compared favourably with those of existing leading assemblers.

3. FSG: Fast String Graph

The new integrated assembler has been assessed on a standard benchmark, showing that fast string graph (FSG) is significantly faster than SGA while maintaining a moderate use of main memory, and showing practical advantages in running FSG on multiple threads. Moreover, we have studied the effect of coverage rates on the running times.

4.  BASE https://github.com/dhlbh/BASE

It enhances the classic seed-extension approach by indexing the reads efficiently to generate adaptive seeds that have high probability to appear uniquely in the genome. Such seeds form the basis for BASE to build extension trees and then to use reverse validation to remove the branches based on read coverage and paired-end information, resulting in high-quality consensus sequences of reads sharing the seeds. Such consensus sequences are then extended to contigs. BASE is a practically efficient tool for constructing contig, with significant improvement in quality for long NGS reads. It is relatively easy to extend BASE to include scaffolding.

5. Fermi https://github.com/lh3/fermi/

Fermi is a de novo assembler with a particular focus on assembling Illumina short sequence reads from a mammal-sized genome. In addition to the role of a typical assembler, fermi also aims to preserve heterozygotes which are often collapsed by other assemblers. Its ultimate goal is to find a minimal set of unitigs to represent all the information in raw reads.

If you want to learn about String Graph assembler, please read the following papers -

i) The Fragment Assembly String Graph - E. W. Myers

This paper describes the String Graph concept.

ii) Efficient construction of an assembly string graph using the FM-index - Jared T. Simpson and Richard Durbin

This earlier paper from Simpson and Durbin

iii) Efficient de novo assembly of large genomes using compressed data structures - Jared T. Simpson and Richard Durbin

(Under Tenure-Track Model)

Open Advt. No IIITA/DIC/16012014

IIIT-Allahabad has several Openings for the Faculty positions at the Assistant Professor level.

It is a regular tenure-track faculty positions for 3-5 years in teaching and research. A regular faculty is expected to engage heavily in research and teaching. The eligibility criteria for regular faculty positions are similar as in IITs. For an Assistant Professor position, a candidate must have a PhD (in IT/Computer Science &/or Engineering/Electrical, Electronics &/or Communication Engineering/etc; for interdisciplinary areas the PhD may be in an appropriate field), plus three years experience. However, for PhDs from a well known University/Institute (e.g. IITs/IISc/TIFR/ISI in India or well known research universities across the world), and a good research/academic record, the 3 years experience requirement may be waived.

The pay scale for faculty is same as in IITs. Other benefits include initiation research grant, travel support, book grant, professional society membership, etc., and personal benefits such as medical/LTC, on campus subsidized family housing with excellent modern infrastructural facilities.

Areas of Interest

IIIT-Allahabad aims to build strong research groups in important and emerging areas in CS/IT/ECE as well as in emerging interdisciplinary areas, and applications are invited in all these areas. Some of the areas of special interest, besides strengthening the existing research areas, are : Software Engineering, Theoretical Computer Science, Cyber Physical Systems, Robotics, Network science, Digital Media, Computational neuroscience, Machine learning, Healthcare informatics, Computational Biology, Communications networks (both at hardware and protocol levels), Circuits (including VLSI, analog, low power, etc), Energy systems and technologies, Biomedical electronics and systems, Computer Architecture, signal/image processing, Embedded and control systems.

Application Process

Interested candidates can apply by sending their detailed CV with list of publications clearly mentioning Journal names and citation index with three references through email entitled “Faculty positions at IIIT Allahabad” to faculty.applications@iiita.ac.in. Do not send your applications in any other email addresses. Applications will be considered regularly, hence there is no deadline for applying.

Important Clarifications on Eligibility

A PhD in CS/IT (or other disciplines, as announced) is the minimum expected requirement for an Assistant Professor.

Advertisement: http://iiita.ac.in/pub/Faculty-Position-IIITA1.pdf

if you map reads to GRCh38 or hg38, use the following:

ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCA/000/001/405/GCA_000001405.15_GRCh38/seqs_for_alignment_pipelines.ucsc_ids/GCA_000001405.15_GRCh38_no_alt_analysis_set.fna.gz

There are several other versions of GRCh37/GRCh38. What’s wrong with them? Here are a collection of potential issues:

More at http://lh3.github.io/2017/11/13/which-human-reference-genome-to-use

Address of the bookmark: http://lh3.github.io/2017/11/13/which-human-reference-genome-to-use

An Autonomous National Institute under Government of India,
Ministry of Science & Technology
Department of Biotechnology

No: RGCB/ Advt./2014/1
January 24, 2014

Scientist Positions

Group Leader in Computational Biology/Bioinformatics
A 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 analyze complex data sets generated by RGCB scientists and collaborators.

Location: Thiruvananthapuram (Trivandrum)

The above positions will be at the E-II, F or equivalent levels. For senior applicants with an outstanding track record, an option of a contract career path for research excellence at Scientist G or H equivalent level can also be discussed. All positions will initially be for 5 years. Essential and desired qualifications as well as other relevant details for all the above positions are posted on the RGCB website (http://www.rgcb.res.in). The last date for receiving applications is March 14, 2014.

Sd/-
Director

Rajiv Gandhi Centre for Biotechnology
Thycaud, P.O., Poojappura,
Thiruvananthapuram, Kerala, India-695 014
Ph.: 91-471-2529400 (30 Lines), 2347975, 2348104, 2348753, 2345899
Fax: 91-471-2348096, 2346333

More at http://rgcb.res.in/jobs.html

The tool should be compatible with most UNIX flavors and has been successfully tested on the following operating systems:

• Mac OS X 10.11.1
• Mac OS X 10.10.3
• Ubuntu 14.04 LTS

Address of the bookmark: https://github.com/mculinovic/EAGLER

The C-DAC will launch Param Bio Blaze, a supercomputing facility, to address the challenges in bioinformatics on Tuesday at a three-day symposium, titled: 'Accelerating biology: Computing life'. The supercomputing facility will be inaugurated on Tuesday by Ramakrishna Ramaswamy, vice-chancellor, Central University of Hyderabad at the Yashada. The new C-DAC's facility will have a capacity of 10 teraflop and will be able to analyse human cells and its functions.

Param Bio Blaze will help have a larger storage space and better computing facility for the bioinformatics sector. The facility will help capture the movement of molecules and also interaction between two molecules and the effects.

Applications of Param BioBlaze

- Collaboration with National Centre for Cell Science for research on Malaria and understanding how the disease spreads

- Collaborative work with Tata Memorial hospital on breast cancer and find out the difference between normal tissues and tissues from breast cancer patients

- Designing anti-cancer molecules, a collaborative research with the University of Pune

Reference:

Times of India

Image:datacenterjournal.com

I'm running "OPERA-LG_v2.0.5/bin/preprocess_reads.pl" and have the following error:

fail to open file './temporarySam'

[bwa_aln_core] write to the disk... 0.09 sec
[bwa_aln_core] 70778880 sequences have been processed.
[bwa_aln_core] calculate SA coordinate... 161.35 sec
[bwa_aln_core] write to the disk... 0.06 sec
[bwa_aln_core] 70989574 sequences have been processed.
[main] Version: 0.7.15-r1140
[main] CMD: bwa aln -t 30 all_p_ctg.fa -
[main] Real time: 2402.523 sec; CPU: 53429.488 sec
[E::hts_open_format] Failed to open file temporarySam
samtools sort: can't open "temporarySam": No such file or directory
[bwa_aln_core] convert to sequence coordinate... 1.00 sec
[bwa_aln_core] refine gapped alignments... 6.07 sec
[bwa_aln_core] print alignments... PREPROCESS:
Fastq format is recognized
[Thu Jun 14 18:16:47 2018] Building bwa index...
bwa index -p all_p_ctg.fa /home/urbe/Tools/OPERA-LG_v2.0.6/all_p_ctg.fa
[Thu Jun 14 18:18:35 2018] Finding the SA coordinates of the reads using BWA aln...
[Thu Jun 14 18:58:37 2018] Generate alignments of reads using bwa sampe...
bwa samse -n 1 all_p_ctg.fa read.sai - | grep '\(^@\|XT:A:U\)' | /usr/local/bin/samtools view -S -h -b -F 0x4 - | /usr/local/bin/samtools sort -@ 20 -no - temporarySam > FALCON-Unzip-Scaff.bam
Mapping long-reads using blasr...
/home/urbe/Tools/SSpace/SSPACE-LongRead_v1-1/blasr -nproc 40 -m 1 -minMatch 5 -bestn 10 -noSplitSubreads -advanceExactMatches 1 -nCandidates 1 -maxAnchorsPerPosition 1 -sdpTupleSize 7 /media/urbe/MyDDrive/ONTdata/allONT/allONT.fasta /home/urbe/Tools/OPERA-LG_v2.0.6/all_p_ctg.fa | cut -d ' ' -f1-5,7-12 | sed 's/ /\t/g' > FALCON-Unzip-Scaff.map
sh: 1: /home/urbe/Tools/SSpace/SSPACE-LongRead_v1-1/blasr: Permission denied
Sorting mapping results...
sort -k1,1 -k9,9g FALCON-Unzip-Scaff.map > FALCON-Unzip-Scaff.map.sort
Analyzing sorted results...
Extracting linking information...
i3 2000 5000
i2 1000 2000
i4 5000 15000
i0 -200 300
i5 15000 40000
i1 300 1000
Repeat detection...
/home/urbe/Tools/OPERA-LG_v2.0.6/bin//filter_conflicting_edge.pl pairedEdges_i0 contig_length.dat 100 2
Illegal division by zero at /home/urbe/Tools/OPERA-LG_v2.0.6/bin//filter_conflicting_edge.pl line 93.
readline() on closed filehandle FILE at bin/OPERA-long-read.pl line 250.
rm anchor_contig_info.dat contig_length.dat filtered_edges.dat filtered_edges_cov.dat *.sai
rm: cannot remove 'anchor_contig_info.dat': No such file or directory
mv FALCON-Unzip-Scaff.bam FALCON-Unzip-Scaff-with-repeat.bam
/home/urbe/Tools/OPERA-LG_v2.0.6/bin//filter_repeat.pl FALCON-Unzip-Scaff-with-repeat.bam repeat.dat | /usr/local/bin/samtools view - -h -S -b > FALCON-Unzip-Scaff.bam
rm FALCON-Unzip-Scaff-with-repeat.bam
/home/urbe/Tools/OPERA-LG_v2.0.6/bin/OPERA-LG config > log
Analyzing 1 library: FALCON-Unzip-Scaff.bam
min library mean : 0
minimum contig length is 500
Current library: 1 out of 7
Analyzing file: pairedEdges_no_repeat_i0
Analyzing file: pairedEdges_no_repeat_i1
Analyzing file: pairedEdges_no_repeat_i2
Analyzing file: pairedEdges_no_repeat_i3
Analyzing file: pairedEdges_no_repeat_i4
Analyzing file: pairedEdges_no_repeat_i5
ln -s results/scaffoldSeq.fasta scaffoldSeq.fasta

To resolve this, try downloading blasr version 1.3 above and re-run :)

However, the ancient invasions and migrations to slavery and trade, history is embroidered with events that led to interactions between previously separate populations. Early humans migrated due to many factors such as changing climate and landscape and inadequate food supply. Historical migration of human populations begins with the movement of Homo erectus out of Africa across Eurasia about a million years ago. Homo sapiens appear to have occupied all of Africa about 150,000 years ago, moved out of Africa 70,000 years ago, and had spread across Australia, Asia and Europe by 40,000 years BC. Indo-Aryan migration from the Indus Valley to the plain of the River Ganges in Northern India is presumed to have taken place in the Middle to Late Bronze Age, contemporary to the Late Harappan phase in India (ca. 1700 to 1300 BC). From 180 BC, a series of invasions from Central Asia followed, including those led by the Indo-Greeks, Indo-Scythians, Indo-Parthians and Kushans in the northwestern Indian subcontinent.

Using the recent advance technologies researchers have created a historical atlas of instances of such mixing. They use a sophisticated statistical method for making inferences about human history and infer populations interbredings ( happen over the past 4,000 years) with an ease.

The study published the findings and presented with an interactive map. http://admixturemap.paintmychromosomes.com/

These sort of genomic study have some limilation. It is hard to precisely define sources of mixing when it occurred between genetically similar groups, and scenarios involving multiple waves of mixing over time or between multiple groups can be difficult to tease apart. But it is believed that larger sample sizes will improve resolution. These high resolution will provide a deeper understanding of human history.

Reference:

http://www.sciencemag.org/content/early/2014/01/28/science.1245938

http://www.ncbi.nlm.nih.gov/pubmed/21390129?dopt=Abstract&holding=npg

http://www.csulb.edu/~kmacd/paper-ethnicity.html

Image: Wikipedia

Candidate Requirements: A PhD in computer science, computational biology, Bioinformatics, or a related field, along with sufficient programming skills for prototyping. Experience with next generation sequencing data analysis is required. Candidates with MS degree but with relevant work experience can also be considered. The successful candidate must be motivated and capable of working independently as well as in team environment.

Eligible and interested candidates can email your resumes to rohit at strandls dot com

About Strand Life Sciences: Strand was founded in 2000 by computer science and mathematics professors who recognized the need to automate and integrate life science data analysis through an algorithmic and computational approach. Strand’s solutions for life sciences research are robust and easy to use by the most novice user while powerful and configurable for the bioinformatician. Using its award-winning application development platform, AVADIS®, Strand builds innovative products that enable fast and cutting-edge analysis for basic and clinical research, drug discovery and development.

http://www.avadis-ngs.com/careers