BOL: Surabhi Chaudhary's bookmarks

VCF2PopTree: a client-side software to construct population phylogeny from genome-wide SNPs

Surabhi Chaudhary — Sat, 25 Dec 2021 00:13:31 -0600

VCF2PopTree is a client-side software written in Javascript and it runs purely within the user’s computer/browser. VCF2PopTree is compatible with all population browsers including Chrome, Opera, Edge and Firefox and works equally efficient in Mac, Windows and Linux (Ubuntu).

Furthermore, it displays the tree in a mobile phone (iPhone and Android) if the input file size is small. CITATION: Subramanian, S., Ramasamy, U. and Chen, D. (2019). VCF2PopTree: a client-side software to construct population phylogeny from genome-wide SNPs. Peer J. x:yy.

Address of the bookmark: https://github.com/sansubs/vcf2pop

vcf2maf convert !

Surabhi Chaudhary — Fri, 17 Dec 2021 03:20:01 -0600

To convert a VCF into a MAF, each variant must be mapped to only one of all possible gene transcripts/isoforms that it might affect. But even within a single isoform, a Missense_Mutation close enough to a Splice_Site, can be labeled as either in MAF format, but not as both. This selection of a single effect per variant, is often subjective. And that's what this project attempts to standardize. The vcf2maf and maf2maf scripts leave most of that responsibility to Ensembl's VEP, but allows you to override their "canonical" isoforms, or use a custom ExAC VCF for annotation. Though the most useful feature is the extensive support in parsing a wide range of crappy MAF-like or VCF-like formats we've seen out in the wild.

Address of the bookmark: https://github.com/mskcc/vcf2maf

maftools

Surabhi Chaudhary — Fri, 17 Dec 2021 03:18:28 -0600

With advances in Cancer Genomics, Mutation Annotation Format (MAF) is being widely accepted and used to store somatic variants detected. The Cancer Genome Atlas Project has sequenced over 30 different cancers with sample size of each cancer type being over 200. Resulting data consisting of somatic variants are stored in the form of Mutation Annotation Format. This package attempts to summarize, analyze, annotate and visualize MAF files in an efficient manner from either TCGA sources or any in-house studies as long as the data is in MAF format.

https://www.bioconductor.org/packages/devel/bioc/vignettes/maftools/inst/doc/maftools.html

Address of the bookmark: https://github.com/PoisonAlien/maftools

Pandoc: a universal document converter

Surabhi Chaudhary — Thu, 24 Jun 2021 01:33:47 -0500

If you need to convert files from one markup format into another, pandoc is your swiss-army knife. Pandoc can convert almost all formats

https://pandoc.org/index.html

Address of the bookmark: https://pandoc.org/

gggenomes: A grammar of graphics for comparative genomics

Surabhi Chaudhary — Mon, 01 Feb 2021 14:47:32 -0600

gggenomes is a versatile graphics package for comparative genomics. It extends the popular R visualization packageggplot2 by adding dedicated plot functions for genes, syntenic regions, etc. and verbs to manipulate the plot to, for example, quickly zoom in into gene neighborhoods.

Address of the bookmark: https://github.com/thackl/gggenomes

Phytozome v12.1: plant science community hub for accessing palnts genomic data

Surabhi Chaudhary — Tue, 17 Mar 2020 07:30:17 -0500

Phytozome, the Plant Comparative Genomics portal of the Department of Energy's Joint Genome Institute, provides JGI users and the broader plant science community a hub for accessing, visualizing and analyzing JGI-sequenced plant genomes, as well as selected genomes and datasets that have been sequenced elsewhere. As of release v12.1.6, Phytozome hosts 93 assembled and annotated genomes, from 82 Viridiplantae species. More than half of these genomes have been sequenced, assembled and/or annotated with JGI Plant Science program resources. By integrating this large collection of plant genomes into a single resource and performing comprehensive and uniform annotation and analyses, Phytozome facilitates accurate and insightful comparative genomics studies.

Address of the bookmark: https://phytozome.jgi.doe.gov/pz/portal.html

RAINBOWR: Reliable Association INference By Optimizing Weights with R (R package for SNP-set GWAS and multi-kernel mixed model)

Surabhi Chaudhary — Fri, 28 Feb 2020 23:27:37 -0600

RAINBOWR(Reliable Association INference By Optimizing Weights with R) is a package to perform several types of GWAS as follows.

Single-SNP GWAS with RGWAS.normal function
SNP-set (or gene set) GWAS with RGWAS.multisnp function (which tests multiple SNPs at the same time)
Check epistatic (SNP-set x SNP-set interaction) effects with RGWAS.epistasis (very slow and less reliable)

https://github.com/KosukeHamazaki/RAINBOWR

https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1007663

https://cran.r-project.org/web/packages/RAINBOWR/index.html

Address of the bookmark: https://github.com/KosukeHamazaki/RAINBOWR

PRICE (Paired-Read Iterative Contig Extension), a de novo genome assembler implemented in C++.

Surabhi Chaudhary — Mon, 11 Jun 2018 03:08:26 -0500

We are pleased to release PRICE (Paired-Read Iterative Contig Extension), a de novo genome assembler implemented in C++. Its name describes the strategy that it implements for genome assembly: PRICE uses paired-read information to iteratively increase the size of existing contigs. Initially, those contigs can be individual reads from a subset of the paired-read dataset, non-paired reads from sequencing technologies that provide non-paired data, or contigs that were output from a prior run of PRICE or any other assembler. http://derisilab.ucsf.edu/software/price/

Address of the bookmark: http://derisilab.ucsf.edu/software/price/

SEX-DETector: A Probabilistic Approach to Study Sex Chromosomes in Non-Model Organisms

Surabhi Chaudhary — Wed, 30 May 2018 15:57:31 -0500

SEX-DETector is a probabilistic method that relies on RNAseq data from a cross (parents and progeny of each sex) to infer autosomal and sex-linked genes (genes located on the non recombining part of sex chromosomes).

How does SEX-DETector work?

SEX-DETector does not require prior sequencing of a reference genome: the same sequencing data can be used for the assembly and for the mapping of the reads. A full documentation on the pipeline can be found here.

we recommend Trinity for the assembly.
Trinity components should be merged with cap3. Our code to perform the merging is available here.
We recommend BWA for mapping of the reads.
When the mapping has been perfomed, the individuals need to be genotyped; SEX-DETector takes files produced by Reads2snp (which is available for download on the PopPhyl website) as input.

Address of the bookmark: http://lbbe.univ-lyon1.fr/-SEX-DETector-.html?lang=eg

%MinMax: A versatile tool for calculating and comparing synonymous codon usage and its impact on protein folding.

Surabhi Chaudhary — Thu, 24 May 2018 02:53:31 -0500

%MM calculates whether a given gene sequence encodes amino acids using the most common codons possible, the least common codons possible, or (most typically) some combination of these extremes. See our PLoS ONE paper for more details on how the %MinMax algorithm works. %MinMax results are averaged over an 18-codon sliding window; hence the result for "codon window = 1" is the average codon usage for codons 1-18, codon window 2 = codons 2-19, etc.

Address of the bookmark: http://www.codons.org/