計算等位基因頻率有兩種方式,第一種用vcftool計算:
/path/to/vcftools --vcf file.vcf --freq --chr 1 --out filefreq
很簡單的一個命令行,file.vcf指的是你要輸入的vcf文件,--freq表示計算等位基因頻率,--chr后面的1表示你要計算的區域在1號染色體,當然,你也可以選擇你想計算的染色體區域,filefreq指的是輸出的文件名。
結果如下圖所示:

第二種用plink計算:
/path/to/plink-1.07-x86_64/plink --noweb --bfile file --freq --out filefrquency
效果如下圖所示:

.frq (basic allele frequency report)
Produced by --freq. Valid input for --read-freq.
A text file with a header line, and then one line per variant with the following six fields:
| CHR |
Chromosome code |
| SNP |
Variant identifier |
| A1 |
Allele 1 (usually minor) |
| A2 |
Allele 2 (usually major) |
| MAF |
Allele 1 frequency |
| NCHROBS |
Number of allele observations |
