使用 PLINK 把 vcf 的0/0,0/1,1/1转为字母格式的基因型(比如AA,AG,GG)

file.vcf文件如下所示，包含两个样本、四个变异位点：

#CHROM  POS    ID    REF     ALT     QUAL    FILTER  INFO    FORMAT  ID1     ID2
chr1    123    rs1   G       A       .       PASS    .       GT      0|0     1|0
chr1    124    rs2   C       A       .       PASS    .       GT      1|1     1|1
chr1    125    rs3   G       T       .       PASS    .       GT      0|0     0|0
chr1    126    rs4   G       A       .       PASS    .       GT      1|1     1|1

现在我想把数字基因型变成字母基因型，比如对于rs1 ，我希望0|0变成GG，1|0变成AG；

对于此需求，要用到compound-genotypes参数；

废话不多说，直接给命令：

plink --vcf file.vcf --recode compound-genotypes --out recode

注意，这里的plink是 1.9 版本哦（https://s3.amazonaws.com/plink1-assets/plink_linux_x86_64_20210606.zip）；

运行完以上命令后，会得到 recode.ped 和 recode.map 两个文件；

recode.ped文件如下所示：

ID1 ID1 0 0 0 -9 GG AA GG AA
ID2 ID2 0 0 0 -9 AG AA GG AA

recode.map文件如下所示：

chr1       rs1       0       123
chr1       rs2       0       124
chr1       rs3       0       125
chr1       rs4       0       126

督促我写下这篇教程的原因是，我今天搜遍了我的博客（https://www.cnblogs.com/chenwenyan/）、用了很多关键词，都找不到相应的参数，但是我清楚记得PLINK是有这个功能的。如果我以前有记录，今天就不需要大费周章寻找了，可见记录是一件多么重要的事；

我知道很多人看不上这种没有技术含量的教程，但这些小技巧有时候能帮人省下很多不必要的coding时间；

致谢橙子牛奶糖（陈文燕），请用参考模版：We thank the blogger (orange_milk_sugar, Wenyan Chen) for XXX

感谢小可爱们多年来的陪伴，我与你们一起成长~

使用 PLINK 把 vcf 的0/0,0/1,1/1转为字母格式的基因型(比如AA,AG,GG)相关推荐

libgstreamer-1.0.so.0: cannot open shared object file: No such file or directory
1. 问题现象 error while loading shared libraries: libgstreamer-1.0.so.0: cannot open shared object file: ...
c+语言+null,C/C++语言中NULL、'\0’和0的区别
NULL.'\0'和0的值是一样的,都是0,不过它们的表现形式不一样: 1. NULL: 即空指针,不过在C和C++中并不一样.在VS 2013的库文件string.h中可以看到如果定义. 1 /* ...
Ubuntu14.04 64位机上安装cuda8.0+cudnn5.0操作步骤
查看Ubuntu14.04 64位上显卡信息,执行: lspci | grep -i vga lspci -v -s 01:00.0 nvidia-smi 第一条此命令可以显示一些显卡的相关信息:如果 ...
Spring Cloud Alibaba 基础教程：Nacos 生产级版本 0.8.0
Spring Cloud Alibaba 基础教程:Nacos 生产级版本 0.8.0 昨晚Nacos社区发布了第一个生产级版本:0.8.0.由于该版本除了Bug修复之外,还提供了几个生产管理非常重要 ...
Silverlight 3发布新版3.0.50106.0
微软1月19日发布Silverlight 3新版本3.0.50106.0. 该版本主要修复以下几个问题: 问题一: 当使用图形硬件加速功能(GPU)的时候,如果GPU驱动报错,Silverlight ...
AS1.0(2.0)中的XML示例
虽然Flash早就升级为AS3.0,但是FMS的服务端编程依然仅支持AS1.0(2.0),服务端与.net通讯的最简单方式莫过于请求一个RESTful的webService或wcf,通过它们返回的xm ...
多数编程语言里的0.1+0.2≠0.3?
作者 | Parul Malhotra 译者 | Raku 出品 | AI科技大本营(ID:rgznai100) 我们从小就被教导说0.1+0.2=0.3,但是在奇妙的计算机编程世界里面,事情变得不一 ...
flannel原理初探针对0.1.0版本
flannel flannel是针对k8s设计的三层的网络解决方案.在k8s中为了使pod之间能够使用一种偏平的网络架构,从而完成跨Pod的网络通信. 官网给的原理图如下: flannel 使用TUN ...
GTX 1080Ti + cuda8.0 + cuDNN6.0 安装及测试
GPU 显卡厂商已经安装好了,直接安装 cuda8.0 + cuDNN6.0 我这里的显卡是 GTX 1080 Ti cuda安装我下载的是cuda8.0的是deb格式的1.9个G地址:https: ...

使用 PLINK 把 vcf 的0/0,0/1,1/1转为字母格式的基因型(比如AA,AG,GG)

使用 PLINK 把 vcf 的0/0,0/1,1/1转为字母格式的基因型(比如AA,AG,GG)相关推荐

最新文章

热门文章