BOL: Related items

Troyanskaya Lab

Tue, 04 Feb 2020 06:40:36 -0600

The goal of our research is to interpret and distill this complexity through accurate analysis and modeling of molecular pathways, particularly those in which malfunctions lead to the manifestation of disease. We are inventing integrative methods for systems-level pathway modeling through integrative analysis of genome-scale datasets. We apply these approaches in studying challenging biological problems, such as how pathways function in diverse cell types and how they change dynamically.

https://function.princeton.edu/

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

UCSC Genome Browser and Blat software !

Jit — Thu, 18 Feb 2016 03:18:57 -0600

This directory contains Genome Browser and Blat application binaries built for standalone
command-line use on various supported Linux and UNIX platforms. To determine which set of binaries
to download, type "uname -a" on the command line to display your machine type. In most cases the
usage statement for the application can be viewed by running the binary with no arguments.

The UCSC Genome Browser and Blat software are free for academic, nonprofit, and personal use. A
license is required for commercial download and installation of these binaries, with the exception
of items built from the following source code directories, which are freely available for all uses:

- kent/src/utils (includes big* tools)
- kent/src/lib
- kent/src/hg/autoSql
- kent/src/hg/autoXml

For information about commercial licensing of the Genome Browser software, see
http://genome.ucsc.edu/license/. The Blat and In-Silico PCR software may be commercially
licensed through Kent Informatics (http://www.kentinformatics.com).

More at http://hgdownload.cse.ucsc.edu/admin/exe/

Address of the bookmark: http://hgdownload.cse.ucsc.edu/admin/exe/

BBMap/BBTools package: Multipurpose tool designed for converting reads or other nucleotide data between different formats.

Jit — Mon, 13 Jun 2016 05:47:21 -0500

Reformatis a member of the BBMap/BBTools package. It is a multipurpose tool designed for converting reads or other nucleotide data between different formats. It supports, and can inter-convert:

fastq
fasta
fasta+qual
sam
scarf (an old Illumina format)
bam (if samtools is installed)
gzip
zip
ascii-33 (sanger)
ascii-64 (old Illumina)
paired files
interleaved files

It is multithreaded and can process data at over 500 megabytes per second, and can accept streams from standard in and write to standard out, allowing it to be easily dropped into the middle of a pipeline for format conversion. Reformat autodetects formats based on file extensions and content, making it very easy to use; and the autodetection can be overridden, allowing flexibility for people who don't like to follow naming conventions, or out-of-spec fastq files with qualities values like -17 or 120.

The program has been gradually expanded, and can now perform various other functions. None of these will break pairing, if the input is paired.

Quality trimming (either or both ends)
Quality filtering
Fixed-length trimming
Generation of histograms (base composition, quality, etc)
Subsampling (to a fraction of input reads, or an exact number of reads or bases)
Changing fasta line-wrapping length
Reverse-complementing (all reads or only read 2)
Adding /1 and /2 suffix to read names
GC-content filtering
Length-filtering
Testing for corrupted interleaved files

Reformat is compatible with any platform that supports Java 1.7 or higher. It also has a bash shellscript for simpler invocation. Typical usage examples:

Reformat fastq into fasta:
reformat.sh in=x.fq out=y.fa

Interleave paired reads:
reformat.sh in1=x1.fq in2=x2.fq out=y.fq

Note - you can actually use a shortcut if paired read files have the same name with a 1 and a 2. This is equivalent to the above command:
reformat.sh in=x#.fq out=y.fq

De-interleave reads:
reformat.sh in=x.fq out1=y1.fq out2=y2.fq

Verify that interleaving appears correct, assuming Illumina namimg conventions:
reformat.sh in=x.fq vint

Convert ASCII-33 to ASCII-64:
reformat.sh in=x.fq out=y.fq qin=33 qout=64

Quality-trim paired reads to Q10 on the left and right ends and discard reads shorter than 50bp after trimming:
reformat.sh in1=x1.fq in2=x2.fq out1=y1.fq out2=y2.fq outsingle=singletons.fq qtrim=rl trimq=10 minlength=50

Subsample 10% of the first 20000 pairs in an interleaved file:
reformat.sh in=x.fq out=y.fq reads=20000 samplerate=0.1 int=t
(in this case "int=t" overrides interleaving autodetection, to ensure reads are treated as pairs)

Pipe in a gzipped sam file and pipe out fasta:
reformat.sh in=stdin.sam.gz out=stdout.fa

Reverse-complement reads:
reformat.sh in=x.fq out=y.fq rcomp

For reformatting a file with very long sequences, Reformat will need more memory; just add the additional flag "-Xmx2g". For example, to change the line-wrapping length on the human genome (which has individual sequences over 200Mbp long) to 70 characters:
reformat.sh -Xmx2g in=HG19.fa.gz out=HG19_wrapped.fa.gz fastawrap=70

For additional functions, please run the shellscript with no arguments, or just read it with a text editor. If you have any questions, please post them in this thread.

For people using a non-bash terminal, you may need to type "bash reformat.sh" instead of just "reformat.sh".
For users of Windows or other platforms that do not support bash shellscripts, replace "reformat.sh" with "java -ea -Xmx200m /path/to/bbmap/current/ jgi.ReformatReads"
for example,
java -ea -Xmx200m C:\bbmap\current\ jgi.ReformatReads in=x.fq out=y.fa

Reformat can be downloaded with BBTools here:
https://sourceforge.net/projects/bbmap/

PLINK2

Jit — Thu, 27 Oct 2016 11:24:45 -0500

This is a comprehensive update to Shaun Purcell's PLINK command-line program, developed by Christopher Chang with support from the NIH-NIDDK's Laboratory of Biological Modeling, the Purcell Lab at Mount Sinai School of Medicine, and others. (What's new?) (Credits.) (Methods paper.)

Address of the bookmark: https://www.cog-genomics.org/plink2/

Statistics Using R with Biological Examples

Neel — Thu, 03 Nov 2016 04:55:41 -0500

This book is a manifestation of my desire to teach researchers in biology a bit more about statistics than an ordinary introductory course covers and to introduce the utilization of R as a tool for analyzing their data. My goal is to reach those with little or no training in higher level statistics so that they can do more of their own data analysis, communicate more with statisticians, and appreciate the great potential statistics has to offer as a tool to answer biological questions.

This is necessary in light of the increasing use of higher level statistics in biomedical research. I hope it accomplishes this mission and encourage its free distribution and use as a course text or supplement.

K Seefeld, May 2007

Information Officer at IIAR, Gujarat

Fri, 04 Nov 2016 05:19:13 -0500

Walk in interview at 10.30 am on Nov 11th, 2016 for the following position at Distributed Information Sub-Centre (DISC) established by Dept. Of Biotechnology, Govt. of India at Indian Institute of Advanced Research, Gandhinagar, Gujarat

Position (scale), qualifications and experience

1. Information Officer (Rs 8000-275-13500): one post
Qualifications and experience: MCA, Post-Graduate in any field of biosciences, bioinformatics with at least two years of experience in working in a bioinformatics setup and good knowledge of linux operating system and computer networking.

General terms and conditions:
1. The above engagements is presently till March 31st, 2017. However on extension of project grant and satisfactory performance of the candidate, your services can be extended beyond March 31st, 2017 based on the terms and conditions of extension of project grant.

2. It is not an appointment with the institute and will not confer any right to the incumbent to any claim implicit or explicit on any position.

3.No TA/DA will be paid for attending the interview. Outstation candidates have to make their own arrangement for their stay.

4.Candidates, appearing in the walk-in interview are requested to bring the hard copy of application addressed to Dr. Anju Pappachan with latest photograph, CV mentioning qualifications, work experience and name of two referees and one page write up as to why you would like to join the project.

Address:
Dr. Anju Pappachan
Indian Institute of Advanced Research (IIAR),
University and Institute of Advanced Research,
The Puri Foundation for Education in India,
Koba Institutional Area,
Gandhinagar- 382 007, Gujarat, India. Contact no. 079-30514152
e-mail- anju@iiar.res.in

RECORD

Bulbul — Fri, 25 Nov 2016 08:23:36 -0600

Background. Next-generation sequencing technologies are now producing multiple times the genome size in total reads from a single experiment. This is enough information to reconstruct at least some of the differences between the individual genome studied in the experiment and the reference genome of the species. However, in most typical protocols, this information is disregarded and the reference genome is used. Results. We provide a new approach that allows researchers to reconstruct genomes very closely related to the reference genome (e.g., mutants of the same species) directly from the reads used in the experiment. Our approach applies de novo assembly software to experimental reads and so-called pseudoreads and uses the resulting contigs to generate a modified reference sequence. In this way, it can very quickly, and at no additional sequencing cost, generate new, modified reference sequence that is closer to the actual sequenced genome and has a full coverage. In this paper, we describe our approach and test its implementation called RECORD. We evaluate RECORD on both simulated and real data. We made our software publicly available on sourceforge. Conclusion. Our tests show that on closely related sequences RECORD outperforms more general assisted-assembly software.

More at https://sourceforge.net/projects/record-genome-assembler/files/

Address of the bookmark: https://www.ncbi.nlm.nih.gov/pubmed/26558255

Statistics and probability

Jit — Tue, 08 Nov 2016 07:34:25 -0600

Topics

Displaying and describing data

Modeling distributions of data

Describing relationships in quantitative data

Designing studies

Probability

Random variables

Sampling distributions

Confidence intervals (one sample)

Significance tests (one sample)

Significance tests and confidence intervals (two samples)

Inference for categorical data (chi-square tests)

Advanced regression (inference and tran

Address of the bookmark: https://www.khanacademy.org/math/statistics-probability

Bioistats Online course

Abhimanyu Singh — Thu, 10 Nov 2016 04:22:51 -0600

One of our primary focuses will be to develop an understanding of the various ways in which we can assign a probability to some chance event. We'll also learn the fundamental properties of probability, investigate how probability behaves, and learn how to calculate the probability of a new chance event.

This book is handy understanding basic concepts.

Address of the bookmark: https://onlinecourses.science.psu.edu/stat414/node/287