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
|
#include <seqan/vcf_io.h>
using namespace seqan2;
int main()
{
// Open input file.
VcfFileIn vcfIn(toCString(getAbsolutePath("demos/tutorial/vcf_io/example.vcf")));
// Attach to standard output.
VcfFileOut vcfOut(vcfIn);
open(vcfOut, std::cout, Vcf());
// Copy over header.
VcfHeader header;
readHeader(header, vcfIn);
writeHeader(vcfOut, header);
// Copy the file record by record.
VcfRecord record;
while (!atEnd(vcfIn))
{
readRecord(record, vcfIn);
writeRecord(vcfOut, record);
}
return 0;
}
|
