Is There A Tool That Would Directly Encode Minor Alllele Freq In The Info Column Of A Vcf?
0
0
Entering edit mode
10.8 years ago
rickyflintoff ▴ 100

Is there a tool that would directly encode minor alllele freq in the INFO column of a VCF? Or write the minor allele frequency to a given file with chrom pos pair?

vcf • 1.9k views
ADD COMMENT
0
Entering edit mode

you wanted to calculate the frequency against 1000genomes/hapmap/dbsnp or ...?
you can write a script to calculate the frequency by using the info, if available, in your vcf file; by dividing AN (allele number) by AC (allele count)

ADD REPLY
0
Entering edit mode

I think the OP wants a tool that will pull the MAFs from dbSNP for each rs, which on first pass might be useful for filtering. If so, use this with caution as many disease-causing alleles are nevertheless in dbSNP. AFAIK, neither of the 2 tools I am most familiar with for vcf annotation (annovar and SeattleSeq) proviide MAF values, and the reason is probably that information is not as useful as it appears to be at first. See http://www.openbioinformatics.org/annovar/annovar_filter.html#dbsnp for a discussion.

ADD REPLY

Login before adding your answer.

Traffic: 1710 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