BOL: Related items

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

Bioinformatics Algorithm Demonstrations and Tutorials

Jitendra Narayan — Thu, 29 Aug 2013 09:23:51 -0500

Abstract

This project presents demonstrations of selected computer science algorithms important in bioinformatics, implemented in the spreadsheet program Microsoft Excel. Spreadsheets provide an interesting platform for demonstration of algorithms, since various steps of the calculations can be exposed in a manner that is easily comprehensible to users with little programming experience. The algorithms demonstrated include two approaches to approximate string matching (dynamic programming and Shift-AND numeric approximate matching), Hierarchical Clustering (used in phylogenetic studies and microarray analysis of gene expression), a Naive Bayes Classifier for simulated microarray gene expression data, and a simple Neural Network. These demonstrations are designed to serve as instructional aids in bioinformatics courses.

Tutorial @ http://www.cybertory.org/downloads/bae/BioinformaticsAlgorithmsInExcel.zip

One of the best resource for online bioinformatics learning is https://stepic.org/Bioinformatics-Algorithms-2 Enjoy the online learning.

Reference : cybertory

" Please add your favourite bioinformatics algorithms and tutorial links below in the comment section, for the benefit of bioinformatics and computational biology community ".

Address of the bookmark: http://www.cybertory.org/downloads/bae/BioinformaticsAlgorithmsExcelDoc.pdf

4273π: Bioinformatics education on low cost ARM hardware

Jitendra Narayan — Mon, 02 Sep 2013 07:02:43 -0500

Are you teaching bioinformatics at universities and found it complicated by typical computer classroom settings. As well as running software locally and online, students should gain experience of systems administration. Hmm don't worry there is one new OS for the rescue. 4273π, an operating system image for Raspberry Pi based on Raspbian Linux. It provides an attractive, general-purpose computing environment, within which the course 4273π Bioinformatics for Biologists is embedded.

Though far slower than current desktop and laptop computers, the Raspberry Pi is notably faster than the Cray 1 supercomputer, a marvel of computer speed in its day. The Raspberry Pi approach includes all the benefits of the laptop approach, above, but at lower cost. In addition, the Raspberry Pi is a new and exciting computer system, which in itself can add interest to the course.

As the Raspbian operating system, Raspberry Pi firmware and hardware and 4273π Bioinformatics for Biologists teaching material develop, further releases of 4273π will be made available. It is anticipated that there will be a minimum of two releases per year during the next four years.

4273π is a means to teach bioinformatics, including systems administration tasks, to undergraduates at low cost.

Descriptive paper @ http://www.biomedcentral.com/1471-2105/14/243

Image source: BMC Bioinformatics

P_RNA_scaffolder: a fast and accurate genome scaffolder using paired-end RNA-sequencing reads

BioStar — Fri, 07 Sep 2018 05:19:06 -0500

P_RNA_scaffolder is a novel scaffolding tool using Pair-end RNA-seq to scaffold genome fragments. The method is suitable for most genomes. The program could utilize Illumina Paired-end RNA-sequencing reads from target speciesies. Our method provides another practical alternative to existing mate-pair_based approaches or other Protein-based approaches (for instance, PEP_scaffolder ) for scaffolding genome sequences. The most important feature of this method is to improve the completeness of gene regions and long-coding gene regions (for instance, circRNA).

Address of the bookmark: http://www.fishbrowser.org/software/P_RNA_scaffolder/#

Huber Lab

Mon, 09 Sep 2013 21:57:03 -0500

The Huber group develops computational and statistical methods to design and analyse novel experimental approaches in genetics and cell biology.

Future projects and goals

Large-scale systematic maps of gene-gene and gene-environment interactions by automated phenotyping, using image analysis, machine learning, sparse model building and causal inference.
DNA-, RNA- and ChIP-Seq and their applications to gene expression regulation: statistical and computational foundations.
Cancer genomics, genomes as biomarkers, cancer phylogeny.
Image analysis for systems biology: measuring the dynamics of cell cycle and of cell migration of individual cells under normal conditions and many different perturbations (RNAi, drugs).

More @ http://www.embl.de/research/units/genome_biology/huber/index.html

Synima: a Synteny imaging tool for annotated genome assemblies

Abhimanyu Singh — Tue, 30 Oct 2018 10:49:13 -0500

Synima written in Perl, which uses the graphical features of R. Synima takes orthologues computed from reciprocal best BLAST hits or OrthoMCL, and DAGchainer, and outputs an overview of genome-wide synteny in PDF. Each of these programs are included with the Synima package, and a pipeline for their use. Synima has a range of graphical parameters including size, colours, order, and labels, which are specified in a config file generated by the first run of Synima – and can be subsequently edited. Synima runs quickly on a command line to generate informative and publication quality figures. Synima is open source and freely available from https://github.com/rhysf/Synima under the MIT License.

Address of the bookmark: https://github.com/rhysf/Synima

Asst. PROF IN BIOINFORMATICS at JAIPUR NATIONAL UNIVERSITY

Thu, 12 Sep 2013 07:18:02 -0500

