Published October 5, 2018 | Version 2.013 (alpha)

vcfSampleCompare 2 (alpha) - Sort and filter variants whose states differ between sample groups

Authors/Creators

  • 1. Princeton University

Description

This utility sorts and (optionally) filters the rows/variants of a VCF file (containing data for 2 or more samples) based on the differences in the variant data between samples or sample groups. Degree of "difference" is determined by either the best possible degree of separation of sample groups by genotype calls or the difference in average allelic frequency of each sample or sample group (with a gap size threshold). The pair of samples or sample groups used to represent the difference for a variant row is the one leading to the greatest difference in consistent genotype or average allelic frequencies (i.e. observation ratios, e.g. AO/DP) of the same variant state. If sample groups are not specified, the pair of samples leading to the greatest difference is greedily discovered and chosen to represent the variant/row.

Currently, there are 2 separate output files and filtered lines of the VCF file are omitted. The next version will instead incorporate sample comparison data into the VCF format and set the filter column value.

This incremental release fixes an issue with the handling of how no data is represented among samples combined using `bcftools merge` when the number of comma-delimited ALT states differs. Previous versions of this tool die with a fatal error in this case.

Files

Files (1.4 MB)

Name Size Download all
md5:44699b072543d3ac633bf6cae28ad907
435.8 kB Download
md5:55ee2a35dbffb5cb8e52bb4bc55b4ee7
436.0 kB Download
md5:1af0dc8b17f5a8a67d38e7c33110496b
553.7 kB Download