请问老师,使用annovar进行注释出现报错,如下:Error: invalid record found in exonic_variant_function file (exonic format error)

之前使用的时候发现此类报错,添加nopolish参数就可以运行,这一次添加完报错信息是:在输出文件时无法识别格式invalid record found in annovar outputfile:。是不是软件没有更新?annovar/2019Oct24/annotate_variation.pl 我用的是这一版。还是说有一些低质量没有过滤的snp导致了识别错误,我上一步的过滤是GATK硬标记完事之后,只用vcftools过滤了maf和maxmiss。

请先 登录 后评论

1 个回答

shidandan

Use of uninitialized value $1 in hash element at /mnt/data/home/tycloud/annovar/annovar/2019Oct24/annotate_variation.pl line 1445, <$QUERY> line 3623855.

NOTICE: Finished analyzing 2000000 query variants

Use of uninitialized value $1 in hash element at /mnt/data/home/tycloud/annovar/annovar/2019Oct24/annotate_variation.pl line 1445, <$QUERY> line 3623855.

NOTICE: Finished analyzing 3000000 query variants

NOTICE: Reading FASTA sequences from /mnt/data/home/tycloud/my_reseq/ref/unknown_refGeneMrna.fa ... Done with 20994 sequences

WARNING: A total of 49 sequences will be ignored due to lack of correct ORF annotation


NOTICE: Running with system command <coding_change.pl  /mnt/data/home/tycloud/my_reseq/annovar/snp.refGene.exonic_variant_function.orig /mnt/data/home/tycloud/my_reseq/ref/unknown_refGene.txt /mnt/data/home/tycloud/my_reseq/ref/unknown_refGeneMrna.fa -alltranscript -out /mnt/data/home/tycloud/my_reseq/annovar/snp.refGene.fa -newevf /mnt/data/home/tycloud/my_reseq/annovar/snp.refGene.exonic_variant_function>

Error: invalid record found in exonic_variant_function file (exonic format error):

这是不是软件版本问题啊?老师,很奇怪,上个月用的时候还可以运行,这个月一样报错,输入文件是同一个vcf文件,应该可以排除输入文件的问题。

请先 登录 后评论