BOL: Related items

Queensland Centre for Medical Genomics, Grimmond Lab

Sun, 14 Jul 2013 11:58:34 -0500

Queensland Centre for Medical Genomics

Research Area:
pancreatic cancer; ovarian cancer; prostate cancer; bowel cancer; brain cancer; endometrial cancer; breast cancer; personalised medicine; high-throughput genomics

Link @ http://www.imb.uq.edu.au/sean-grimmond

Complex Systems: From Physics to Biology October 15-16 2013 at JNU Convention Center

Mon, 23 Sep 2013 10:17:17 -0500

The symposium intents to focus on complex systems arising in a variety of settings in physics and biology. In particular, applications of the concepts of physics to biological sciences will be the major theme of this meeting.

Selected Topics:

Cluster Dynamics
Non-equilibrium Statistical Mechanics
Forced Systems
Hamiltonian Dynamics
Synchronization & Control
Genomics & Systems Biology
Computational Neuroscience
Econophysics

More @ http://www.jnu.ac.in/Conference/SCS2013/

Rolland-Lagan lab

Sun, 14 Jul 2013 12:57:57 -0500

The Rolland-Lagan lab at the University of Ottawa is specializing in computational and developmental biology. We use a combination of experimental work, microscopy, image analysis and computer simulations to explore developmental mechanisms in two and three dimensions.

Research Area

Developmental biology, Computational biology, Simulation modeling, Image data analysis

Link @ http://mysite.science.uottawa.ca/arolland/index.html

MashMap: a fast and approximate software for mapping long reads (PacBio/ONT) or assembly to reference genome(s)

Jit — Tue, 12 Dec 2017 17:23:31 -0600

MashMap is a fast and approximate software for mapping long reads (PacBio/ONT) or assembly to reference genome(s). It maps a query sequence against a reference region if and only if its estimated alignment identity is above a specified threshold. It does not compute the alignments explicitly, but rather estimates a k-mer based Jaccard similarity using a combination of Winnowing and MinHash. This is then converted to an estimate of sequence identity using the Mash distance. An appropriate k-mer sampling rate is automatically determined given minimum local alignment length and identity thresholds. The efficiency of the algorithm improves as both of these thresholds are increased.

Address of the bookmark: https://github.com/marbl/MashMap

Frequent Paired-end reads (PE 2x100) mapping command lines

Jit — Tue, 15 May 2018 08:59:29 -0500

bowtie2 -x hs37m -X 650 -q -1 r1.fq -2 r2.fq -S r12.bowtie2.sam

bwa aln hs37m.fa r1.fq > r1.sai && bwa aln hs37m.fa r2.fq > r2.sai \
&& bwa sampe hs37m r1.sai r2.sai r1.fq r2.fq > r12.bwa.sam

bwa bwasw ../index/bwa/hs37m.fa r12.fq > r12.bwasw.sam

gsnap -A sam -d hs37m r1.fq r2.fq > r12.gsnap.sam

novoalign -r Random -o SAM -f r1.fq r2.fq -i 500 50 -d hs37m-k14s3.novo > r12.novo.sam

smalt map -f samsoft -i 650 -o r12.smalt-k20s13.sam hs37m-k20s13 r1.fq r2.fq

stampy.py -g hs37m -h hs37m -o r12.stampy.sam -M r1.fq,r2.fq

soap -D hs37m.fa.index -a r1.fq -b r2.fq -l 32 -g 3 -u dummy -2 dummy -o r12.soap

COSINE: non-seeding method for mapping long noisy sequences

Jit — Fri, 26 Oct 2018 00:41:59 -0500

Third generation sequencing (TGS) are highly promising technologies but the long and noisy reads from TGS are difficult to align using existing algorithms. Here, we present COSINE, a conceptually new method designed specifically for aligning long reads contaminated by a high level of errors.

Address of the bookmark: https://github.com/SUwonglab/COSINE

NCBI Magic-BLAST

Jit — Tue, 14 Aug 2018 18:11:11 -0500

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/

Common steps for reads mapping !

BioStar — Thu, 09 Mar 2023 02:48:02 -0600

Mapping reads to a reference genome is an essential step in many types of genomic analysis, such as variant calling and gene expression analysis. Here are some general steps to follow for mapping reads to a genome:

Choose a read mapper: There are many read mappers available, such as BWA, Bowtie, and HISAT2. Choose a mapper that is appropriate for your type of data and research question.
Index the reference genome: Before mapping reads, the reference genome needs to be indexed. This involves creating an index of the genome sequence that allows the mapper to quickly find matches to the reads. Most mappers have their own indexing tools.
Prepare the read data: The reads should be in a format that is compatible with the mapper. Most mappers accept FASTQ or BAM files. Depending on the quality of the data, it may need to be filtered or trimmed before mapping.
Run the mapper: The mapper is run with the command-line interface or using a graphical user interface. The specific command depends on the mapper being used, but typically involves specifying the input data, reference genome, and output file format.
Evaluate the mapping results: After the mapping is complete, the results should be evaluated. This includes assessing the quality of the mapping, such as the mapping rate, the number of mapped reads, and the mapping quality score.
Post-processing: Depending on the analysis being performed, post-processing of the mapped reads may be necessary. This can include filtering reads based on quality, removing duplicate reads, and calling variants.

Overall, mapping reads to a reference genome is a complex process that requires careful consideration of the type of data, the research question, and the specific mapper being used.

PHD SCHOLARSHIP DENMARK

Fri, 13 Jun 2014 13:44:07 -0500

ne PhD position is available at the Bioinformatics Center, Department of Biology, University of Copenhagen, Denmark. The PhD position concerns protein structure prediction, and will be in the Structural Bioinformatics group of Associate professor Thomas Hamelryck. The group is an integrated part of the Bioinformatics Center, which is headed by Professor Anders Krogh, employs around sixty scientists (including PhD students) and focuses on non-coding RNA, eukaryotic gene regulation and protein structure prediction. The center provides a modern, pleasant, international working environment with excellent modern facilities, in the heart of Copenhagen.
The project will be supervised by Associate Professor Thomas Hamelryck.

The protein folding problem is of enormous practical, theoretical and medical importance - and in addition forms a fascinating intellectual challenge. The aim of this project is to develop and implement a probabilistic method to infer the structure of proteins, building on various probabilistic models of protein structure developed by the Hamelryck group. The method will also take the dynamic nature of proteins into account, and involves a close collaboration with the statistics department at the university within the interdisciplinary project "Dynamical Systems: Mathematical Modeling and Statistical Methods for the Social, Health, and Natural Sciences" (http://dsin.ku.dk/).

Qualifications
Knowledge of programming (C++) and statistics or machine learning. Knowledge of biology, physics or biophysics is a plus, but not a requirement.

The deadline for applications is June 15, 2014

More at : https://job.jobnet.dk/CV/FindJob/details.aspx/3695051%20

Gblocks: eliminates poorly aligned positions and divergent regions of a DNA or protein alignment

Poonam Mahapatra — Sat, 02 Jun 2018 07:36:05 -0500

Gblocks eliminates poorly aligned positions and divergent regions of a DNA or protein alignment so that it becomes more suitable for phylogenetic analysis. This server implements the most important features of the Gblocks program to make its use as simple as possible without loosing the functionality that it is necessary in most of the cases. Other options can be changed in the stand-alone program. You can see here an example output file showing the blocks selected from a protein alignment. Further information can be found in the online documentation.

Address of the bookmark: http://molevol.cmima.csic.es/castresana/Gblocks_server.html