Require Genotypes In Vcf File In Order To Output Frequency Statistics
1
0
Entering edit mode
22 months ago
ja4123 • 0

I have a merged .vcf file, that although wasn't merged using vcftools seems almost identical to vcf-merge files. If I try and use --freq to calculate allele frequencies from this file I get the error message "Require Genotypes in VCF file in order to output Frequency Statistics". This seems odd as the file has the same GT:GQ:PL format as other .vcf files.

I was wondering if anyone knows what vcftools uses to calculate --freq and why I am getting this error message. I have gone back to source files, merged these in vcftools then performed --freq with this file, this time the analysis "works" but I get incorrect allele frequencies (compared to working them out by hand).

Regards

vcf vcftools AF allele frequency • 983 views
ADD COMMENT
0
Entering edit mode

Can you grep the CHROM... header line and the first couple lines of the VCF file and paste it here? IE the output from: grep -m 1 "^#CHROM" file.vcf and grep -m 5 -v "^#" file.vcf

ADD REPLY
0
Entering edit mode
22 months ago

Is the header information for the VCF file correct? If the GT entry is missing in the header then vcftools may assume it's also absent from the body of the file.

ADD COMMENT

Login before adding your answer.

Traffic: 1082 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6