该问题已被锁定!
2
关注
1108
浏览

利用GATK的HaplotypeCaller生成GVCF一直报错怎么回事

查看全部 2 个回答

pagu 初级会员 用户来自于: 北京市
2018-09-27 20:49
我在运行HyplotypeCaller的时候遇到过类似的问题,当时的问题在于: [b]bam文件中没有加样本信息[/b]   如果是DNA的测序分析数据,在使用bwa进行比对的时候可以直接用-R参数来添加这些信息 软件中的说明如下:[code]-R STRComplete read group header line. ’\t’ can be used in STR and will be converted to a TAB in the output SAM. The read group ID will be attached to every read in the output. An example is ’@RG\tID:foo\tSM:bar’. [null] 中文解释: -R STR 完整的read group的头部,可以用 '\t' 作为分隔符, 在输出的SAM文件中被解释为制表符TAB. read group 的ID,会被添加到输出文件的每一个read的头部。[/code]   使用方法举例:[code]bwa mem -t 4 -M -R "\@RG\tID:{library}\tLB:{library}\tPL:Illumina\tPU:{sample}\tSM:{sample}\" ref.fa read1.fastq read2.fastq > mem-pe.sam 2> ./mem-pe.log[/code]   如果是RNA比对的数据,或者是使用bwa比对时候忘了添加头部文件的数据,可以使用gatk中的AddOrReplaceReadGroups来添加对应信息,具体使用方法参考[url]https://software.broadinstitute.org/gatk/documentation/tooldocs/4.0.0.0/picard_sam_AddOrReplaceReadGroups.php[/url]   加了头文件再试试看看能不能运行。     参考内容: 1. BWA 命令详解   http://starsyi.github.io/2016/05/24/BWA-%E5%91%BD%E4%BB%A4%E8%AF%A6%E8%A7%A3/  

关于作者

问题动态

发布时间
2018-09-25 20:48
更新时间
2018-09-28 22:27
关注人数
2 人关注

相关问题

求助GATK 说明书
GATK call snp跳过了好大一块
在仅有每个样本(有3个repeat)所有转录本的fpkm的情况下,如何利用cuffdiff进行差异基因分析。
使用deeptools生成ChIP-seq信号热图与谱图
GATK和Samtools联合寻找变异
如何利用利用TPM或者FPKM完成DESeq2完成的工作?
chipseq分析,利用deeptools的computeMatrix reference-point画样本的peak center信号图
运行gatk HaplotypeCaller的时候提示没有AVX加速

推荐内容

使用tophat2和bowtie1寻找环形RNA时报错
linux系统中使用fastqc报错
R进行GO时出现No gene can be mapped....
求助GATK 说明书
如何快速将基因名称转换为基因ID(拟南芥)
分析CRISPR 高通量筛选数据
科研/读博
bowtie2使用报错,求助啦
画heatmap出现错误提示
运行gatk HaplotypeCaller的时候提示没有AVX加速
All Rights Reserved Powered BY WeCenter V4.1.0 © 2024