BOL: Related items

Mugsy: multiple whole genome alignment tool

Jit — Fri, 08 Dec 2017 17:41:14 -0600

Mugsy is a multiple whole genome aligner. Mugsy uses Nucmer for pairwise alignment, a custom graph based segmentation procedure for identifying collinear regions, and the segment-based progressive multiple alignment strategy from Seqan::TCoffee. Mugsy accepts draft genomes in the form of multi-FASTA files and does not require a reference genome.

To cite Mugsy, use:

Angiuoli SV and Salzberg SL. Mugsy: Fast multiple alignment of closely related whole genomes.Bioinformatics 2011 27(3):334-4

Address of the bookmark: http://mugsy.sourceforge.net/

MATHomics Lab

Tue, 19 Nov 2013 18:17:32 -0600

Mathomics is a collaborative research group of the Center for Mathematical Modeling and the Center for Genome Regulation at University of Chile, created to play a central role in the development of biotechnological projects, providing state of the art bioinformatics and mathematical modeling tools, allowing to face these problems from the point of view of Systems Biology.

Lab page @ http://www.mathomics.cl/

Magic-BLAST: a tool for mapping large next-generation RNA or DNA sequencing runs against a whole genome or transcriptome.

Jit — Tue, 26 Dec 2017 22:23:39 -0600

Magic-BLAST is a tool for mapping large next-generation RNA or DNA sequencing runs against a whole genome or transcriptome. Each alignment optimizes a composite score, taking into account simultaneously the two reads of a pair, and in case of RNA-seq, locating the candidate introns and adding up the score of all exons. This is very different from other versions of BLAST, where each exon is scored as a separate hit and read-pairing is ignored.

