BOL: Related items

EAGLER: a scaffolding tool for long reads.

Jit — Mon, 04 Jun 2018 05:26:03 -0500

EAGLER is a scaffolding tool for long reads. The scaffolder takes as input a draft genome created by any NGS assembler and a set of long reads. The long reads are used to extend the contigs present in the NGS draft and possibly join overlapping contigs. EAGLER supports both PacBio and Oxford Nanopore reads.

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

https://bib.irb.hr/datoteka/844447.Diplomski_2015_Luka_terbi.pdf

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

JRF in Bioinformatics @ INMAS, DRDO,Delhi

Wed, 01 Oct 2014 07:01:07 -0500

Institute of Nuclear Medicine and Allied Sciences (INMAS), Delhi under the aegis of Defence Research and Development Organisation (DRDO), is engaged in research and developmental work in radiation sciences, Neuro-Computing and Medical Image Processing. INMAS is looking for meritorious young researchers for pursuing research in the frontier areas at INMAS. The Institute invites applications from young and meritorious Indian nationals who are creative, have passion and desire to pursue R&D in frontier areas. INMAS possesses ambience of a research cum academic institute coupled with an advanced R&D infrastructure in a mission mode. It provides the best infrastructure, motivation and personality development prospects for talented students, dreaming of unparalleled success in their professional endeavors. INMAS provides state of the art research facilities for undertaking pioneering research with defence applications.

JRF (Maximum Tenure‐ Five Years: 2yrs as JRF and 3yrs as SRF)
A first class Master’s Degree in Bioinformatics (likely 2 posts)
Around Rs 16,000/ Plus 30% HRA (as per rules of funding agency)

Applications are invited from candidates possessing the above qualifications. The upper age limit is as on the last date for receipt of application. (5 years relaxation to SC/ST candidates, 3 years to OBC candidates, and other entitled categories as per Govt rules). Actual No. of vacancies may vary.

Application form can be download from the website www.drdo.gov.in and E Mailed to inmashrd@gmail.com.
Last date to apply by email is 1700 hrs on 15 Oct 2014
Incomplete applications are liable to be rejected.
Confirmation will be sent to short-listed candidates through email only
Antecedents of selected candidates will be verified.
Written Test will be conducted from 0930-1030 hrs. Latecomers will not be considered.
Candidates will be required to produce certificates/testimonials in original at the time of interview.
It may please be noted that offer of Fellowship does not confer on fellows any right for absorption in DRDO.
Candidates should carry photocopy of Application form sent by email with them.
No TA/DA will be paid for attending interview & on joining.
Last date to apply by email is 1700 hrs on 15 Oct 2014

More at http://drdo.gov.in/drdo/English/jrf29092014.pdf
http://drdo.gov.in/drdo/English/index.jsp?pg=inmas29092014.jsp

Download blasr 1.3 version

Jit — Fri, 15 Jun 2018 03:01:20 -0500

DOWNLOAD LINK: https://github.com/BioInf-Wuerzburg/proovread/raw/master/util/blasr-1.3.1/blasr

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 :)

Internship Program for Bioinformatics / Biotechnology Professionals (No. Of Vacancy: 2)

Wed, 08 Oct 2014 01:10:08 -0500

ArrayGen is offering an Internship Program for Post graduate Bioinformatics / Biotechnology students and professionals. ArrayGen Technologies provide an excellent opportunity to gain research experience and explore if a scientific career is right for you. Currently we offer positions to outstanding students interested in Next Generation Sequencing (NGS) data analysis. Applications are accepted throughout the year. Accepted students will be listed on web with their schedules. Accepted students can attend our future workshops and trainings freely at the specified venue.

How far can bioinformatics go creating organisms used for testing?

Fri, 17 Oct 2014 02:08:16 -0500

"I think you can get very far on a technical level. The problem is that a human body is more complex than just one cell." ... "At some point we still need clinical tests on animals and humans before we use it for real treatment. But we will likely be able to remove 99 % of animal tests in the future." Erik Lindahl, Professor of Theoretical and Computational Biophysics at KTH Royal Institute of Technology is telling us about his work. From the episode "Science for life – mapping the building blocks of the human body". Watch the rest of the talk, and other talks at www.crosstalks.tv Crosstalks is an academic talkshow produced by KTH Royal Institute of Technology and Stockholm University.

Converting a VCF into a FASTA given some reference !

Jit — Fri, 20 Jul 2018 10:03:53 -0500

