有些参考基因组里面有一些 ambiguous 碱基 ,ANNOVAR处理不了,这是作者的回答:https://annovar.openbioinformatics.org/en/latest/misc/faq/
解决办法就是将 ambiguous 碱基根据IUPAC :(https://www.omicsclass.com/article/409 )编码替换成一单个碱基;
具体报错信息如下:
NOTICE: Running with system command <coding_change.pl ./indel_AT899.refGene.exonic_variant_function.orig /data/home/zhangzhi/genome-9-8//unknown_refGene.txt /data/home/zhangzhi/genome-9-8//unknown_refGeneMrna.fa -alltranscript -out ./indel_AT899.refGene.fa -newevf ./indel_AT899.refGene.exonic_variant_function>
Died at /data/home/zhangzhi/sft/pkg/annovar/coding_change.pl line 677.
Error running system command: <coding_change.pl ./indel_AT899.refGene.exonic_variant_function.orig /data/home/zhangzhi/genome-9-8//unknown_refGene.txt /data/home/zhangzhi/genome-9-8//unknown_refGeneMrna.fa -alltranscript -out ./indel_AT899.refGene.fa -newevf ./indel_AT899.refGene.exonic_variant_function>
Error running system command: </data/home/zhangzhi/sft/pkg/annovar/table_annovar.pl ./indel_AT899.avinput /data/home/zhangzhi/genome-9-8/ -buildver unknown -outfile ./indel_AT899 -remove -protocol refGene -operation g -nastring . -otherinfo>