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

Post-Doc Positions at the Institute of Evolution, University of Haifa, Haifa, Israel

Thu, 04 Sep 2014 03:59:38 -0500

We are looking for independent, motivated, diligent, laborious, dedicated Bioinformaticians as post-doctorate fellows for a project aimed at revealing the mechanisms of cancer-resistance and anti-cancer activity of the hypoxia-tolerant subterranean, blind mole-rat, Spalax along its underground evolutionary adaptations. Our project has captured the interest of the scientific community and we have ample financial support for the studies. Generous fellowships ($30K to $40K according to qualifications and performance) are available, immediately, for Post-Docs experts in bioinformatics with a background of good understanding biological questions. That is that can independently handle raw output data of RNA-seq / miR seq/ Genomic, analyze it and can interpret intelligently the relevant biological background. Outstanding candidates for PhD experienced in Bioinformatics will also be considered. Familiarity with cancer research is an advantage. Experience of writing manuscripts for publication and a publication record in relevant journals are expected. English skills both oral and written are required. American, Western-European or Israeli education is a significant benefit.

Our present objectives is to identify and isolate the substances secreted by Spalax cells, resolve with which components they interact that are active only on cancer cells, in order to unravel the biological mechanisms and pathways that evolved in Spalax cell machinery and ultimately lead to the death of cancer-cells. The study could attest to be a breakthrough in cancer research, using the long lived, hypoxia- and cancer-tolerant Spalax as a significant biological resource for biomedical research that hopefully could open new horizons in treatment and prevention of cancer in humans.

