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

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

为什么被折叠? 0 个回复被折叠
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-28 22:27
加了还是不行,这次貌似是另外一种错误

关于作者

问题动态

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

相关问题

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

推荐内容

使用tophat2和bowtie1寻找环形RNA时报错
运行gatk HaplotypeCaller的时候提示没有AVX加速
如何处理fastqc报告中duplication level报错的问题
分析CRISPR 高通量筛选数据
bowtie2使用报错,求助啦
如何快速将基因名称转换为基因ID(拟南芥)
使用cutadapt去掉adaptor后,fastqc结果的sequence distribution报错
R进行GO时出现No gene can be mapped....
linux系统中使用fastqc报错
GATK call snp跳过了好大一块
All Rights Reserved Powered BY WeCenter V4.1.0 © 2024