BOL: Related items

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.

Painless package development for R

Abhi — Tue, 03 May 2016 05:31:06 -0500

Devtools makes package development a breeze: it works with R’s existing conventions for code structure, adding efficient tools to support the cycle of package development. With devtools, developing a package becomes so easy that it will be your default layout whenever you’re writing a significant amount of code.

Before you get started be sure to check out:

Address of the bookmark: https://www.rstudio.com/products/rpackages/devtools/

SOAP2 : Short Oligonucleotide Analysis Package 2

Poonam Mahapatra — Wed, 23 May 2018 10:09:22 -0500

SOAPaligner/soap2 is a member of the SOAP (Short Oligonucleotide Analysis Package). It is an updated version of SOAP software for short oligonucleotide alignment. The new program features in super fast and accurate alignment for huge amounts of short reads generated by Illumina/Solexa Genome Analyzer. Compared to soap v1, it is one order of magnitude faster. It require only 2 minutes aligning one million single-end reads onto the human reference genome. Another remarkable improvement of SOAPaligner is that it now supports a wide range of the read length. SOAPaligner benefitted in time and space efficiency by a revolution in the basic data structures and algorithms used.The core algorithms and the indexing data structures (2way-BWT) are developed by the algorithms research group of the Department of Computer Science, the University of Hong Kong (T.W. Lam, Alan Tam, Simon Wong, Edward Wu and S.M. Yiu).

Address of the bookmark: http://soap.genomics.org.cn/soapaligner.html

Genobuntu: Package for Next Generation Sequencing and Genome Assembly

BioStar — Mon, 18 May 2020 16:47:56 -0500

Genobuntu is a software package containing more than 70 software and packages oriented towards NGS. In its current version, Genobuntu supports pre assembly tools, genome assemblers as well as post assembly tools.

Commonly used biological software and example script files for different assembly pipelines have also been provided, where the example script files can be updated to suit one’s experimental needs. Genobuntu attempts to reduce the amount of time and energy needed to build software workstations and it can also act as a good teaching source for a class room setting.

Therefore, Genobuntu offers a well-tailored environment for both novices and experts working in the field of genome assembly.

Features

Velvet
MiB
SSAKE
EULER
VCAKE
ABySS
ALLPATHS
Celera
SHARCGS
Allpaths
IDBA
TAIPAN
Edena
SOAPdenovo
Maq
IDBA-UD
No. of Reads present in the Ref. Seq.
ART NGS Reads Simulator
HiTEC, FASTQC
Minimum Description Length
SOAPaligner
Sequencing Read Archive Toolkit

Address of the bookmark: https://sourceforge.net/projects/genobuntu/

scikit-bio™ is an open-source, BSD-licensed, python package providing data structures, algorithms, and educational resources for bioinformatics.

Jit — Fri, 02 Mar 2018 04:29:47 -0600

scikit-bio is currently in beta. We are very actively developing it, and backward-incompatible interface changes can and will arise. To avoid these types of changes being a surprise to our users, our public APIs are decorated to make it clear to users when an API can be relied upon (stable) and when it may be subject to change (experimental). See the API stability docs for more details, including what we mean by stable and experimental in this context.

Address of the bookmark: http://scikit-bio.org/

GeneOverlap: An R package to test and visualize gene overlaps

Jit — Thu, 27 Dec 2018 19:45:52 -0600

Overlapping gene lists can reveal biological meanings and may lead to novel hypotheses. For example, histone modification is an important cellular mechanism that can pack and re-pack chromatin. By making the chromatin structure more dense or loose, the gene expression can be turned on or off. Tri-methylation on lysine 4 of histone H3 (H3K4me3) is associated with gene activation and its genome-wide enrichment can be mapped by using ChIP-seq experiments. Because of its activating role, if we overlap the genes that are bound by H3K4me3 with the genes that are highly expressed, we should expect a positive association. Similary, we can perform such kind of overlapping between the gene lists of different histone modifications with that of various expression groups and establish each histone modification’s role in gene regulation.

Address of the bookmark: https://bioconductor.org/packages/release/bioc/vignettes/GeneOverlap/inst/doc/GeneOverlap.pdf

IgBLAST: a popular NCBI package for classifying and analyzing immunoglobulin (IG) and T cell receptor (TCR) variable domain sequences

BioJoker — Thu, 23 Jan 2020 11:34:37 -0600

NCBI team released a new version of IgBLAST with four new improvements. IgBLAST is a popular NCBI package for classifying and analyzing immunoglobulin (IG) and T cell receptor (TCR) variable domain sequences. Improvements are:

1. Support for the new FWR4 annotation feature in the AIRR format, both in standard format and in the AIRR alignment format.

2. The previous “-penalty” parameter was renamed as -V_penalty to be consistent with other IgBLAST penalty options.

3. Restored constant internal BLAST search parameters for domain annotation (i.e., FWR/CDR) such that this process is not influenced by user parameters.

4. Corrected FWR/CDR annotations for certain mouse VK and rat VH germline genes.

IgBLAST 1.15.0 is available for download from the BLAST FTP area. See the the new manual on GitHub for information about setting up and running IgBLAST.

If you have any questions or concerns, please contact blast-help@ncbi.nlm.nih.gov

PheWAS: R package is designed to provide an accessible interface to the phenome wide association study

BioStar — Thu, 30 Jul 2020 22:06:37 -0500

The PheWAS R package is designed to provide an accessible interface to the phenome wide association study. For a description of the methods available and some simple examples, please see the package vignette or the R documentation. For installation help, see below. ##Installing the PheWAS Package The PheWAS package can be installed using the devtools package. The following code when executed in R will get you started:

install.packages("devtools")
#It may be necessary to install required as not all package dependencies are installed by devtools:
install.packages(c("dplyr","tidyr","ggplot2","MASS","meta","ggrepel","DT"))
devtools::install_github("PheWAS/PheWAS")
library(PheWAS)

Address of the bookmark: https://github.com/PheWAS/PheWAS

Luigi: a Python package that helps you build complex pipelines of batch jobs.

Neel — Thu, 24 Jun 2021 05:43:31 -0500

Luigi is a Python (3.6, 3.7, 3.8, 3.9 tested) package that helps you build complex pipelines of batch jobs. It handles dependency resolution, workflow management, visualization, handling failures, command line integration, and much more.

Run pip install luigi to install the latest stable version from PyPI. Documentation for the latest release is hosted on readthedocs.

Run pip install luigi[toml] to install Luigi with TOML-based configs support.

Address of the bookmark: https://github.com/spotify/luigi