JAIPUR NATIONAL UNIVERSITY, SCHOOL OF LIFE SCIENCES (SIILAS CAMPUS) URGENTLY REQUIRES

Asst. PROF IN BIOINFORMATICS.

QUALIFICATION: AS PER UGC

DESIRABLE: 1 YEAR EXPERIENCE IN ACADEMICS

CONTACT immediately

Prof D.S.Bhatia
Director
9351288070

Last date within 7 days of the publication.

Find more @ http://jnujaipur.ac.in/downloads/AdvtDec2012.jpg

ANItools web: a web tool for fast genome comparison within multiple bacterial strains

Jit — Wed, 14 Nov 2018 04:34:23 -0600

ANItools is a software package written by PERL scripts that can be run in a Linux/Unix system. If you want to compare bacterial genomes and calculate their average nucleotide identity (ANI), you could download and run this program directly. Or you could send us the genome sequence by email. Then we will do the analysis work for you.

https://academic.oup.com/database/article/doi/10.1093/database/baw084/2630454

Address of the bookmark: http://ani.mypathogen.cn/

Research Associate@ School of Computational and Integrative Sciences under Jawaharlal Nehru University

Fri, 22 Nov 2013 19:06:44 -0600

School of Computational and Integrative Sciences under Jawaharlal Nehru University, New Delhi invited applications for filling up 4 posts of Research Associates (RA) and Junior Research Fellow (JRF) (2 posts each) purely on temporary basis, liable to be terminated at any time without prior notice or ceased/withdrawn by the funding agency. The vacancies are for a Department of Biotechnology, Government of India funded project entitled "Computational Core for Plant Metabolomics" (Project ID: 632) being administered by Prof Indira Ghosh. Interested candidates should send their applications till 13 December 2013.
Important Dates
Last Date for receipt of applications: 13 December 2013
Vacancy Details
Total Vacancies: 4 posts
Type of recruitment: Temporary
Sl. No.: 01
Name of the Post: Research Associate
No of Posts: 1 post
Remuneration: Rs. 23000 + 30%
Qualifications: PhD in Bioinformatics / computational biology / Biophysics / Physical Chemistry / Computer Science. Computational experience, proven by paper published, is a necessary qualification.
Sl. No.: 02
Name of the Post: Research Associate
No of Posts: 1 post
Remuneration: Rs. 23000 + 30%
Qualifications: PhD in Computational Biology / Bioinformatics & related subjects. Computational experience, proven by paper published, is a necessary qualification.
Sl. No.: 03
Name of the Post: Junior Research Fellow
No of Posts: 1 post
Remuneration: Rs. 12000 + 30%
Qualifications: M. Sc. / B. Tech. preferably in Computational Biology /Bioinformatics and related fields with experience in Website designing & maintenance of Database.
Sl. No.: 04
Name of the Post: Junior Research Fellow
No of Posts: 1 post
Remuneration: Rs. 12000 + 30%
Qualifications: M. Sc. / MCA / B. Tech. preferably in Computational Biology / Computer science with experience in Programming in Java / Python, C++ etc & designing of Database.
Selection Procedure: Selection will be done on the basis of candidates’ performance in the Interview.
Candidates short-listed / selected for Interview will be informed through email only.
How to Apply: Interested eligible candidates should send their applications, in the prescribed format, along with their current CV by post to “Prof Indira Ghosh, Project Investigator, Hall#6, School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi-110 067” so as to reach the concerned authority by 13 December 2013.
Name of the post applied for’ must be superscripted on the envelope containing the application.
NOTE: For the post of Research Associates, only those candidates who have submitted thesis are eligible to apply. However, salary will be provided as per DBT / DST guidelines (i.e. candidates who have qualified NET /BET / BINC will have higher pay scale).
Candidates interested to register for PhD may not apply for JRF.
More @ http://www.jnu.ac.in/Career/currentjobs.htm

Purge Haplotigs: Pipeline to help with curating heterozygous diploid genome assemblies

Rahul Nayak — Mon, 17 Dec 2018 03:17:20 -0600

Some parts of a genome may have a very high degree of heterozygosity. This causes contigs for both haplotypes of that part of the genome to be assembled as separate primary contigs, rather than as a contig and an associated haplotig. This can be an issue for downstream analysis whether you're working on the haploid or phased-diploid assembly.

Identify pairs of contigs that are syntenic and move one of them to the haplotig 'pool'. The pipeline uses mapped read coverage and Minimap2 alignments to determine which contigs to keep for the haploid assembly. Dotplots are optionally produced for all flagged contig matches, juxtaposed with read-coverage, to help the user determine the proper assignment of any remaining ambiguous contigs. The pipeline will run on either a haploid assembly (i.e. Canu, FALCON or FALCON-Unzip primary contigs) or on a phased-diploid assembly (i.e. FALCON-Unzip primary contigs + haplotigs). Here are two examples of how Purge Haplotigs can improve a haploid and diploid assembly.

Address of the bookmark: https://bitbucket.org/mroachawri/purge_haplotigs