varianttools icon indicating copy to clipboard operation
varianttools copied to clipboard

--from_stat VTOOLS v 2.7.0

Open BoPeng opened this issue 8 years ago • 1 comments

I have 49 individuals in one sample but actually, using --from_stat to define the number of WT genotype for example I obtain numbers like 98.. My sample is well define.

When I use:

vtools update variant --from_stat 'N_heterozygalt=#(het)' --samples "pheno='0'" -j8

Vtools INFO is: INFO: 49 samples are selected which is the correct number of individuals but, when I do an output, I have:

chr   ref   alt   #(wtGT)  #(het)  #(hom) #(alt)
1        T       A             4       24      70      164

but I have only 49 individuals..

When I used vtools v2.6.2 I did not have this issue.

BoPeng avatar Jan 05 '17 15:01 BoPeng

Could you send me the exact command for your do an output? It would be best if you could send me, even in private email, a sample vcf file and commands to reproduce the problem. Thanks.

BoPeng avatar Jan 05 '17 15:01 BoPeng