BOL: Related items

Bioinformatician’s Pocket Reference !!

RAJESH DETROJA — Sun, 08 Jun 2014 09:56:58 -0500

It is amusing how brain of bioinformaticians work! Learning a new programming language for days feels so much of fun that making 5 minute discussion with neighbours (unless under special circumstances!) in our own mother-tongue. Today every bioinformatician keeps more than few languages and core IT toolkits on their plate. It has become mandatory to be able to mould different code snippets to build our own custom workflows, and thus keeping syntax at our fingertips has become essential.Although Google is best way to get syntax problem solved, it is not a bad idea to keep reference sheets is our smartphones or stick out some printed sheets on the back of your door, in the old fashion way!!

Address of the bookmark: http://infoplatter.wordpress.com/2014/04/06/bioinformaticians-pocket-reference/

Perl Special Vars Quick Reference

Abhimanyu Singh — Tue, 07 Feb 2017 05:08:47 -0600

`$_`	The default or implicit variable.
`@_`	Subroutine parameters.
`$a` `$b`	sort comparison routine variables.
`@ARGV`	The command-line args.
Regular Expressions
`$`	Regexp parenthetical capture holders.
`$&`	Last successful match (degrades performance).
`${^MATCH}`	Similar to `$&` without performance penalty. Requires /p modifier.
$`	Prematch for last successful match string (degrades performance).
`${^PREMATCH}`	Similar to $` without performance penalty. Requires `/p` modifier.
`$'`	Postmatch for last successful match string (degrades performance).
`${^POSTMATCH}`	Similar to `$'` without performance penalty. Requires `/p` modifier.
`$+`	Last paren match.
`$^N`	Last closed paren match (last submatch).
`@+`	Offsets of ends of successful submatches in scope.
`@-`	Offsets of starts of successful submatches in scope.
`%+`	Like `@+`, but for named submatches.
`%-`	Like `@-`, but for named submatches.
`$^R`	Last regexp (?{code}) result.
`${^RE_DEBUG_FLAGS}`	Current value of regexp debugging flags. See `use re 'debug';`
`${^RE_TRIE_MAXBUF}`	Control memory allocations for RE optimizations for large alternations.
Encoding
`${^ENCODING}`	The object reference to the Encode object, used to convert the source code to Unicode.
`${^OPEN}`	Internal use: \0 separated Input / Output layer information.
`${^UNICODE}`	Read-only Unicode settings.
`${^UTF8CACHE}`	State of the internal UTF-8 offset caching code.
`${^UTF8LOCALE}`	Indicates whether UTF8 locale was detected at startup.
IO and Separators
`$.`	Current line number (or record number) of most recent filehandle.
`$/`	Input record separator.
`$\|`	Output autoflush. 1=autoflush, 0=default. Applies to currently selected handle.
`$,`	Output field separator (lists)
`$\`	Output record separator.
`$"`	Output list separator. (interpolated lists)
`$;`	Subscript separator. (Use a real multidimensional array instead.)
Formats
`$%`	Page number for currently selected output channel.
`$=`	Current page length.
`$-`	Number of lines left on page.
`$~`	Format name.
`$^`	Name of top-of-page format.
`$:`	Format line break characters
`$^L`	Form feed (default "\f").
`$^A`	Format Accumulator
Status Reporting
`$?`	Child error. Status code of most recent system call or pipe.
`$!`	Operating System Error. (What just went 'bang'?)
`%!`	Error number hash
`$^E`	Extended Operating System Error (Extra error explanation).
`$@`	Eval error.
`${^CHILD_ERROR_NATIVE}`	Native status returned by the last pipe close, backtick (`` ) command, successful call to wait() or waitpid(), or from the system() operator.
ID's and Process Information
`$$`	Process ID
`$<`	Real user id of process.
`$>`	Effective user id of process.
`$(`	Real group id of process.
`$)`	Effective group id of process.
`$0`	Program name.
`$^O`	Operating System name.
Perl Status Info
`$]`	Old: Version and patch number of perl interpreter. Deprecated.
`$^C`	Current value of flag associated with -c switch.
`$^D`	Current value of debugging flags
`$^F`	Maximum system file descriptor.
`$^I`	Value of the -i (inplace edit) switch.
`$^M`	Emergency Memory pool.
`$^P`	Internal variable for debugging support.
`$^R`	Last regexp (?{code}) result.
`$^S`	Exceptions being caught. (eval)
`$^T`	Base time of program start.
`$^V`	Perl version.
`$^W`	Status of -w switch
`${^WARNING_BITS}`	Current set of warning checks enabled by `use warnings;`
`$^X`	Perl executable name.
`${^GLOBAL_PHASE}`	Current phase of the Perl interpreter.
`$^H`	Internal use only: Hook into Lexical Scoping.
`%^H`	Internaluse only: Useful to implement scoped pragmas.
`${^TAINT}`	Taint mode read-only flag.
`${^WIN32_SLOPPY_STAT}`	If true on Windows `stat()` won't try to open the file.
Command Line Args
`ARGV`	Filehandle iterates over files from command line (see also `<>`).
`$ARGV`	Name of current file when reading <>
`@ARGV`	List of command line args.
`ARGVOUT`	Output filehandle for -i switch
Miscellaneous
`@F`	Autosplit (-a mode) recipient.
`@INC`	List of library paths.
`%INC`	Keys are filenames, values are paths to modules included via `use, require,` or `do`.
`%ENV`	Hash containing current environment variables
`%SIG`	Signal handlers.
`$[`	Array and substr first element (Deprecated!).

