使用自建数据库。用ANNOVAR的annotate_variation.pl进行注释时,没有报错。
[*****@head ****]$ perl ../tessDSS/annovar/annotate_variation.pl -geneanno -buildver HP ./bin/M5-31.avinput HPdb/NOTICE: Output files are written to ./bin/M5-31.avinput.variant_function, ./bin/M5-31.avinput....
...8.0_linux_x86_64, 503.4GB RAM, 128 cores
我看了你们示例输出文件中的Homology.db_blast_results.txt.fnolowm50seqmGC,发现输出都是:
A C G T N totalnoN total00000000 00000000 00000000 00000000 00000000 00000000 00000000AT 00000000 GC 00000000
明显你们的示例数据也计...
#去掉序列ID中注释信息,避免错误:no ASV ids overlap between input FASTA and sequence abundance table
get_fa_by_id.pl $workdir/5.pick_otu_qiime/pick_de_novo_otus/otu_table_clean.txt $workdir/5.pick_otu_qiime/pick_de_novo_otus/rep_set/qiime_rep_set.fasta rep_set.farm -rf picrust2_out...