首页 > 其他分享 >gatk 实现对vcf文件的合并

gatk 实现对vcf文件的合并

时间:2022-12-01 19:33:51浏览次数:46  
标签:rw vcf 合并 seg1 --. gatk Dec root

 

001、测试数据

[root@PC1 test]# ls                            ## 测试数据
seg1_1.vcf  seg1_2.vcf  seg1_3.vcf
[root@PC1 test]# ll -h
total 1.2G
-rw-r--r--. 1 root root 392M Dec  1 18:14 seg1_1.vcf
-rw-r--r--. 1 root root 410M Dec  1 18:15 seg1_2.vcf
-rw-r--r--. 1 root root 386M Dec  1 18:14 seg1_3.vcf

 

 

002、对vcf文件进行合并

gatk MergeVcfs -I seg1_1.vcf -I seg1_2.vcf -I seg1_3.vcf -O xxx.vcf      ## 合并命令

 

合并结果如下:

[root@PC1 test]# ls
seg1_1.vcf  seg1_2.vcf  seg1_3.vcf  xxx.vcf  xxx.vcf.idx
[root@PC1 test]# ll -h
total 2.4G
-rw-r--r--. 1 root root 392M Dec  1 18:14 seg1_1.vcf
-rw-r--r--. 1 root root 410M Dec  1 18:15 seg1_2.vcf
-rw-r--r--. 1 root root 386M Dec  1 18:14 seg1_3.vcf
-rw-r--r--. 1 root root 1.2G Dec  1 19:16 xxx.vcf
-rw-r--r--. 1 root root 6.2K Dec  1 19:16 xxx.vcf.idx

 

 

003、将需要合并的vcf文件生成列表

[root@PC1 test]# ls
seg1_1.vcf  seg1_2.vcf  seg1_3.vcf  vcf.list
[root@PC1 test]# cat vcf.list            ## 待合并的vcf列表文件
seg1_1.vcf
seg1_2.vcf
seg1_3.vcf
[root@PC1 test]# java -jar /home/software/picard/picard.jar MergeVcfs I=vcf.list O=xxx.vcf       ## 合并命令

 

 

004、合并结果

[root@PC1 test]# ls
seg1_1.vcf  seg1_2.vcf  seg1_3.vcf  vcf.list  xxx.vcf  xxx.vcf.idx
[root@PC1 test]# ll -h
total 2.4G
-rw-r--r--. 1 root root 392M Dec  1 18:14 seg1_1.vcf
-rw-r--r--. 1 root root 410M Dec  1 18:15 seg1_2.vcf
-rw-r--r--. 1 root root 386M Dec  1 18:14 seg1_3.vcf
-rw-r--r--. 1 root root   33 Dec  1 19:24 vcf.list
-rw-r--r--. 1 root root 1.2G Dec  1 19:24 xxx.vcf
-rw-r--r--. 1 root root 6.2K Dec  1 19:24 xxx.vcf.idx

 

标签:rw,vcf,合并,seg1,--.,gatk,Dec,root
From: https://www.cnblogs.com/liujiaxin2018/p/16942450.html

相关文章

  • 北理工45. 【字符】合并字符串
    45.【字符】合并字符串 输入两个已经按从小到大顺序排列好的字符串,编写一个合并两个字符串的函数,使合并后的字符串,仍然是从小到大排列。输入:          ......
  • 2、通知签名(把两个文件合并成一个文件)
    如果通知有两份文件的话,那么领导若不在办公室,在家办公,那么他签名的时候只会签第二页的文件,他反馈回来的是一张图片,那么如何把一个【pdf文件】和【图片】合并成一张图片用......
  • 区间合并
    给定 nn 个区间 [l,r],要求合并所有有交集的区间。注意如果在端点处相交,也算有交集。输出合并完成后的区间个数。#include<iostream>#include<algorithm>#includ......
  • C# List合并相同元素
     C#List合并相同元素,对相同属性的列进行求和原数据合并后     publicclassStudent{publicstring?Name{get;set;}pu......
  • 合并两个有序链表
    21.合并两个有序链表根据题目要求可以知道这两个链表是升序的,需要将这两个链表合成一个新的升序链表并返回方法一:迭代当list1和list2都不是空链表时,判断list1和l......
  • layui合并单元格
    在别人的基础上解决了多列合并和同一个页面多个表格的问题1//合并单元格2functionmerge(id,res,columsName,columsIndex){34var......
  • git的回退以及合并,删除什么的
    有时候不小心合并了别的分支中的commit.我们需要回退某些提交记录。可以通过reset来操作,reset会回退到指定commit。这种方式会删除记录,我们最好使用revert命令来操作git......
  • bat脚本查找指定目录下的文件并将所有文件内容合并输出到一个文件里 写软著专用提取代
    1.将下面的代码复制到 并另存为【查找文件并合并内容.bat】里echooff&color0A::判断输入路径是不是文件夹,如果是,则获取文件夹下的所有文件名(包括子文件夹下的):......
  • 【算法训练营day20】LeetCode654. 最大二叉树 LeetCode617. 合并二叉树 LeetCode700.
    LeetCode654.最大二叉树题目链接:654.最大二叉树初次尝试和昨天最后一题的思路很像,本质上都是递归构建二叉树。classSolution{public:TreeNode*constructMa......
  • 合并两个不同远程仓库的Git命令
    需求场景描述同一个项目的代码,在coding上有一份,在GitHub上也有一份。但coding上的代码较新,领先了几个commit,GitHub上的比较旧。这是想要将GitHub上的代码与coding上的代码......