搜索
查看: 3363|回复: 1

使用GATK的GenotypeGVCFs工具处理HaplotypeCaller结果为何会丢失id?

[复制链接]

1

主题

5

帖子

110

积分

注册会员

Rank: 2

积分
110
发表于 2017-7-26 18:06:00 | 显示全部楼层 |阅读模式
这是HaplotypeCaller的结果,里面有rs的id这是GenotypeGVCFs处理的结果,但是rs的编号却丢失了,有人知道这是怎么回事吗?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x



上一篇:[提问]biocLite函数安装各种库,出现SSL connect error
下一篇:从UCSC上下载人类基因组TSS(转录起始位点)区域bed文件实战
回复

使用道具 举报

1

主题

4

帖子

58

积分

注册会员

Rank: 2

积分
58
发表于 2017-9-15 15:44:40 | 显示全部楼层
本帖最后由 小龙 于 2017-9-15 15:52 编辑

可能是在GenotypeGVCFs中没有添加--dbsnp 参数吧。

[Shell] 纯文本查看 复制代码
 java -jar GenomeAnalysisTK.jar \
   -T GenotypeGVCFs \
   -R reference.fasta \
   --variant sample1.g.vcf \
   --variant sample2.g.vcf \
   --dbsnp dbsnp_138.b37.vcf\
   -o output.vcf
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树 ( 粤ICP备15016384号  

GMT+8, 2019-10-14 22:17 , Processed in 0.029774 second(s), 28 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.