// Calculate r2, lf, hd double nhd[] = {0,0,0,0,0,0,0,0,0}; double ab = 0, aa = 0, bb = 0, a = 0, b = 0; int nab = 0, ndiff = 0; int an_tot = 0, bn_tot = 0; for (i=0; i<arec->n_sample; i++)...
Generalization and better support for custom functions that allow adding new INFO tags based on arbitrary-i, --includetype of expressions. For example, to calculate a missing INFO/DP annotation from FORMAT/AD, it is possible to use: -t 'DP:1=int(sum(FORMAT/AD))' ...
(see *<<expressions,EXPRESSIONS>>*) %POS0 POS in 0-based coordinates %PBINOM(TAG) Calculate phred-scaled binomial probability, the allele index is determined from GT %SAMPLE Sample name %TAG{INT} Curly brackets to print a subfield (e.g. INFO/TAG{1}, the indexes are 0-based) %TBCSQ ...