gatk-4.1.9.0/gatk VariantFiltration warnings

15:55:13.094 WARN JexlEngine - ![0,14]: 'ReadPosRankSum < -8.0;' undefined variable ReadPosRankSum

15:55:13.094 WARN  JexlEngine - ![0,9]: 'MQRankSum < -12.5;' undefined variable MQRankSum

15:55:13.094 WARN  JexlEngine - ![0,14]: 'ReadPosRankSum < -8.0;' undefined variable ReadPosRankSum

15:55:13.094 WARN  JexlEngine - ![0,9]: 'MQRankSum < -12.5;' undefined variable MQRankSum



你的vcf文件中,有的行INFO那一列没有“MQRankSum” or "ReadPosRankSum"信息,所以才会出现这样的警告,

可以用“grep -v MQRankSum your_vcf”查找出来,这不是运行错误,可以不管继续进行后续分析;

若想要去除这个警告可以这样做

/home/software/gatk-4.1.9.0/gatk VariantFiltration --variant GenotypeGVCFs_SNP.vcf.gz --output GenotypeGVCFs_filter_SNP.vcf.gz -R /home/definegene/test/yelvti/call_snp/00.data/ref.fa -filter "(vc.hasAttribute('QD')&&QD<2.0)||MQ<40.0||(vc.hasAttribute('ReadPosRankSum')&&ReadPosRankSum<-8.0)||(vc.hasAttribute('MQRankSum')&&MQRankSum<-12.5)" --filter-name filter

/home/software/gatk-4.1.9.0/gatk VariantFiltration --variant GenotypeGVCFs_SNP.vcf.gz --output GenotypeGVCFs_filter_SNP.vcf.gz -filter "QD<2.0" --filter-name QD  -filter "MQ<40.0" --filter-name MQ  -filter "FS>60.0" --filter-name FS -filter "SOR>4.0" --filter-name SOR -filter "(vc.hasAttribute('MQRankSum')&&MQRankSum < -12.5)" --filter-name MQRankSum -filter "QUAL<30.0" --filter-name QUAL -filter "(vc.hasAttribute('ReadPosRankSum')&&ReadPosRankSum < -8.0)" --filter-name ReadPosRankSum

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • ####下载GRCm38 snp和indel vcf文件 #### ftp://ftp-mouse.sanger....
    njmujjc阅读 5,998评论 1 1
  • 转载GATK4.0和全基因组数据分析实践(下) 原创:矿工碱基矿工2018-03-23 前言 在上一篇文章中我已经...
    商乙农林科技阅读 10,653评论 1 11
  • Part0背景知识 Q:什么是外显子测序呢?A:外显子组测序是指利用序列捕获或者靶向技术将全基因组外显子区域DNA...
    天秤座的机器狗阅读 13,643评论 5 63
  • Part1 数据下载 先去Korean Personal Genome Project下载了编号为KPGP-00...
    天秤座的机器狗阅读 19,289评论 5 97
  • 1. 要分析的序列名称中,一般不要有空格2. 准备Reference文件(需为fasta格式)及比对1) Inde...
    晓佥阅读 7,349评论 0 4

友情链接更多精彩内容