Question: Vcftools --Freq
1
gravatar for kai.clarkson
6.4 years ago by
kai.clarkson20
kai.clarkson20 wrote:

I'm getting some really confusing results using vcftools and I was wondering if anyone could shed some light. I have merged 5 individual SNP call .vcf files using vcf-merge. This merged file seems fine until I try to calculate allele frequency from it (--freq). Some of the allele freq are correct (compared to calculating by hand) but some are wrong. The Fst calculator (--fst-pop) also throws up some weird highly negative scores.

The genotype scores for the individuals all seem to be correct (the same as the source files) in the rows that give wrong allele frequencies, so how is vcftools getting an incorrect value?

Any help would be greatly appreciated,

Regards,

Kai

vcf frequency allele vcftools • 3.3k views
ADD COMMENTlink written 6.4 years ago by kai.clarkson20

I have never worked with calculating AF from a merged vcf file using --freq option.What I have noticed is the merged file you get from vcftools usually takes AF from the first vcf input.It doesn't give you an average of the AF over all the samples.I am not sure if this would help,but just an idea.

ADD REPLYlink modified 6.4 years ago • written 6.4 years ago by Kssr110

Could you post few lines of the input and output?

ADD REPLYlink written 6.4 years ago by sa9800
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1024 users visited in the last hour