Contact: The applications should be submitted, together with extended CV and bibliography, summary of past accomplishments, and contact information of 3 referees, to Prof of Research Aaron Avivi (aaron@research.haifa.ac.il) AND Dr. Imad Shams (imadshams@gmail.com). (http://bit.ly/1lywShk) aaron@research.haifa.ac.il

More at http://evolution.haifa.ac.il/index.php/29-people/personal-websites/77-personal-site-avivi

Software engineer/Computational Biologist - Equinome Ltd., Dublin, Ireland

Thu, 04 Sep 2014 19:21:26 -0500

Equinome (www.equinome.com) is the world leader in the research and
development of state-of-the-art novel genomic tools to inform the breeding,
selection and training of Thoroughbred racehorses. Since its launch in 2010,
Equinome has successfully commercialised three performance-related genetic
tests, with a pipeline of further genetic tests in development. We work with
many of the world's leading racehorse trainers and breeders in Europe,
Australasia, USA and South Africa. The company has been featured on CNN,
Bloomberg, RTE, BBC, The Guardian, Discovery Channel and Channel 4, among
others.

The Role

We are looking for a Software Engineer - Computational Biologist with 3+
years' experience in a similar role to design and implement a backend system
to support an online individualised genomics interface. This position is a
great opportunity for an ambitious, self-motivated individual to work in a
demanding, challenging and interesting role.

Position Description:
. Participate in planning, design, and implementation of Equinome back
end systems and technologies.
. Implement interfaces and management tools for back end services.
. Manage, analyse, interpret and visualise large genomics data sets.
. Work closely with scientific team to develop new features and
application enhancements
. Design, develop and manage a genomics research database.

Qualification/Experience:
. Minimum MSc in Computer Science, Genetics, Bioinformatics or in a
related field (A Ph.D qualification would be an advantage).
. Proven 3+ years of experience in similar role.
. Highly proficient in Python, SQL, MySQL.
. Excellent knowledge of mammalian genomics, bioinformatics and
statistical/population genetics.
. Hands-on experience working with large data sets.
. Experience with front-end technologies (HTML/CSS/Javascript) an
advantage.
. Experience in rapid web application development: e.g. Django.
. Knowledge or experience of Unix Scripting and R statistical
programming language would be an advantage.
. Ability to work with minimum supervision to deliver high-quality
code on time.
. Fluency in English and good written and communication skills.
. Meticulous attention to detail.

Applications should be submitted before Friday, 26 September 2014 using the
following link:
http://bit.ly/WgbhxS

Note: Full information and application procedure is available at this link:
http://bit.ly/WgbhxS

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

Project Assistant Position at JMI

Fri, 12 Sep 2014 00:37:44 -0500

Project Assistant Position (@ Rs.10,000/pm Fixed) is available for one year ina research project funded by the Department of Science and Technology entitled, "Folding and stability of naturally truncated photosynthetic pigment,C- phycoerythrin from cyanobacterium Phormidium tenue", at Centre forInterdisciplinary Research in Basic Sciences, lamia Millia Islamia, New Delhi-110025 under' the supervision of Dr. Md. Imtaiyaz Hassan (PrincipalInvestigator).

Eligibility:
M.Sc. in any stream of Life Sciences with minimum 55% marks.

Desirable:
Candidates having experience in Molecular Spectroscopy, Protein Folding and Bioinformatics will be preferred.

Interested candidate may appear in the walk in Interview conducted on September 16, 2014 (Tuesday) 11:00 AM in the Director's Office, Centre for Interdisciplinary Research in Basic Sciences, lamia Millia Islamia, New Delhi-110025.
Candidates are required to bring a set of Xerox copy of their recent CV and qualifying degree (certificate/mark sheet) along with original documents. NoTA/DA will be paid.

For any further information you may e-mail to: mihassan@jmLac.in

Read more at http://jmi.ac.in/upload/advertisement/jobs_cirbs_2014september8.pdf

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

JAMIA HAMDARD Bioinformatics Faculty Jobs 2014

Sat, 20 Sep 2014 21:00:05 -0500

JAMIA HAMDARD

(Deemed University)

Hamdard Nagar, New Delhi – 110 062

R E C R U I T M E N T

(Advertisement No. 5/2014)

Applications on prescribed form are invited for filling up the following teaching positions in the Department of Biotechnology, Faculty of Science in the university. Eligible candidates may apply on or before 30.09.2014.

1. Professor/Associate Professor - One in Pay Band of Rs. 37400-67000+ AGP Rs.10000/9000

2. Assistant Professor - Two in Pay Band of Rs. 15600-39100+ AGP Rs. 6000/-

ASSISTANT PROFESSOR – 02 (including 01 SFS)

Specialization : Bioinformatics

Qualification and Experience :

Ph.D. in Biotechnology or an allied discipline with M.Sc. in Biotechnology/ Biochemistry in the First division or equivalent grade from a recognized University/ Institute.

NET in Life Science or allied discipline in addition to the above qualification.

Experience : At least two years of Post-doctoral teaching and/or research experience in Bioinformatics or relevant field in a UGC recognized Institution of repute or international research institute/ University. Proof of research to be evidenced by publications in SCI-indexed journals of high impact factor as the first or corresponding author.

Note : University may consider exempting candidates from NET, who has been awarded Ph.D. degree from ‘A’ Grade accredited University following the procedure as notified by the UGC in its Regulations of 2009 and adopted by Jamia Hamdard.

For more information: http://www.jamiahamdard.ac.in/PDF/Online%20application%20form%20_Teaching_1.pdf
http://www.jamiahamdard.ac.in/PDF/PBAS.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/

Joao Pedro de Magalhaes Lab

Fri, 26 Sep 2014 19:08:34 -0500

Ageing has a profound impact on human society and modern medicine, yet it remains a major puzzle of biology. The goal of my work is to help understand the genetic, cellular, and molecular mechanisms of ageing. In the long term, I would like my work to help ameliorate age-related diseases and preserve health. No other biomedical field has so much potential to improve human health as research on the basic mechanisms of ageing. Please see our lab website for further details about our work and publications.

Functional and Comparative Genomics

http://jp.senescence.info/
http://www.senescence.info/
http://www.liv.ac.uk/integrative-biology/staff/joao-de-magalhaes/

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