An automated workflow for single nucleotide polymorphisms detection is developed. This is known as OVarFlow . This workflow also helps in the identification of insertions and deletions.
OVarFlow provides a wide range of applications in sequence annotation of model and non-model organisms . The workflow is based on GATK best practices. It makes the task of variant calling easy by reducing the computation time and resources.
A reference genome is provided as input in FASTA format along with reference annotation in a general feature format (.gff.gz). After that, second-generation sequencing data is provided in fastq format. Using this data, functionally annotated variants are generated. The workflow also consists of an extended workflow that generates variants by base quality score recalibration (BQSR).
OVarFlow is written in Python & Snakemake, and Bash languages and is available for Linux on GitLab.
For more information, read here.
- Bathke, J., Lühken, G. (2021). OVarFlow: a resource optimized GATK 4 based Open source Variant calling workFlow. BMC Bioinformatics 22, 402.
How to execute matlab from terminal in Ubuntu (Linux)?
How to install Kpax on Ubuntu (Linux)?
Kpax is a bioinformatics program to search and align protein structures . It is currently available for Linux platforms only. In this article, we are going to install the latest version of Kpax (5.1.3) on Ubuntu (Linux). (more…)
How to run do_dssp command (mkdssp) in Gromacs 2022?
In the latest version of GROMACS (2022) , there are some issues regarding the gmx do_dssp command. Apparently, this command either does not run displaying a fatal error, or if it runs then it does not read any frame from MD simulation files. In this article, we are going to run the same command for GROMACS 2022. (more…)