在SNPcalling结束后进行合并gvcf文件的时候发现文件里只有双亲的信心,核查发现两个子代的gvcf文件的表头中两个池的名字都是感亲的名字这个怎么修改

在SNPcalling结束后进行合并gvcf文件的时候发现文件里只有双亲的信心,核查发现两个子代的gvcf文件的表头中两个池的名字都是感亲的名字这个怎么修改?或者是哪里出现了问题?

四个文件并不一样大小

attachments-2020-12-OpwBq1fz5fcd8f034d59e.pngattachments-2020-12-PLXrNG4a5fcd8f04e2992.png

感病亲本

attachments-2020-12-sfXPyktu5fcd8daac9e5d.png

attachments-2020-12-evVTfYym5fcd8d985a8e6.png

子代感池,名称那一栏,子代的名字变为了Sparent

attachments-2020-12-LytE9KMS5fcd8dd43b08b.pngattachments-2020-12-xKYJrIao5fcd8e513b6bc.png


子代抗池,里面抗池名称那一栏变为了Sparent

attachments-2020-12-W4rCfP1h5fcd8e02b6399.png

attachments-2020-12-WC34tNOA5fcd8e2b84f10.png

请先 登录 后评论

最佳答案 2020-12-07 22:41

GATK是自动从bam文件中识别样本ID的,你做bwa比对的时候注意设置 -R 选项设置样本名字:

-R "@RG\tID:$sample\tLB:$sample\tSM:$sample\tPL:illumina\tPU:$sample"
请先 登录 后评论

其它 0 个回答

  • 1 关注
  • 0 收藏,2647 浏览
  • 杨辉 提出于 2020-12-07 10:10