建库用的参考基因组，和vcf里面的参考基因组不一致，你检查检查，要用同一套基因组，版本也有相同；
dadi
计算dadi时,用这个命令运行后出现错误：运行命令为：cd annovar
## 根据基因组创建annovar库
gtfToGenePred -genePredExt $gtf unknown_refGene.txt
retrieve_seq_from_fasta.pl --format refGene --seqfile $fa --outfile unknown_refGeneMrna.fa unknown_refGene.txt
table_annovar.pl $data/goat.recode.vcf ./ -buildver unknown -out snp -remove -protocol refGene -operation g -nastring . -vcfinput
