plink, vcftool计算等位基因频率(allele frequency,vcf)

摘要:
计算等位基因频率有两种方法。第一种方法是使用vcftool:/path/to/vcftools--vcffile Vcf--freq--chr1--outfilefreq一个简单的命令行,文件Vcf引用要输入的Vcf文件。--freq表示计算等位基因频率。——chr之后的1表示要计算的区域在1号染色体上。当然,您也可以选择要计算的染色体区域。Filefreq指的是要输出的文件名。

计算等位基因频率有两种方式,第一种用vcftool计算:

/path/to/vcftools --vcf file.vcf --freq --chr 1 --out filefreq

  很简单的一个命令行,file.vcf指的是你要输入的vcf文件,--freq表示计算等位基因频率,--chr后面的1表示你要计算的区域在1号染色体,当然,你也可以选择你想计算的染色体区域,filefreq指的是输出的文件名。

结果如下图所示:

plink, vcftool计算等位基因频率(allele frequency,vcf)第1张

第二种用plink计算:

/path/to/plink-1.07-x86_64/plink --noweb --bfile file --freq --out filefrquency

  效果如下图所示:

plink, vcftool计算等位基因频率(allele frequency,vcf)第2张

.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

免责声明:文章转载自《plink, vcftool计算等位基因频率(allele frequency,vcf)》仅用于学习参考。如对内容有疑问,请及时联系本站处理。

上篇OpenShift实战(一):OpenShift安装delphi 文件的操作下篇

宿迁高防,2C2G15M,22元/月;香港BGP,2C5G5M,25元/月 雨云优惠码:MjYwNzM=

相关文章

全基因组关联分析(Genome-Wide Association Study,GWAS)流程

全基因组关联分析流程: 一、准备plink文件 1、准备PED文件 PED文件至少有六列,内容如下: Family ID Individual ID Paternal ID Maternal ID Sex (1=male; 2=female; other=unknown) Phenotype(-9 missing 0 missing 1 unaffecte...

GWAS学习笔记(一) | 质量控制(QC)

本系列文章采用的数据集与代码来自https://github.com/MareesAT/GWA_tutorial。 该教程获得了许多人的推荐,是一份很详细的step-by-step guide。 本文将介绍该教程中的QC部分(1_QC_GWAS.zip),后续或将继续添加有关QC的其他细节。 0. 目录 目录 0. 目录 1. 准备 2. 数据简介...