Samtools/BCFtools (Heng Li) provides a Perl script vcfutils.pl which does this, the function vcf2fq (lines 469-528)

This script has been modified by others to convert InDels as well, e.g. this by David Eccles

./vcf2fq.pl -f <input.fasta> <all-site.vcf> > <output.fastq>

https://github.com/gringer/bioinfscripts/blob/master/vcf2fq.pl

https://github.com/lh3/samtools/blob/master/bcftools/vcfutils.pl

Graduate research assistantships @ University of Nebraska-Lincoln (UNL)

Wed, 22 Oct 2014 10:05:31 -0500

Graduate research assistantships in quantitative genetics are available with Gota Morota in the Department of Animal Science at the University of Nebraska-Lincoln (UNL).

Current projects in the Morota lab include developing kernel-based whole-genome prediction and kernel-based genome-wide association models, polygenic modeling of binary traits, reexamining the results from quantitative genetics analysis in light of functional annotation, and extending kernel methods (such as GBLUP and RKHS) specifically tailored for diverse types of emerging omics data.

In addition, candidates will be expected to leverage opportunities to interact with faculty in animal genetics and biometrics at the UNL in the areas of bioinformatics, breeding, functional genomics, quantitative genetics, and molecular genetics.

Candidates should have a B.S. or M.S. degree in quantitative disciplines with strong background and interest in statistical computing.
The starting date is Fall 2015.
For more information about research in the Morota lab at the UNL, visit: http://www.morotalab.org

A letter of interest in the position, C.V., and contact information for
three references should be emailed to Gota Morota at .
Review of applications will begin immediately, and continue until the
positions are filled. Informal inquiries are also welcome.

Also, please see: http://animalscience.unl.edu/anscprospectivegraduatestudents

LR_Gapcloser: a tiling path-based gap closer that uses long reads to complete genome assembly

Rahul Nayak — Thu, 14 May 2020 15:09:52 -0500

LR_Gapcloser is a gap closing tool using long reads from studied species. The long reads could be downloaed from public read archive database (for instance, NCBI SRA database ) or be your own data. Then they are fragmented and aligned to scaffolds using BWA mem algorithm in BWA package. In the package, we provided a compiled bwa, so the user needn't to install bwa. LR_Gapcloser uses the alignments to find the bridging that cross the gap, and then fills the long read original sequence into the genomic gaps.

Address of the bookmark: https://github.com/CAFS-bioinformatics/LR_Gapcloser

JRF/SRF at University of Calcutta

Fri, 31 Oct 2014 08:53:10 -0500

Applications are invited to appear at a walk-in-interview for one post of Junior Research Fellow in the DBT(DBT Twinning NER) sponsored project entitled “Protein folding kinetics is a selection force on shaping codon usage bias in the high expression genes” in the room of the HOD, Department of Biotechnology and the Coordinator, DR. B. C. Guha Centre for Genetic Engineering and Biotechnology, University College of Science, 35 Ballygunge Circular Road, Kolkata 700019 on the 12th November, 2014 at 3:00 p.m.

Essential qualifications: First class M. Sc. in any branch of life sciences and qualified CSIR-UGC NET/GATE Examination.

Desirable qualifications: Practical experience in biochemical and biophysical studies of proteins

Emoluments: as per DBT norms

The project is tenable for two years, initially for one year.

Age: Below 28 years (relaxable in the case of SC/ST/OBC/women candidates)

Candidates are requested to bring two sets of complete applications on plain paper furnishing bio-data and copies of attested certificates along with originals (for verification) on the date of interview.

No TA/DA is admissible for candidates appearing at the interview.

Dr. Rajat Banerjee
Assistant Professor
Department of Biotechnology and
Dr. B. C. Guha Centre for Genetic Engineering and Biotechnology
University College of Science
35, Ballygunge Circular Road
Kolkata 700019

Advertisement: www.caluniv.ac.in/news/jrf_biotech_2.pdf

Long read assembly workshop !

Rahul Nayak — Thu, 04 Oct 2018 17:23:18 -0500

This is a tutorial for a workshop on long-read (PacBio) genome assembly.

It demonstrates how to use long PacBio sequencing reads to assemble a bacterial genome, and includes additional steps for circularising, trimming, finding plasmids, and correcting the assembly with short-read Illumina data.

Please comment if you know any other long read addembly tutorial.

Address of the bookmark: http://sepsis-omics.github.io/tutorials/modules/cmdline_assembly_v2/