为什么不用GATK的best practise 呢? 这个soapsnp只有BGI那些人才用的吧
命令:~/biosoft/SOAPsnpZ/soapsnp -i chr01.sorted -d /disk3/wangshaobo/refnippon/IRGSPb5.fa -o chr01.snp.cns
错误信息:Reading Chromosome and dbSNP information Done.
terminate called after throwing an instance of 'std::out_of_range' what(): basic_string::substr: __pos (which is 63) > this->size() (which is 45) Aborted (core dumped)