See perlvar for detailed descriptions of each of these (and a few more) special variables.

Consed--A Finishing Package (BAM File Viewer, Assembly Editor, Autofinish, Autoreport, Autoedit, and Align Reads To Reference Sequence)

Neel — Fri, 07 Feb 2020 07:16:22 -0600

Supports Illumina, 454, other Next-Gen and Sanger Reads and allows mixtures of these read types
Consed includes BamScape which can view bam files with unlimited numbers of reads. BamScape can bring up consed to edit reads and the reference sequence in targeted regions.
Consed is compatible with Newbler, Cross_match, Phrap, MIRA, Velvet and PCAP output.
Quickly takes the user to each variant site for viewing (also available as an automated report)
Overview of assembly can help detect and fix misassemblies
Editing time reduced by the program's ability to pin-point problem areas
Editing is guided by error probabilities

Address of the bookmark: http://www.phrap.org/consed/consed.html

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

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

What Junk DNA? It’s an Operating System

Rahul Agarwal — Mon, 19 Aug 2013 15:24:26 -0500

The report adds to growing experimental support for the idea that all that extra stuff in the human genes, once referred to as “junk DNA,” is more than functionless, space-filling material that happens to make up nearly 98% of the genome. The paper adds to a growing body of knowledge establishing a considerable role for this material in the regulation of gene expression and its potential role in human disease.

Address of the bookmark: http://www.genengnews.com/keywordsandtools/print/3/32115/

RNA Sequencing Helps Identify Functional Variants from GWAS

Andaleeb — Fri, 22 Nov 2013 21:33:33 -0600

For Alzheimer’s and other complex disorders, mining the genome for disease-associated variants is no longer the obstacle. The challenge nowadays is figuring out how the identified loci relate to disease. As reported last month in Nature and its associated journals, advances in high-throughput RNA sequencing are providing new tools for understanding how disease loci influence gene expression—a starting point for understanding their connection to pathogenesis.

Address of the bookmark: http://schizophreniaforum.org/new/detail.asp?id=1953

Pathway Analysis

Rahul Agarwal — Fri, 03 Oct 2014 08:51:13 -0500

Pathway Analysis is usually performed with aim to enrich the genes with their functional information and reveal the underlying biological mechanisms pursue by genes. Pathway Analysis is not only limited to what biological pathways a particular set of expressed genes follow but also to disclose the relationships between these genes. With availability of more genomics, transcriptomics and proteomics data, interactions between genes involve in multiple pathways become more clear and also relationships between the genes, their transcripts, and their gene products. However, existing tools and dbs mainly based on knowledge driven approach in which pathways will be identified by finding the correlation between the information in one of the pathway knowledge databases (KEGG,Reactome,Panther,BioCarta, Panther,GO,NCI,WikiPathways,etc) and gene expression result for a specific conditions for instance tumor, obesity , cold resistant crops/plants, etc.

Introductory Articles/ppt/sources:

http://www.ploscompbiol.org/article/info%3Adoi%2F10.1371%2Fjournal.pcbi.1002375

http://bioinformatics.mdanderson.org/MicroarrayCourse/Lectures09/Pathway%20Analysis.pdf

http://gettinggeneticsdone.blogspot.de/2012/03/pathway-analysis-for-high-throughput.html

http://davetang.org/muse/tag/pathway/

https://www.biostars.org/p/42219/

http://bioinformatics.ca//files/public/Pathways_2014_Module4_v2.pdf

http://bioinformatics.ca//files/public/Pathways_2014_Module2.pdf

Impotant Database and Tools:

GeneMANIA, Cytoscape, IPA and Metacore (Commerical ), Pathway Commons, Reactome ,Panther, BioCyc, WikiPathways, Pathvisio, KEGG, NCI, Stringdb, Amigo, WebGestalt ,ConsensusPathDB ,GSEA,Blast2go

Popular R based tools:

Reactome.db, ReactomePA, ClusterProfiler, Gage, SPIA, topGO, Pathview,DOSE,GOStat

More:

http://www.bioconductor.org/help/search/index.html?q=Enrichment+analysis+

Hoffman Lab

Tue, 12 Jan 2016 02:47:41 -0600

They develop machine learning techniques to better understand chromatin biology. These models and algorithms transform high-dimensional functional genomics data into interpretable patterns and lead to new biological insight.

https://www.pmgenomics.ca/hoffmanlab/

Salzberg lab

Mon, 15 May 2017 05:14:01 -0500

We are a computational biology lab that develops novel methods for analysis of DNA and RNA sequences. Our research includes software for aligning and assembling RNA-seq data, whole-genome assembly, and microbiome analysis. We work closely with biomedical scientists to apply these methods to current problems arising in a broad spectrum of biological and medical research areas. We’re also part of the Center for Computational Biology, a group of 20+ faculty members and their labs at Johns Hopkins working on computational, statistical, and mathematical methods that can turn massive genomic data sets into biologically and clinically useful information.

https://salzberg-lab.org/