LCVCFtools is a simple C++ program designed for working with VCF v4.2 files generated from low-coverage whole genome sequencing. This tools is useful for data optimization before genotype imputation steps by removing samples or variantes with excessive missing data.
You can easily compile this program using LCVCFtools.pro file and qmake tool, just run the following command:
sudo apt install qt5-qmake libboost-all-dev libz-dev
git clone https://github.com/marcusnizalvarez/LCVCFtools.git
cd LCVCFtools/
qmake && make
- --vcf Read from VCF file. Use - to read from stdin.
- --gzvcf Read from Gzip compressed VCF file. Use - to read from stdin.
- --minGQ Minimum genotype quality in PhredScale. [Default=20]
- --minDP Minimum depth. [Default=5]
- --MAF Minor allele frequency, based on allele depth (AD). [Default=0.1]
- --minGCR Minimum genotype call rate [Default=0].
- --minDPR Minimum DP rate. Can be defined multiple times.
- --minGQR Minimum GQ rate. Can be defined multiple times.
- --remove Remove samples listed in a file.
- --keep Keep samples listed in a file, after --remove.
- --sample-stats Output sample statistics to 'stats.tsv'.
- --keep-multiallelic Don't skip multiallelic variants.
- --ID Generate generic ID, useful for programs like Plink.
- --verbose Verbose mode.
- --help Print this message.
./LCVCFtools --gzvcf example.vcf.gz - --minGQ 20 --minDP 5 --minGCR 0.25 --minDPR 5 0.5 --MAF 0.1 --sample-stats | gzip -c > output.vcf.gz
Cite as: Alvarez MVN. LCVCFtools v1.0.2‑alpha. 2022. https://doi.org/10.5281/zenodo.5259931.
This software was originally developed for this paper: Alvarez, MVN. et al. Nyssorhynchus darlingi genome-wide studies related to microgeographic dispersion and blood-seeking behavior. Parasites & Vectors. 2022. 15(1):106.
Author: Marcus Vinicius Niz Alvarez
Email: marcus.alvarez@unesp.br
São Paulo State University, UNESP - Biotechnology Institute and Bioscience Institute, Botucatu, 18618-689, Brazil.
LCVCFtools is under GNU GPLv3.0 license and Boost Library is under Boost Software License v1.0.