1# vcflib
2
3### A C++ library for parsing and manipulating VCF files.
4
5![Github-CI](https://github.com/vcflib/vcflib/workflows/CI/badge.svg) [![Travis-CI](https://travis-ci.org/vcflib/vcflib.svg?branch=master)](https://travis-ci.org/vcflib/vcflib) [![AnacondaBadge](https://anaconda.org/bioconda/vcflib/badges/installer/conda.svg)](https://anaconda.org/bioconda/vcflib) [![DL](https://anaconda.org/bioconda/vcflib/badges/downloads.svg)](https://anaconda.org/bioconda/vcflib) [![BrewBadge](https://img.shields.io/badge/%F0%9F%8D%BAbrew-vcflib-brightgreen.svg)](https://github.com/brewsci/homebrew-bio) [![GuixBadge](https://img.shields.io/badge/gnuguix-vcflib-brightgreen.svg)](https://www.gnu.org/software/guix/packages/V/) [![DebianBadge](https://badges.debian.net/badges/debian/testing/libvcflib-dev/version.svg)](https://packages.debian.org/testing/libvcflib-dev)
6[![C++0x](https://img.shields.io/badge/Language-C++0x-steelblue.svg)](https://www.cprogramming.com/c++11/what-is-c++0x.html)
7[![Gitter](https://badges.gitter.im/ekg/vcflib.svg)](https://gitter.im/ekg/vcflib?utm_source=badge&utm_medium=badge&utm_campaign=pr-badge&utm_content=badge)
8
9## overview
10
11The [Variant Call Format
12(VCF)](http://www.1000genomes.org/wiki/Analysis/Variant%20Call%20Format/vcf-variant-call-format-version-41)
13is a flat-file, tab-delimited textual format that
14describes reference-indexed variations between individuals.  VCF
15provides a common interchange format for the description of variation
16in individuals and populations of samples, and has become the
17*de facto* standard reporting format for a wide array of genomic
18variant detectors.
19
20vcflib provides methods to manipulate and interpret sequence variation
21described by VCF.  It is both:
22
23 * an API for parsing and operating on records of genomic variation as
24   it can be described by the VCF format
25 * a collection of command-line utilities for executing complex
26   manipulations on VCF files
27
28vclib is both a library (with an API) and a collection of useful
29tools. The API provides a quick and extremely permissive method to
30read and write VCF files.  Extensions and applications of the library
31provided in the included utilities (*.cpp) comprise the vast bulk of
32the library's utility.
33
34---
35
36Short index:
37
38- [Install](#INSTALL)
39- [Usage](#USAGE)
40- [TOOLS](#TOOLS)
41  * [Filter](#filter)
42  * [Metrics](#metrics)
43  * [Phenotype](#phenotype)
44  * [Genotype](#genotype)
45  * [Transformation](#transformation)
46  * [Statistics](#statistics)
47  * [Scripts](#scripts)
48- [Link library](#link-library)
49- [Build from source](#build-from-source)
50- [Development](#Development)
51- [LICENSE](#LICENSE)
52
53---
54
55## INSTALL
56
57### [Bioconda](https://bioconda.github.io/user/install.html)
58
59Conda installs in user land without root access
60
61```sh
62conda install -c bioconda vcflib
63```
64
65### [Homebrew](https://brew.sh)
66
67Homebrew installs on Linux and Mac OSX
68
69```sh
70brew install brewsci/bio/vcflib
71```
72
73### [Debian](https://debian.org/)
74
75For Debian and Ubuntu
76
77```sh
78apt-get install libvcflib-tools libvcflib-dev
79```
80
81### [GNU Guix](https://guix.gnu.org/)
82
83We develop against guix
84
85```sh
86guix package -i vcflib
87```
88
89## USAGE
90
91Users are encouraged to drive the utilities in the library in a
92streaming fashion, using Unix pipes to fully utilize resources on
93multi-core systems.  Piping provides a convenient method to interface
94with other libraries (vcf-tools, BedTools, GATK, htslib,
95[bio-vcf](https://github.com/vcflib/bio-vcf), bcftools,
96[freebayes](https://github.com/freebayes)) which interface via VCF
97files, allowing the composition of an immense variety of processing
98functions. Examples can be found in the scripts,
99e.g. [script](./scripts/vcfgtcompare.sh).
100
101
102# TOOLS
103
104<!--
105
106    >>> from pytest.rtest import run_stdout, head, cat
107
108-->
109
110
111<!--
112  Created with ./scripts/bin2md.rb --index
113-->
114
115
116## filter
117
118| filter command | description |
119| :-------------- | :---------- |
120 | [vcfuniq](./doc/vcfuniq.md) | List unique genotypes. Like GNU uniq, but for VCF records. Remove records which have the same positon, ref, and alt as the previous record. |
121 | [vcfuniqalleles](./doc/vcfuniqalleles.md) | List unique alleles For each record, remove any duplicate alternate alleles that may have resulted from merging separate VCF files. |
122 | [vcffilter](./doc/vcffilter.md) | VCF filter the specified vcf file using the set of filters |
123
124## metrics
125
126| metrics command | description |
127| :-------------- | :---------- |
128 | [vcfcheck](./doc/vcfcheck.md) | Validate integrity and identity of the VCF by verifying that the VCF record's REF matches a given reference file. |
129 | [vcfhethomratio](./doc/vcfhethomratio.md) | Generates the het/hom ratio for each individual in the file |
130 | [vcfhetcount](./doc/vcfhetcount.md) | Calculate the heterozygosity rate: count the number of alternate alleles in heterozygous genotypes in all records in the vcf file |
131 | [vcfdistance](./doc/vcfdistance.md) | Adds a tag to each variant record which indicates the distance to the nearest variant. (defaults to BasesToClosestVariant if no custom tag name is given. |
132 | [vcfentropy](./doc/vcfentropy.md) | Annotate VCF records with the Shannon entropy of flanking sequence. Anotates the output VCF file with, for each record, EntropyLeft, EntropyRight, EntropyCenter, which are the entropies of the sequence of the given window size to the left, right, and center of the record. Also adds EntropyRef and EntropyAlt for each alt. |
133
134## phenotype
135
136| phenotype command | description |
137| :-------------- | :---------- |
138 | [permuteGPAT++](./doc/permuteGPAT++.md) | **permuteGPAT++** is a method for adding empirical p-values to a GPAT++ score. |
139
140## genotype
141
142| genotype command | description |
143| :-------------- | :---------- |
144 | [normalize-iHS](./doc/normalize-iHS.md) | normalizes iHS or XP-EHH scores. |
145 | [hapLrt](./doc/hapLrt.md) | HapLRT is a likelihood ratio test for haplotype lengths. The lengths are modeled with an exponential distribution. The sign denotes if the target has longer haplotypes (1) or the background (-1). |
146 | [abba-baba](./doc/abba-baba.md) | **abba-baba** calculates the tree pattern for four indviduals. This tool assumes reference is ancestral and ignores non **abba-baba** sites. The output is a boolian value: 1 = true , 0 = false for abba and baba. the tree argument should be specified from the most basal taxa to the most derived. |
147
148## transformation
149
150| transformation command | description |
151| :-------------- | :---------- |
152 | [vcfinfo2qual](./doc/vcfinfo2qual.md) | Sets QUAL from info field tag keyed by [key]. The VCF file may be omitted and read from stdin. The average of the field is used if it contains multiple values. |
153 | [vcfsamplediff](./doc/vcfsamplediff.md) | Establish putative somatic variants using reported differences between germline and somatic samples. Tags each record where the listed sample genotypes differ with <tag>. The first sample is assumed to be germline, the second somatic. Each record is tagged with <tag>={germline,somatic,loh} to specify the type of variant given the genotype difference between the two samples. |
154 | [vcfaddinfo](./doc/vcfaddinfo.md) | Adds info fields from the second file which are not present in the first vcf file. |
155 | [vcfremoveaberrantgenotypes](./doc/vcfremoveaberrantgenotypes.md) | strips samples which are homozygous but have observations implying heterozygosity. Remove samples for which the reported genotype (GT) and observation counts disagree (AO, RO). |
156 | [vcfglxgt](./doc/vcfglxgt.md) | Set genotypes using the maximum genotype likelihood for each sample. |
157 | [dumpContigsFromHeader](./doc/dumpContigsFromHeader.md) | Dump contigs from header |
158 | [vcfevenregions](./doc/vcfevenregions.md) | Generates a list of regions, e.g. chr20:10..30 using the variant density information provided in the VCF file to ensure that the regions have even numbers of variants. This can be use to reduce the variance in runtime when dividing variant detection or genotyping by genomic coordinates. |
159 | [vcfcat](./doc/vcfcat.md) | Concatenates VCF files |
160 | [vcfannotategenotypes](./doc/vcfannotategenotypes.md) | Examine genotype correspondence. Annotate genotypes in the first file with genotypes in the second adding the genotype as another flag to each sample filed in the first file. annotation-tag is the name of the sample flag which is added to store the annotation. also adds a 'has_variant' flag for sites where the second file has a variant. |
161 | [vcfafpath](./doc/vcfafpath.md) | Display genotype paths |
162 | [vcfclassify](./doc/vcfclassify.md) | Creates a new VCF where each variant is tagged by allele class: snp, ts/tv, indel, mnp |
163 | [vcfallelicprimitives](./doc/vcfallelicprimitives.md) | If multiple allelic primitives (gaps or mismatches) are specified in a single VCF record, split the record into multiple lines, but drop all INFO fields. Does not handle genotypes (yet). MNPs are split into multiple SNPs unless the -m flag is provided. Records generated by splits have th |
164 | [vcfqual2info](./doc/vcfqual2info.md) | Puts QUAL into an info field tag keyed by [key]. |
165 | [vcfcreatemulti](./doc/vcfcreatemulti.md) | If overlapping alleles are represented across multiple records, merge them into a single record. Currently only for indels. |
166 | [vcfgeno2alleles](./doc/vcfgeno2alleles.md) | modifies the genotypes field to provide the literal alleles rather than indexes |
167 | [vcfsample2info](./doc/vcfsample2info.md) | Take annotations given in the per-sample fields and add the mean, median, min, or max to the site-level INFO. |
168 | [vcfld](./doc/vcfld.md) | Compute LD |
169 | [vcfnumalt](./doc/vcfnumalt.md) | outputs a VCF stream where NUMALT has been generated for each record using sample genotypes |
170 | [vcfstreamsort](./doc/vcfstreamsort.md) | Sorts the input (either stdin or file) using a streaming sort algorithm. Guarantees that the positional order is correct provided out-of-order variants are no more than 100 positions in the VCF file apart. |
171 | [vcfinfosummarize](./doc/vcfinfosummarize.md) | Take annotations given in the per-sample fields and add the mean, median, min, or max to the site-level INFO. |
172 | [vcflength](./doc/vcflength.md) | Add length info field |
173 | [vcfkeepgeno](./doc/vcfkeepgeno.md) | Reduce file size by removing FORMAT fields not listed on the command line from sample specifications in the output |
174 | [vcfcombine](./doc/vcfcombine.md) | Combine VCF files positionally, combining samples when sites and alleles are identical. Any number of VCF files may be combined. The INFO field and other columns are taken from one of the files which are combined when records in multiple files match. Alleles must have identical ordering to be combined into one record. If they do not, multiple records will be emitted. |
175 | [vcfprimers](./doc/vcfprimers.md) | For each VCF record, extract the flanking sequences, and write them to stdout as FASTA records suitable for alignment. |
176 | [vcfflatten](./doc/vcfflatten.md) | Removes multi-allelic sites by picking the most common alternate. Requires allele frequency specification 'AF' and use of 'G' and 'A' to specify the fields which vary according to the Allele or Genotype. VCF file may be specified on the command line or piped as stdin. |
177 | [vcf2dag](./doc/vcf2dag.md) | Modify VCF to be able to build a directed acyclic graph (DAG) |
178 | [vcfcleancomplex](./doc/vcfcleancomplex.md) | Removes reference-matching sequence from complex alleles and adjusts records to reflect positional change. |
179 | [vcfbreakmulti](./doc/vcfbreakmulti.md) | If multiple alleles are specified in a single record, break the record into multiple lines, preserving allele-specific INFO fields. |
180 | [vcfindex](./doc/vcfindex.md) | Adds an index number to the INFO field (id=position) |
181 | [vcfkeepinfo](./doc/vcfkeepinfo.md) | To decrease file size remove INFO fields not listed on the command line |
182 | [vcfgeno2haplo](./doc/vcfgeno2haplo.md) | Convert genotype-based phased alleles within --window-size into haplotype alleles. Will break haplotype construction when encountering non-phased genotypes on input. |
183 | [vcfintersect](./doc/vcfintersect.md) | VCF set analysis |
184 | [vcfannotate](./doc/vcfannotate.md) | Intersect the records in the VCF file with targets provided in a BED file. Intersections are done on the reference sequences in the VCF file. If no VCF filename is specified on the command line (last argument) the VCF read from stdin. |
185 | [smoother](./doc/smoother.md) | smoothes is a method for window smoothing many of the GPAT++ formats. |
186 | [vcf2fasta](./doc/vcf2fasta.md) | Generates sample_seq:N.fa for each sample, reference sequence, and chromosomal copy N in [0,1... ploidy]. Each sequence in the fasta file is named using the same pattern used for the file name, allowing them to be combined. |
187 | [vcfsamplenames](./doc/vcfsamplenames.md) | List sample names |
188 | [vcfleftalign](./doc/vcfleftalign.md) | Left-align indels and complex variants in the input using a pairwise ref/alt alignment followed by a heuristic, iterative left realignment process that shifts indel representations to their absolute leftmost (5') extent. |
189 | [vcfglbound](./doc/vcfglbound.md) | Adjust GLs so that the maximum GL is 0 by dividing all GLs for each sample by the max. |
190 | [vcfcommonsamples](./doc/vcfcommonsamples.md) | Generates each record in the first file, removing samples not present in the second |
191 | [vcfecho](./doc/vcfecho.md) | Echo VCF to stdout (simple demo) |
192 | [vcfkeepsamples](./doc/vcfkeepsamples.md) | outputs each record in the vcf file, removing samples not listed on the command line |
193 | [vcf2tsv](./doc/vcf2tsv.md) | Converts VCF to per-allelle or per-genotype tab-delimited format, using null string to replace empty values in the table. Specifying -g will output one line per sample with genotype information. When there is more than one alt allele there will be multiple rows, one for each allele and, the info will match the 'A' index |
194 | [vcfoverlay](./doc/vcfoverlay.md) | Overlay records in the input vcf files with order as precedence. |
195 | [vcfgenosamplenames](./doc/vcfgenosamplenames.md) | Get samplenames |
196 | [vcfremovesamples](./doc/vcfremovesamples.md) | outputs each record in the vcf file, removing samples listed on the command line |
197 | [vcfremap](./doc/vcfremap.md) | For each alternate allele, attempt to realign against the reference with lowered gap open penalty. If realignment is possible, adjust the cigar and reference/alternate alleles. Observe how different alignment parameters, including context and entropy-dependent ones, influence variant classification and interpretation. |
198 | [vcffixup](./doc/vcffixup.md) | Generates a VCF stream where AC and NS have been generated for each record using sample genotypes |
199
200## statistics
201
202| statistics command | description |
203| :-------------- | :---------- |
204 | [vcfgenosummarize](./doc/vcfgenosummarize.md) | Adds summary statistics to each record summarizing qualities reported in called genotypes. Uses: RO (reference observation count), QR (quality sum reference observations) AO (alternate observation count), QA (quality sum alternate observations) |
205 | [vcfcountalleles](./doc/vcfcountalleles.md) | Count alleles |
206 | [meltEHH](./doc/meltEHH.md) |  |
207 | [genotypeSummary](./doc/genotypeSummary.md) | Generates a table of genotype counts. Summarizes genotype counts for bi-allelic SNVs and indel |
208 | [vcfrandomsample](./doc/vcfrandomsample.md) | Randomly sample sites from an input VCF file, which may be provided as stdin. Scale the sampling probability by the field specified in KEY. This may be used to provide uniform sampling across allele frequencies, for instance. |
209 | [pVst](./doc/pVst.md) | **pVst** calculates vst, a measure of CNV stratification. |
210 | [vcfrandom](./doc/vcfrandom.md) | Generate a random VCF file |
211 | [segmentFst](./doc/segmentFst.md) | **segmentFst** creates genomic segments (bed file) for regions with high wcFst |
212 | [sequenceDiversity](./doc/sequenceDiversity.md) | The **sequenceDiversity** program calculates two popular metrics of haplotype diversity: pi and extended haplotype homozygoisty (eHH). Pi is calculated using the Nei and Li 1979 formulation. eHH a convenient way to think about haplotype diversity. When eHH = 0 all haplotypes in the window are unique and when eHH = 1 all haplotypes in the window are identical. |
213 | [segmentIhs](./doc/segmentIhs.md) | Creates genomic segments (bed file) for regions with high wcFst |
214 | [vcfgenotypes](./doc/vcfgenotypes.md) | Report the genotypes for each sample, for each variant in the VCF. Convert the numerical represenation of genotypes provided by the GT field to a human-readable genotype format. |
215 | [vcfaltcount](./doc/vcfaltcount.md) | count the number of alternate alleles in all records in the vcf file |
216 | [plotHaps](./doc/plotHaps.md) | **plotHaps** provides the formatted output that can be used with 'bin/plotHaplotypes.R'. |
217 | [vcfsitesummarize](./doc/vcfsitesummarize.md) | Summarize by site |
218 | [vcfgenotypecompare](./doc/vcfgenotypecompare.md) | adds statistics to the INFO field of the vcf file describing the amount of discrepancy between the genotypes (GT) in the vcf file and the genotypes reported in the <other-genotype-tag>. use this after vcfannotategenotypes to get correspondence statistics for two vcfs. |
219 | [vcfstats](./doc/vcfstats.md) | Prints statistics about variants in the input VCF file. |
220 | [wcFst](./doc/wcFst.md) | **wcFst** is Weir & Cockerham's Fst for two populations. Negative values are VALID, they are sites which can be treated as zero Fst. For more information see Evolution, Vol. 38 N. 6 Nov 1984. Specifically **wcFst** uses equations 1,2,3,4. |
221 | [permuteSmooth](./doc/permuteSmooth.md) | **permuteSmooth** is a method for adding empirical p-values smoothed wcFst scores. |
222 | [bFst](./doc/bFst.md) | **bFst** is a Bayesian approach to Fst. Importantly **bFst** accounts for genotype uncertainty in the model using genotype likelihoods. For a more detailed description see: `A Bayesian approach to inferring population structure from dominant markers' by Holsinger et al. Molecular Ecology Vol 11, issue 7 2002. The likelihood function has been modified to use genotype likelihoods provided by variant callers. There are five free parameters estimated in the model: each subpopulation's allele frequency and Fis (fixation index, within each subpopulation), a free parameter for the total population's allele frequency, and Fst. |
223 | [vcfroc](./doc/vcfroc.md) | Generates a pseudo-ROC curve using sensitivity and specificity estimated against a putative truth set. Thresholding is provided by successive QUAL cutoffs. |
224 | [vcfparsealts](./doc/vcfparsealts.md) | Alternate allele parsing method. This method uses pairwise alignment of REF and ALTs to determine component allelic primitives for each alternate allele. |
225 | [pFst](./doc/pFst.md) | **pFst** is a probabilistic approach for detecting differences in allele frequencies between two populations. |
226 | [iHS](./doc/iHS.md) | **iHS** calculates the integrated haplotype score which measures the relative decay of extended haplotype homozygosity (EHH) for the reference and alternative alleles at a site (see: voight et al. 2006, Spiech & Hernandez 2014). |
227 | [popStats](./doc/popStats.md) | General population genetic statistics for each SNP |
228
229See also [vcflib.md](./doc/vcflib.md).
230
231## scripts
232
233The vcflib source repository contains a number of additional scripts.
234Click on the link to see the source code.
235
236| script | description |
237| :-------------- | :---------- |
238| [vcfclearinfo](./scripts/vcfclearinfo) | clear INFO field |
239| [vcfqualfilter](./scripts/vcfqualfilter) | quality filter |
240| [vcfnulldotslashdot](./scripts/vcfnulldotslashdot) | rewrite null genotypes to ./. |
241| [vcfprintaltdiscrepancy.r](./scripts/vcfprintaltdiscrepancy.r) | show ALT discrepancies in a table |
242| [vcfremovenonATGC](./scripts/vcfremovenonATGC) | remove non-nucleotides in REF or ALT |
243| [plotSmoothed.R](./scripts/plotSmoothed.R) | smooth plot of wcFst, pFst or abba-baba |
244| [vcf_strip_extra_headers](./scripts/vcf_strip_extra_headers) | strip headers |
245| [plotHapLrt.R](./scripts/plotHapLrt.R) | plot results of pFst |
246| [vcfbiallelic](./scripts/vcfbiallelic) | remove anything that is not biallelic |
247| [vcfsort](./scripts/vcfsort) | sort VCF using shell script |
248| [vcfnosnps](./scripts/vcfnosnps) | remove SNPs |
249| [vcfmultiwayscripts](./scripts/vcfmultiwayscripts) | more multiway comparisons |
250| [vcfgtcompare.sh](./scripts/vcfgtcompare.sh) | annotates records in the first file with genotypes and sites from the second |
251| [plotPfst.R](./scripts/plotPfst.R) | plot pFst |
252| [vcfregionreduce_and_cut](./scripts/vcfregionreduce_and_cut) | reduce, gzip, and tabix |
253| [plotBfst.R](./scripts/plotBfst.R) | plot results of pFst |
254| [vcfnobiallelicsnps](./scripts/vcfnobiallelicsnps) | remove biallelic SNPs |
255| [vcfindels](./scripts/vcfindels) | show INDELS |
256| [vcfmultiway](./scripts/vcfmultiway) | multiway comparison |
257| [vcfregionreduce](./scripts/vcfregionreduce) | reduce VCFs using a BED File, gzip them up and create tabix index |
258| [vcfprintaltdiscrepancy.sh](./scripts/vcfprintaltdiscrepancy.sh) | runner |
259| [vcfclearid](./scripts/vcfclearid) | clear ID field |
260| [vcfcomplex](./scripts/vcfcomplex) | remove all SNPs but keep SVs |
261| [vcffirstheader](./scripts/vcffirstheader) | show first header |
262| [plotXPEHH.R](./scripts/plotXPEHH.R) | plot XPEHH |
263| [vcfregionreduce_pipe](./scripts/vcfregionreduce_pipe) | reduce, gzip and tabix in a pipe |
264| [vcfplotaltdiscrepancy.sh](./scripts/vcfplotaltdiscrepancy.sh) | plot ALT discrepancy runner |
265| [vcfplottstv.sh](./scripts/vcfplottstv.sh) | runner |
266| [vcfnoindels](./scripts/vcfnoindels) | remove INDELs |
267| [bgziptabix](./scripts/bgziptabix) | runs bgzip on the input and tabix indexes the result |
268| [plotHaplotypes.R](./scripts/plotHaplotypes.R) | plot results |
269| [vcfplotsitediscrepancy.r](./scripts/vcfplotsitediscrepancy.r) | plot site discrepancy |
270| [vcfindelproximity](./scripts/vcfindelproximity) | show SNPs around an INDEL |
271| [bed2region](./scripts/bed2region) | convert VCF CHROM column in VCF file to region |
272| [vcfplotaltdiscrepancy.r](./scripts/vcfplotaltdiscrepancy.r) | plot ALT discrepancies |
273| [plot_roc.r](./scripts/plot_roc.r) | plot ROC |
274| [vcfmultiallelic](./scripts/vcfmultiallelic) | remove anything that is not multiallelic |
275| [vcfsnps](./scripts/vcfsnps) | show SNPs |
276| [vcfvarstats](./scripts/vcfvarstats) | use fastahack to get stats |
277| [vcfregionreduce_uncompressed](./scripts/vcfregionreduce_uncompressed) | reduce, gzip and tabix |
278| [plotWCfst.R](./scripts/plotWCfst.R) | plot wcFst |
279| [vcf2bed.py](./scripts/vcf2bed.py) | transform VCF to BED file |
280| [vcfjoincalls](./scripts/vcfjoincalls) | overlay files using QUAL and GT from a second VCF |
281| [vcf2sqlite.py](./scripts/vcf2sqlite.py) | push VCF file into SQLite3 database using dbname |
282
283# Development
284
285## build from source
286
287VCFLIB uses the cmake build system, after a recursive checkout
288of the sources make the files in the ./build directory with:
289
290```sh
291git clone --recursive https://github.com/vcflib/vcflib.git
292cd vcflib
293mkdir -p build && cd build
294cmake ..
295cmake --build .
296cmake --install .
297```
298
299and to run the tests
300
301```sh
302ctest --verbose
303```
304
305Executables are built into the `./build` directory in the repository.
306
307Build dependencies can be viewed in the Travis-CI and github-CI
308scripts (see badges above), as well as [guix.scm](./guix.scm) used by
309us to create the build environment. Essentially:
310
311- C++ compiler
312- htslib
313- tabixpp
314
315For include files add
316
317- libhts-dev
318- libtabixpp-dev
319- libtabixpp0
320
321And for some of the VCF executables
322
323- python
324- perl
325
326## link library
327
328The standard build creates `build/vcflib.a`. Take a hint from the
329[cmake](./CMakeLists.txt) file that builds all the vcflib tools.
330
331## source code
332
333See [vcfecho.cpp](./src/vcfecho.cpp) for basic usage.
334[Variant.h](./src/Variant.h) and [Variant.cpp](./src/Variant.cpp)
335describe methods available in the API.  vcflib is incorporated into
336several projects, such as
337[freebayes](https://github.com/freebayes/freebayes), which may provide
338a point of reference for prospective developers.  Note vcflib contains
339submodules (git repositories) comprising some dependencies. A full
340Guix development environment we use is defined [here](./guix.scm).
341
342
343# LICENSE
344
345This software is distributed under the free software [MIT
346LICENSE](./LICENSE).
347