1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
|
% -*- coding: utf-8 -*-
# VCFLib API
This document describes the VCFLIB API as it is used by the vcflib modules and the [python ffi](./pyvcflib.md).
vcflib follows the [VCF standard](http://samtools.github.io/hts-specs/VCFv4.1.pdf).
VCFLIB contains a lot of functionality, but the basis of going through a VCF file and fetching record (by record) information is straightforward and visible in all modules. A recent example can be found in [vcfwave](../../src/vcfwave.cpp).
## Open the VCF file
```C++
VariantCallFile variantFile;
if (optind < argc) {
string filename = argv[optind];
variantFile.open(filename);
} else {
variantFile.open(std::cin);
}
if (!variantFile.is_open()) {
return 1;
}
```
## Read records
The following will parse the records and you can print out the first two fields with
```C++
Variant var(variantFile);
while (variantFile.getNextVariant(var)) {
cout << var.sequenceName << " " << var.position << endl;
}
```
## Other fields
In the file [Variant.h](https://github.com/vcflib/vcflib/blob/master/src/Variant.h) the Variant class is defined with fields/accessors, such as
```C++
string sequenceName;
long position;
long zeroBasedPosition(void) const;
string id;
string ref;
vector<string> alt; // a list of all the alternate alleles present at this locus
vector<string> alleles; // a list all alleles (ref + alt) at this locus
```
See above read records example to parse a file. Some things to know are that info fields are split into
fields with values in 'info' and flags that are true in 'infoFlags'.
The order of info fields is not kept in the C++ map data structure so we have to keep track of order in an infoKeys vector.
## Output a VCF record
The default string outputter of the Variant class outputs a VCF record using the field that are defined:
```C++
Variant var(variantFile);
while (variantFile.getNextVariant(var)) {
cout << var << endl;
}
```
will output VCF.
## Mirroring in the Python FFI
The Python FFI follows this API though some accessors may be renamed. See
[pythonffi.cpp](pythonffi.cpp) and [pyvcflib.md](pyvcflib.md).
|
