Abstract:Currently, there are many publicly available Next Generation Sequencing tools developed for variant annotation and classification. However, as modern sequencing technology produces more and more sequencing data, a more efficient analysis program is desired, especially for variant analysis. In this study, we updated SNPAAMapper, a variant annotation pipeline by converting perl codes to python for generating annotation output with an improved computational efficiency and updated information for broader applicability. The new pipeline written in Python can classify variants by region (Coding Sequence, Untranslated Regions, upstream, downstream, intron), predict amino acid change type (missense, nonsense, etc.), and prioritize mutation effects (e.g., synonymous > non-synonymous) while being faster and more efficient. Our new pipeline works in five steps. First, exon annotation files are generated. Next, the exon annotation files are processed, and gene mapping and feature information files are produced. Afterward, the python scrips classify the variants based on genomic regions and predict the amino acid change category. Lastly, another python script prioritizes and ranks the mutation effects of variants to output the result file. The Python version of SNPAAMapper accomplished the overall speed by running most annotation steps in a substantially shorter time. The Python script can classify variants by region in 53 s compared to 166 s for the Perl script in a test sample run on a Latitude 7480 Desktop computer with 8GB RAM and an Intel Core i5-6300 CPU @ 2.4Ghz. Steps of predicting amino acid change type and prioritizing mutation effects of variants were executed within 1 s for both pipelines. SNPAAMapper-Python was developed and tested on the ClinVar database, a NCBI database of information on genomic variation and its relationship to human health. We believe our developed Python version of SNPAAMapper variant annotation pipeline will benefit the community by elucidating the variant consequence and speed up the discovery of causative genetic variants through whole genome/exome sequencing. Source codes, test data files, instructions, and further explanations are available on the web at https://github.com/BaiLab/SNPAAMapper-Python.

Ingap: an Integrated Next-Generation Genome Analysis Pipeline

mInDel: a high-throughput and efficient pipeline for genome-wide InDel marker development

Ingap-Sv: a Novel Scheme to Identify and Visualize Structural Variation from Paired End Mapping Data

Pseudo-Sanger Sequencing: Massively Parallel Production of Long and Near Error-Free Reads Using NGS Technology

Ingap-Family:Accurate Detection of Meiotic Recombination Loci and Causal Mutations by Filtering out Artificial Variants Due to Genome Complexities

DNAscan: a fast, computationally and memory efficient bioinformatics pipeline for the analysis of DNA next-generation-sequencing data

GAPIT Version 2: an Enhanced Integrated Tool for Genomic Association and Prediction

A De-Novo Genome Analysis Pipeline (DeNoGAP) for large-scale comparative prokaryotic genomics studies

A Novel Approach to Detect Large Indels from Targeted Sequencing Data in Clinical Cancer Setting

SNPAAMapper-Python: A highly efficient genome-wide SNP variant analysis pipeline for Next-Generation Sequencing data

Metapipeline-DNA: A Comprehensive Germline & Somatic Genomics Nextflow Pipeline

SPANDx: a genomics pipeline for comparative analysis of large haploid whole genome re-sequencing datasets

PAPipe: a pipeline for comprehensive population genetic analysis

GMAP and GSNAP for Genomic Sequence Alignment: Enhancements to Speed, Accuracy, and Functionality

Read Annotation Pipeline for High-Throughput Sequencing Data.

Please Mind the Gap: Indel-Aware Parsimony for Fast and Accurate Ancestral Sequence Reconstruction and Multiple Sequence Alignment including Long Indels

NCBI prokaryotic genome annotation pipeline

InteMAP: Integrated metagenomic assembly pipeline for NGS short reads

TGS-GapCloser: A fast and accurate gap closer for large genomes with low coverage of error-prone long reads

A Fully Integrated End-to-End Genome Analysis Accelerator for Next-Generation Sequencing