VCF header causes NR mismatch

[janxkoci]

I noticed that NR is off by the size of a VCF header, as if NR was set before header is removed (the header normally doesn't print with -c vcf flag).

This means if my VCF has header of 10 lines, the first variant is at NR==11, which is not expected.

[janxkoci]

Another way to demonstrate the behaviour - this code will print two different numbers:

$ bioawk -c vcf '{c++} END{print c, NR}' dogs.LDpruned.vcf.gz 
17651	17657

This is unexpected, because if I just print the file I get this:

$ bioawk -c vcf 1 dogs.LDpruned.vcf.gz | wc -l
17651

Intuitively, I'd expect that number to be the same as NR, but it's not.