Magic-BLAST incorporates within the NCBI BLAST code framework ideas developed in the NCBI Magic pipeline, in particular hit extensions by local walk and jump (http://www.ncbi.nlm.nih.gov/pubmed/26109056), and recursive clipping of mismatches near the edges of the reads, which avoids accumulating artefactual mismatches near splice sites and is needed to distinguish short indels from substitutions near the edges.

Address of the bookmark: https://ncbi.github.io/magicblast/

Research Assistant @ University of Hyderabad

Mon, 25 Nov 2013 10:21:26 -0600

University of Hyderabad
Repository for Tomato Genomic Resources
Department of Plant Sciences
Bioinformatics Position in Tomato Functional Genomics

At the Repository for Tomato Genomics Resources, we are working on Tomato Functional Genomics, using TILLING, Insertional Mutagenesis, proteomics, metabolomics approaches to study fruit ripening in tomato. The current aims of the group include using reverse and forward genetics strategies to isolate tomato mutants delayed in ripening, having high lycopene and folate content in tomato fruits and analysis of light and hormonal signal transduction pathways. For recent publications of the group see (Plant Physiol 161: 2085–2101, Plant Physiol 156: 1424-1438; Molecular Plant 3: 854-869; Plant Methods 6: 3; Plant Methods 5:18; Plant Signaling and Behavior 5:11.).

Currently we have one position available in the projects awarded to Prof. R.P. Sharma funded by Dept of Biotechnology. The qualification for this Position is as follows:

Research Assistant: Applicants should have experience in networking using R language and should be able to develop networks using the transcriptome, proteome and metabolite data sets. M.Tech. in Bioinformatics is required. The selected candidate would be paid Rs. 13,000/-pm- consolidated.

Candidates interested in above positions should send a one page statement clearly explaining how their skills are relevant to the position. The candidates should also enclose detailed CV and the name/email id for three referees. The candidates can send their application by email at rameshwar.sharma@uohyd.ac.in and y.sreelakshmi@uohyd.ac.in on or before December 10th, 2013. The position is purely temporary in nature. Shortlisted candidates would be called for interview. No TA/DA would be provided for attending the interview. We also have openings for CSIR-NET JRF candidates for pursuing PhD in above research areas.

Interested candidates with CSIR-NET JRF can send their CV to the above email
addresses.

http://www.uohyd.ac.in/images/recruitment/tomanet_positions_221113.pdf

Carefully opt for human reference genome

biogeek — Tue, 18 Feb 2020 07:43:32 -0600

Heng Li posted several issues with the human reference genomes given in these resources and suggests the following compressed FASTA file to be used as hg38/GRCh38 human reference genome.

if you map reads to GRCh38 or hg38, use the following:

ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCA/000/001/405/GCA_000001405.15_GRCh38/seqs_for_alignment_pipelines.ucsc_ids/GCA_000001405.15_GRCh38_no_alt_analysis_set.fna.gz

There are several other versions of GRCh37/GRCh38. What’s wrong with them? Here are a collection of potential issues:

More at http://lh3.github.io/2017/11/13/which-human-reference-genome-to-use

Address of the bookmark: http://lh3.github.io/2017/11/13/which-human-reference-genome-to-use

BOKU Lab

Wed, 08 Jan 2014 19:33:12 -0600

We are interested in the study of complex systems in living organisms. Novel views augmenting the classical gene by gene approaches are required to overcome the engineered redundancies and combinatorial effects prevalent in higher eukaryotes. We therefore combine work to establish improved quantitative experimental assays, such as microarrays or differential in-gel electrophoresis, and development of modern computational methods, such as hierarchical probabilistic models or integration of heterogeneous data sources, focussed by biological studies in our laboratory and collaborations.

Highlights of our research include:

Optimization of microarray design, probe signal interpretation
Advanced models and tools for expression profiling
State-of-the-art applications and integrated analyses

Lab page @ http://bioinf.boku.ac.at/

Metassembler: merging and optimizing de novo genome assemblies

Rahul Nayak — Tue, 08 May 2018 04:52:33 -0500

Metassembler combines multiple whole genome de novo assemblies into a combined consensus assembly using the best segments of the individual assemblies.

Genome assembly projects typically run multiple algorithms in an attempt to find the single best assembly, although those assemblies often have complementary, if untapped, strengths and weaknesses. We present our metassembler algorithm that merges multiple assemblies of a genome into a single superior sequence.

Address of the bookmark: https://sourceforge.net/projects/metassembler/?source=directory

Postdoctoral Position (Bioinformatics/Computational Biology)

Thu, 12 Dec 2013 17:58:00 -0600

University College Cork
LAPTI
Cork-Co Cork-Ireland

Postdoctoral position is available for three years to work on development of Bioinformatics resources for the analysis and visualization of ribosome profiling data. Ribosome profiling (ribo-seq) is a technology that allows mapping positions of the ribosomes on the whole transcriptome level with a nucleotide precision. The technology allows obtaining high resolution digital snapshots of gene expression in cells. The position is available starting on the 1st of October, 2013.

Candidate is expected to have Ph.D. in Bioinformatics or Computational Biology. Candidates with the degree in non-Biological disciplines such as Computer Science, Statistics, Applied Mathematics, Physics or Electrical Engineering will also be considered.

The position is available at LAPTI (http://lapti.ucc.ie) that is located in the Western Gate Building (http://www.stwarchitects.com/project-information.php?c=1&p=09993) at University College Cork. Western Gate Building Research Complex hosts several UCC departments and provides ideal environment for interdisciplinary research. Cork (sometimes referenced as “Venice of Ireland”) is the second most populous city in the Republic. It has friendly cosmopolitan atmosphere and vibrant culture. A number of American industrial giants such as Apple , EMC and Pfizer have chosen Cork as a home for their European headquarters.

The details of the application process are given at http://lapti.ucc.ie/jobs.html. To ensure prompt processing of your application use the subject line: ‘Postdoc computational’. All applications received prior to August the 1st are guaranteed equal consideration. However, applications at the later dates will also be considered until the position is filled.

For more info visit http://lapti.ucc.ie

Critical to discoveries in bioinformatics

Mon, 16 Dec 2013 17:13:24 -0600

EMBL-EBI distributes datasets worldwide using the Janet network. This biological data enables the discovery of new drugs, new diagnostics and increasingly new agro-chemicals. Their work, which includes the 1000-genome project, has generated petabytes of data and this growth is showing no signs of abating. On-demand bandwidth over Janet will therefore be critical to their ongoing work.

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