gatk gvcf合并之CombineGVCF

于 2024-03-03 14:46:09 发布
阅读量724 收藏 9
点赞数 5
文章标签： linux python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/jinning3825/article/details/136431651
版权

啊啊啊啊啊啊啊

我发现CombineGVCF很慢，每次都加几个样本，我都是从头再合并，忽然想到为啥我不用合并后的60个样本，再去合并新添的10个呢！！！

肯定省钱省力省时间！！！

金金金鱼
关注
  • 5
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
2020.11.20丨使用GATK CombineGVCFs命令批量合并vcf文件
11-20 1万+
GATK是一款强大的数据处理软件，最近在优化GWAS流程时遇到一个麻烦事，就是要将各样品的VCF文件进行合并，本来GATK里面有一个可以合并VCF数据的命令 CombineGVCFs，可以将所有样品的VCF合并成一个文件。但是这个命令需要一个一个输入文件名。 熟悉GWAS的小伙伴应该清楚，GWAS项目动辄上百个样品，让人一个一个输入还是很繁琐的。因此我写了个shell脚本，能够快速输入样品名称，并执行CombineGVCFs命令。 脚本 Ref_genome="genome.fna"
GATK4注意事项
weixin_30586257的博客
08-01 2567
近期在测试多样品的WES的过程中发现用HC得到gvcf之后，合并多个样品的gvcf文件的过程中，使用CombineGVCFs的过程中很慢，发现官网推荐使用GenomicsDBImport 用法如下： gatk GenomicsDBImport \ -V data/gvcfs/mother.g.vcf \ -V data/gvcfs/father.g.vcf \...
GATK4 流程分析- 从fastq到vcf（样例数据处理）
zhanghongyi_cpp的博客
10-03 1653
参考基因组采用在ensembl下载的鸡的参考基因组（以下命名为1.fa）以下为对参考基因组的预处理，三个步骤一个也不能少大于2G的参考基因组要用参数-a bwtsw生成1.fa.fai文件生成参考基因组的dict⽂件。
java queue GATK_gatk4使用总结
weixin_33673142的博客
03-06 1069
昨天看了gatk的官网，从2018年发布正式版的4.0.0开始，到现在已经更新到4.1.8，在速度和准确度上都有了大幅的提升。gatk4除了整合picard软件之外，在使用上与gatk3基本相同，只不过是在命令运行、功能划分及运行速度上进行了调整。gatk4软件的安装及使用，网上有很多帖子，可以自行查阅。今天主要给大家介绍gatk4与gatk3明显不同的地方，这些使得gatk4的应用更加方便。1....
多样本GVCF文件合并
weixin_59587425的博客
09-07 1882
在千人基因组下载了几百个中国人的GVCF文件，每个样本3个G，用CombineGVCFs合并会报错，后来用bcftool挺好的。个警告，不知道为什么。
多个VCF文件合并处理方法
06-01
多个VCF文件合并成一个VCF文件的处理方法,此方法必定可以解决您的烦恼
bcftools合并vcf文件
weixin_33757609的博客
07-06 1226
见命令： bcftools merge A.vcf.gz B.vcf.gz C.vcf.gz -Oz -o ABC.vcf.gz 　　参考链接：http://vcftools.sourceforge.net/htslib.html#merge 转载于:https://www.cnblogs.com/chenwenyan/p/9273726.html...
GATK4 实用技巧丨线程数和内存大小怎么设置？HaplotypeCaller、GenotypeGVCFs、CombineGVCFs、MarkDuplicates
青笋的博客
08-08 3233
从后台数据来看，阅读量的来源大部分是通过推荐获得，说明现在公众号的推送机制是与内容质量挂钩，好的内容系统会自动推荐给更多适合的人，形成良性循环机制，也会对速度造成影响，当超过20个线程后，随着线程增多速度不会有明显变化，因此在实际运算过程中根据服务器配置，选择20个线程以下，这样效率最高。，虽然文章的内容很简单，在生信圈子里属于小菜一碟，但还是获得了70多位朋友的转发，而且大部分朋友都看完了整篇内容。这是一个分享学习笔记的公众号，创建不到一年，根本算不上什么，能有人愿意来读，已经是莫大的鼓励。
实操 | 合并VCF文件的几种方法及注意事项
热门推荐
悟道西方
08-15 1万+
背 景 在基因组分析领域的很多不同场景中，需要合并VCF文件。VCF文件。简单来说，就是记录样本基因型的文件。但多数VCF文件不只记录了基因型，也包含有关该基因型的来源的细节。其它文件。VCF文件的上游是BAM文件，主要记录Reads与参考基因组的比对信息；更上游的，就是FASTQ测序数据，以及物种的参考基因组。不同类型的VCF文件。VCF文件有单样本的、多样本的；也有普通VCF文件...
GATK使用方法详细介绍
07-31
GATK软件使用方法详细描述，好的资源不容错过！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！
gvcf-call:从gVCF文件调用变体
04-09
gvcf呼叫 从gVCF文件调用变体。 给定一个参考面板（归一化的插入缺失，多等位基因分解）和一个单样本gVCF文件基因型，根据gVCF文件中的证据对参考面板中的位置进行基因分型。 gVCF文件不应分解。 假设条件 如果变体“匹配”，则面板变体用样本gVCF文件中的基因型调用（GT）和似然性（PL）进行注释。 面板变型与参考块重叠。 参照区中的每个面板变体都将具有纯合的参照基因型，并且基于全部<NON> ALT具有相同的PL。 面板变体与样本变体匹配。 当且仅当以下两种情况定义为匹配的两个变体： 它们在POS和REF上匹配 面板ALT是样本变体中最可能的ALT 要求 Python3 cyvcf2
GATK HaplotypeCaller SNP Calling 自动化流程
03-13
基于GATK3.3以上版本的HaplotypeCaller标准流程进行，测试脚本以sra文件为最初输入，将读取文件见内的所有sra文件（视为同一个样本的数据），进行SNP Calling产生GVCF，然后用户可自行进行jointCall
atlas_gatk4:在Atlas HPC上测试gatk4模块
04-07
Atlas HPC上的GATK4模块 在Atlas HPC上测试gatk4模块 示例数据集 重用GATK nextflow管道中的数据集（ ） wget https://iastate.box.com/shared/static/wt85l6s4nw4kycm2bo0gpgjq752osatu.gz tar -xf wt85l6s4nw4...
gatk-sv
03-19
GATK-SV 用于Illumina短读全基因组测序（WGS）数据的结构变异发现管道。 目录 原始呼叫者和证据收集 批量质量控制 -gCNV模型创建 批量证据合并，BAF生成和深度调用者 网站集群 网站指标 过滤 跨批次站点合并 基因...
gatk-cwl-generator:从GATK文档生成CWL文件
04-27
gatk-cwl生成器 从生成文件 安装 一，安装模块 git clone https://github.com/wtsi-hgi/gatk-cwl-generator cd gatk-cwl-generator python setup.py install 您可能还需要安装来运行生成的CWL文件 要求 的Python 3 ...
GATK3.2.2小结
浮生终有醒
08-02 6123
经过几天的摸索和网上资料的查询对GATK软件有点小心得，现总结如下： 1. fasta文件最好用定位到染色体上的数据，可以不用注释VCF文件（GVF），但如果用VCF文件保证以下几个条件： 1）VCF染色体必须和fasta的染色体数目一致，顺序一致 2）VCF的位点必须从小到大排序 3）VCF的碱基有可能有其他符号，如“~”等，要去除干净 2. 做之前分别使用picard
用GATK进行二代测序数据 SNP Calling 流程：（三）GenomicsDBImport 的多样本变异检测
Gossie的博客
10-26 9096
1. Genomics Database 对于群体数据来说，多样本同时时行 SNP Calling 的准确度要优于单个样本的 SNP Calling. GATK3 的多样本 SNP Calling 功能是 CombineGVCFs，GATK4 新出了 GenomicsDBImport功能，官网建议它适合1000个样本以上的 SNP Calling，但是它的另一个优点是可扩展性，即随时可以向 database 里添加新的材料，以扩大群体数量，而不用对旧的数据再从头操作一次。 Tips：GenomicsDBIm
动植物全基因组重测序（GATK）
weixin_47799897的博客
04-05 406
o --outdir:输出路径 -t --threads:线程数。比对软件（bwa+samtools）质控软件（fastqc）gvcf文件转vcf文件。picard 去除重复。建立索引 index。samtools转换。变异检测（gatk）
利用gatk进行calling-snp
hgz2020的博客
12-04 4717
使用gatk软件进行calling-snp
哪里可以直接下载gatk gvcf数据
最新发布
10-29
可以在GATK官方网站上找到直接下载GATK gVCF数据的链接。在官方网站上，你可以找到GATK的下载页面，其中包含了GATK的各个版本和相应的数据文件。在下载页面上，你可以找到与你所需版本相匹配的gVCF数据文件，并通过点击相应的链接进行下载。下载完成后，你可以将下载的文件用于相应的研究或分析工作。 除了官方网站，你还可以在一些公共数据库或资源库中找到可用的GATK gVCF数据。例如，国际人类基因组计划（1000 Genomes Project）提供了一些基因组变异的公共数据，其中包括gVCF格式的数据。你可以访问该项目的网站，查找与你研究相关的数据，并进行下载和使用。 此外，一些科研机构或实验室也可能在其网站上提供可下载的GATK gVCF数据。你可以通过搜索相关的研究机构或实验室的网站，并查看他们是否提供有关数据的下载链接或说明。 需要注意的是，下载GATK gVCF数据时，你应该确保选择与你的研究目的相匹配的数据集，并遵循相应的数据使用规定和许可证。此外，如果数据来自于不同的资源或数据库，你可能还需要考虑数据的一致性和可靠性，并根据需要进行数据处理和清洗。
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值