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

Java utilities for Next Generation Sequencing by Pierre Lindenbaum

Robert M Willioms — Mon, 15 Sep 2014 17:24:49 -0500

Java utilities for Bioinformatics

https://github.com/lindenb/jvarkit

Address of the bookmark: https://github.com/lindenb/jvarkit

KAU Thrissur Biotech/Bioinformatics RA/SRF/JRF/Trainee/Studentships

Fri, 26 Sep 2014 20:07:28 -0500

Applications are invited from eligible candidates for the following posts at Bioinformatics Centre (DIC), IT- BT Complex, College of Horticulture, Kerala Agricultural University, Vellanikkara, Thrissur.

1. Research Associate
Emoluments*: 14880/- + HRA
Qualification needed: Ph.D/M.Sc in Bioinformatics or Ph.D/M.Sc in Agriculture or Biotechnology with advanced Diploma in Bioinformatics
Desirable: 2 year experience in Bioinformatics.

2 Senior Research Fellow
Emoluments*: 10230/-
Qualification needed: M.Sc/ M.tech in Bioinformatics or M.Sc in Agriculture/ Biotechnology with Diploma in Bioinformatics.
Desirable: One year experience in Bioinformatics

3 Junior Research Fellow
Emoluments*: 9300/-
Qualification needed: M.Sc/ M.tech in Bioinformatics or M.Sc in Agriculture/Biotechnology/Plant Sciences with Diploma in Bioinformatics.

4 .Trainee/Studentship Bioinformatics
Emoluments*: 5000/-
Qualification needed: M.Sc in Bioinformatics with good knowledge of Bioinformatics softwares and tools.

5 Trainee/ Studentship Biotechnology
Emoluments*: 5000/-
Qualification needed: M.Sc in Biotechnology, with working knowledge in tissue culture, molecular markers and cloning of genes.

Candidates with the required qualifications and experience may give an application in the prescribed format with attested copies of certificates to prove eligibility on or before 30th November 2014. The applications are to be addressed to The Associate Dean, College of Horticulture and send to "Professor & Coordinator, Bioinformatics Centre (DIC), IT-BT Complex, Kerala Agricultural University, Vellanikkara, Thrissur, Kerala 680 656”. The envelope may be superscribed “Application for the post at Bioinformatics Centre”.

*Emoluments are likely to be revised in 2014-2015

More at http://www.kaubic.in/downloads/Notification_bic.pdf
http://www.kaubic.in/downloads/Application%20form.pdf

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

FACULTY POSITION at NATIONAL INSTITUTE OF TECHNOLOGY, ROURKELA

Sun, 05 Oct 2014 15:45:13 -0500

NIT, Rourkela, an institute of national importance under Ministry of HRD, Govt. of India invites applications from Indian nationals possessing excellent academic background along with commitment to quality teaching and research for faculty positions at the level of Professor, Associate Professor and Assistant Professor in most branches of Engineering, Science, Management and Humanities as per following table:-

1 Asst. Professor (on Pre-Ph.D. Contract)

2 Asst. Professor (on Contract)

3 Asst. Professor

4 Associate Professor

5 Professor

Life Sciences:-

i)Biochemistry and Molecular Biology; ii)Cell and Developmental Biology; iii)Immunology and Molecular Medicine; (iv) Microbiology and Ecology (v)Bioinformatics Group; vi)Biophysical Sciences

HOW TO APPLY:-

a. Candidates willing to apply for one or more posts are requested to apply online at “http://www.nitrkl.ac.in/ JOBS & TENDERS /Faculty Position” .
b. Persons employed in Government and Semi-Government organizations may apply directly against the standing advertisement. For this the application should be completed online. The printout of the application generated online should be submitted through employer if shortlisted for interview.
c. The online application can be filled in multiple sessions.
d. Candidates are required to check the Institute website from time to time for latest information, application status call for interview, change of dates and final results.
e. Applications shall be received online only. Please do not send application or CV against this advertisement by email or letter mail.

More Info: http://nitrkl.ac.in/Jobs_Tenders/1FacultyPosition/Default.aspx

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

Software developed in pevsner lab

Robert M Willioms — Mon, 06 Oct 2014 12:41:26 -0500

DRAGON: Database Referencing of Array Genes Online

SNOMAD: Standardization and Normalization of Microarray Data

SNPduo: SNP Analysis Between Two Individuals

SNPtrio: Analyzing and Visualizing and Inheritance Patterns in Trios

SNPscan: Data Analysis and Visualization of SNP Data

pediSNP: Analyze SNP Data From a Pedigree of Two Generations

kcoeff: Calculate Cotterman Coefficients of SNP Genotype Data

triPOD: Detects chromosomal abnormalities in parent-child trio-based microarray data

Address of the bookmark: http://pevsnerlab.kennedykrieger.org/php/?q=software

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/

Research Fellows at CSIR - Institute of Himalayan Bioresource Technology, Palampur, Himachal Pradesh

Tue, 19 May 2015 07:17:32 -0500

CSIR - Institute of Himalayan Bioresource Technology 2 vacancies of Project Fellow

Name of the Post: Project Fellow

No. of the Post: 02 Two

Salary: Rs. 12000/- per month or Rs. 14000/- per month

Age Limit: Max. 28 years as on 10.06.2015 and relaxation as per rules

Required Job Profile:

Candidate must possess first class B.Tech. in bioinformatics or computational biology OR M.Sc. in bioinformatics or computational biology with fifty five percent marks OR M.Tech. in bioinformatics or computational biology with fifty five percent marks.

How to apply:

Eligible and interested candidates should need to appear for walk-in interview on 10.06.2015 at 9:00 am at the above mentioned address.

Refer to http://www.ihbt.res.in/recruit/AdvtNo7_2015.pdf

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