iscanvcfmerge
Python tool to merge cross-species Illumina iScan genotype data with a reference set of genotypes from a pre-existing source.
Science Score: 23.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
○CITATION.cff file
-
○codemeta.json file
-
○.zenodo.json file
-
✓DOI references
Found 3 DOI reference(s) in README -
✓Academic publication links
Links to: ncbi.nlm.nih.gov -
○Committers with academic emails
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (14.9%) to scientific vocabulary
Keywords
Repository
Python tool to merge cross-species Illumina iScan genotype data with a reference set of genotypes from a pre-existing source.
Basic Info
Statistics
- Stars: 0
- Watchers: 0
- Forks: 0
- Open Issues: 0
- Releases: 0
Topics
Metadata Files
README.md
iScanVCFMerge
iScanVCFMerge is a Python tool to facilitate the cross-species application of Illumina iScan system microarrays. The tool merges VCF genotypes exported from GenomeStudio with a second VCF, comprising genotypes derived from other samples and sources. Merging is based on matches of chromosome, position and certain conditions of major and minor alleles, with matched rows from each VCF concatenated into a single row (comprising all individuals) in the output files. The full algorithm is explained in the accompanying manuscript, where we reported use of the human Infinium Multi-Ethnic Global and Infinium Omni 2.5 arrays (hg19) to genotype great apes, and merged those with the genotypes of conspecifics previously published elsewhere.
What's new since the paper came out?
Version 1.2
- Bug fixed to properly output VCF header, which had spurious line breaks on some Python versions.
- Output VCF contigs are now sorted in the order of the sequence dictionary, as pulled from the
. - The FORMAT column is restored with a GT value; this was causing output VCFs to fail GATK ValidateVariants.
- Upgraded deprecated pandas functions.
Version 1.1
- Bugs fixed to properly handle some multi-allelic sites.
- The reference population VCF file must now be bgzipped and indexed with tabix. This requirement does not apply to the iScan VCF file, which can either be uncompressed or gzip compressed.
- In the prior version, the complete reference population VCF file was read into memory before the relevant records were pulled. This caused issues for some users handling enormous reference VCF files. In this version, we use the Pysam library's lightweight wrapper of the htslib C-API to pull only the relevant records in the first place. The script should now run near-instantaneously, irrespective of input file size.
- Console output is now handled by the Python logging module and is written to a .log file in the output directory.
- Version numbering now follows 1.x vs 0.x format.
Installation
iScanVCFMerge 1.2 requires Python 3.9. It has been successfully tested on MacOS Monterey 12.3.1 and CentOS Linux 7.9.2009.
Install the required packages first, if needed:
python3 -m pip install pandas pysam
Clone from Github and execute:
git clone "https://github.com/baneslab/iScanVCFMerge.git"
cd iScanVCFMerge
python3 iScanVCFMerge.py
Please note that the PyPi package is no longer maintained.
Usage
iScanVCFMerge [-h] -I <iScan_vcf> -R <reference_vcf> -O <output_directory>
Optional arguments:
-h, --help Show the help message
-I, --iScanVCF Path to your iScan VCF file (.vcf or .vcf.gz)
-R, --ReferenceVCF Path to your reference VCF file, with which the iScan file will be merged. This must be bgzip compressed and be indexed with tabix
-O, --output_directory Name of the output directory (will be created if it doesn't exist)
Citation
Please cite the use of this software as follows:
Fountain, E. D., Zhou, L-C., Karklus, A., Liu, Q-X., Meyers, J., Fontanilla, I. K., Rafael, E. F., Yu, J-Y., Zhang, Q., Zhu, X-L., Pei, E-L., Yuan, Y-H. and Banes, G. L. (2021). Cross-species application of Illumina iScan microarrays for cost-effective, high-throughput SNP discovery. Frontiers in Ecology and Evolution, 9:629252, doi: 10.3389/fevo.2021.629252.
The Research Resource Identifier for iScanVCFMerge is RRID:SCR_021193.
GitHub Events
Total
Last Year
Issues and Pull Requests
Last synced: 6 months ago
All Time
- Total issues: 4
- Total pull requests: 2
- Average time to close issues: 2 months
- Average time to close pull requests: 4 days
- Total issue authors: 1
- Total pull request authors: 1
- Average comments per issue: 0.25
- Average comments per pull request: 0.0
- Merged pull requests: 2
- Bot issues: 0
- Bot pull requests: 0
Past Year
- Issues: 0
- Pull requests: 0
- Average time to close issues: N/A
- Average time to close pull requests: N/A
- Issue authors: 0
- Pull request authors: 0
- Average comments per issue: 0
- Average comments per pull request: 0
- Merged pull requests: 0
- Bot issues: 0
- Bot pull requests: 0
Top Authors
Issue Authors
- grahamlbanes (4)
Pull Request Authors
- grahamlbanes (2)
Top Labels
Issue Labels
Pull Request Labels
Packages
- Total packages: 1
-
Total downloads:
- pypi 11 last-month
- Total dependent packages: 0
- Total dependent repositories: 1
- Total versions: 2
- Total maintainers: 1
pypi.org: iscanvcfmerge
Python tool to merge cross-species Illumina iScan genotype data with a reference set of genotypes from a pre-existing source.
- Homepage: https://github.com/baneslab/iScanVCFMerge
- Documentation: https://iscanvcfmerge.readthedocs.io/
- License: MIT
-
Latest release: 1.1
published over 4 years ago