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

National Training on Bioinformatics Computational Tools for Microbial Research Nov 19 to 30, 2013

Fri, 27 Sep 2013 10:49:34 -0500

Agricultural research in modern scientific arena is being represented by proper integration among various research fields of biological, chemical and physical sciences, because this field encompasses many more complexities of biology in nature. In the era of fast accumulating biological data coming out from the research on many crop plants, live stocks and microbes and the impact of changing climate, habitat and other interrelations on these biological entities, bioinformatics has come forward across the globe to solve the problems of analysis, prediction, storage, management, pattern recognition, submission, retrieval and storage of the data to find out a fruitful outcome. This area is becoming increasingly important in the context of systems biology approach where a holistic approach is required to understand the biology and chemistry of the biological entities and their behavior during environmental interactions to resolve the harmful impact of biotic or abiotic causes on crop plants, animals, fishes, livestock sector, beneficial insects as well as microbes. The National Training program on ‘Computational Tools for Microbial Research” is an initiative for the capacity building of NARS scientists/researchers in this most emerging area and fast developing area of i.e. agricultural bioinformatics.

Contact The Director, National Bureau of Agriculturally Important Microorganisms, Kusmaur, Maunath Bhanjan-275101 (U.P.); Phone: 0547-2530080, Fax: 0547-2530358, e mail: nbaimicar@gmail.com; website: www.nbaim.org.in OR

Dr. Dhananjaya P. Singh, Senior Scientist & CCPI, NABG project, NBAIM, Maunath Bhanjan, 275101; Mob.- 09415291703; e mail - dpsfarm@rediffmail.com, nabg.nbaim@gmail.com

More at http://www.nbaim.org.in/Announc.aspx?cd=36

EMBL Postdoc position in Bacterial Gene Gain Loss

Thu, 20 Aug 2015 14:09:21 -0500

A post-doctoral fellowship is available in the research groups of Nick Goldman (EBI) and John Welch (Genetics Department, Cambridge University) under the EMBL-EBI / Cambridge Computational Biomedical Postdoctoral Fellowship scheme.

The project is on bacterial gene gain and loss and emerging pathogenicity, and is described in full here: https://www.ebi.ac.uk/research/postdocs/ebpods/projects/goldman-welch-2015 . The EMBL-EBI / Cambridge Computational Biomedical Postdoctoral (“EBPOD”)

The closing date for applications is 3 September 2015. Nick Goldman EMBL-European Bioinformatics Institute Nick Goldman

More at https://www.ebi.ac.uk/research/postdocs/ebpods/projects/goldman-welch-2015

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/#

Mike Ritchie Lab

Wed, 02 Oct 2013 15:25:45 -0500

Mike Ritchie Lab primary research focus is the detection of susceptibility genes for common diseases such as cancer, diabetes, hypertension, and cardiovascular disease, among others. The approaches will involve the development and application of new statistical methods with a focus on the detection of gene-gene interactions associated with human disease.

Gene expression and protein expression patterns between normal and non-normal tissues is a growing area of research that may lead to the identification of candidate genes for understanding the etiology of common, complex diseases.

Lab homepage @ http://ritchielab.psu.edu/ritchielab/

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

Research Associate @ ICGEB, New Delhi.

Wed, 09 Oct 2013 13:49:20 -0500

Applications are invited for Research Associate position in the DBT Sponsored Bioinformatics Infrastructure Facility at ICGEB, New Delhi.

Essential requirements: Experience of using bioinformatics tools.

Experience of working in Linux. Basic knowledge of computer network administration.

Desirable: Knowledge of Linux installation/administration and proficiency in either of the following:

Shell/PERL/Java/Python/VB/Oracle/MySQL/C/CUDA.

Qualification: PhD. or First class M.Sc degree in Bioinformatics or Biotechnology/life science with specialization in Bioinformatics.

Fellowships: Rs 22,000/- with HRA for PhD qualified, Rs 16000/- with HRA for NET/BET/BINC/GATE qualified and 12000/- with HRA for non NET qualified applicants.

Interested candidates may send their complete biodata along with a write-up of their experience and suitability for the position to Dr. Dinesh Gupta by email only to dinesh@icgeb.res.in within 15 days of publication of this advertisement. Kindly mark the email with subject “Application for BIF-RA-2013”

Closing date for applications: 18 October 2013

Only short listed candidates will be invited for an interview at ICGEB.

No TA/DA will be paid for attending the interview.

Advertisement: http://www.icgeb.org/tl_files/Vacancies/BIF-RA-Advt.pdf

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/

SRF/JRF/RA @ UNIVERSITY OF HYDERABAD

Mon, 14 Oct 2013 07:49:11 -0500

SCHOOL OF CHEMISTRY, UNIVERSITY OF HYDERABAD

Applications on plain paper along with details of CV (relevant photocopies of their
qualifications/experience and reprints of published work to be attached) are invited from qualified candidates for Research Fellowship in CSIR- sponsored research project.

JRF/SRF/RA (one vacancy)

CSIR sponsored “In silico design, identification and in vitro validation of lead molecule inhibitors to Bcr-Abl kinase”

JRF: M.Sc in Chemistry/ Bioinformatics/ Biotechnology with I division and NET or GATE qualified

SRF: M.Sc in chemistry/ Bioinformatics/ Biotechnology with at least two years of post- M.Sc research experience as evidenced from published papers in standard refereed journals in relevant area

RA: PhD in chemistry/ Bioinformatics/ Biotechnology with research experience in
relevant area.

As per CSIR guidelines

Notes:
1) You may visit the University of Hyderabad website www.uohyd.ernet.in to learn more about the University of Hyderabad.
2) Applicants should note that the appointment to be made is purely temporary and there is no right for claiming for any regular appointment in the University.
3) No TA/DA will be paid for attending the interview or at the time of joining the post, if selected.
4) The application should be submitted by post/courier/in-person to the address given below on or before November 1st 2013.

Prof. Lalitha Guruprasad
W-103, Gurbakhsh Singh Building
School of Chemistry
University of Hyderabad
Hyderabad- 500 046
5) Short-listed candidates will be called for interview at a short notice.

Advertisement: http://www.uohyd.ac.in/images/recruitment/chemisry_advt_101013.pdf

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