VCF NC Format

[aew]

Hi,

I ran VARSCAN2 mpileup2snp to call SNPs and generated a vcf file.

However, I noticed a strange thing about the NC flag. This is apparently INFO=<id=nc,number=1,type=integer,description="number of="" samples="" not="" called"="">

However, I see

TR43100|c0_g1_i1    792 .   T   C   .   PASS    ADP=0;WT=1;HET=0;HOM=1;NC=17    GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR    ./.:.:0 ./.:.:0 ./.:.:0 ./.:.:1 ./.:.:0 ./.:.:1 ./.:.:0 ./.:.:0 ./.:.:1 0/0:10:5:5:5:0:0%:1E0:38:0:0:5:0:0  ./.:.:0 ./.:.:1 ./.:.:1 1/1:7:2:2:0:2:100%:1.6667E-1:0:33:0:0:1:1   ./:.:0  ./:.:0  ./:.:0  ./:.:0  ./:.:0

This has NC=17 and I can see from the 19 samples 2 have been called. This seems correct.

TR43100|c1_g1_i1    351 .   T   C   .   PASS    ADP=0;WT=1;HET=1;HOM=0;NC=17    GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR    ./:.:0  ./:.:0  ./:.:0  ./:.:0  ./:.:0  ./:.:1  ./:.:0  ./:.:0  ./:.:0  ./:.:2  ./:.:0  ./:.:0  ./:.:0  0/1:6:4:4:2:2:50%:2.1429E-1:40:38:2:0:1:1   ./:.:0  ./:.:0  ./:.:0  ./:.:0  ./:.:0

This has NC=17 abut I can only see one sample that has been called? What have I missed, why isn't NC=18 here?

And it says WT=1. Where is the WT? This should be 0/0

And also

TR43100|c0_g1_i1 763 . G A . PASS ADP=0;WT=2;HET=0;HOM=1;NC=16 GT:GQ:SDP:DP:RD:AD:FREQ:PVAL:RBQ:ABQ:RDF:RDR:ADF:ADR ./.:.:0 ./.:.:0 ./.:.:0 1/1:13:3:3:0:3:100%:5E-2:0:39:0:0:0:3 ./.:.:0 ./.:.:1 ./.:.:0 ./.:.:0 ./.:.:1 0/0:6:3:3:3:0:0%:1E0:40:0:0:3:0:0 ./.:.:0 ./.:.:1 ./.:.:1 ./.:.:2 ./.:.:0 ./.:.:0 ./.:.:0 ./.:.:0 ./.:.:0

Where are the 2 WT? I can only see 1 0/0

[GenoMax]

For reference/cross-posted: https://www.biostars.org/p/180308/