我build立一个爆炸本地数据库。 但是,当我运行blastn命令时,我得到了这个错误消息:
T0“/home/coremake/release_build/build/PrepareRelease_Linux64-Centos_JSID_01_250088_130.14.22.10_9008__PrepareRelease_Linux64-Centos_1448906370/c++/compilers/unix/../../src/algo/winmask/seq_masker_istat_factory.cpp”,第170行:错误:ncbi :: CSeqMaskerIstatFactory :: DiscoverStatType() – 无法打开
T0“/home/coremake/release_build/build/PrepareRelease_Linux64-Centos_JSID_01_250088_130.14.22.10_9008__PrepareRelease_Linux64-Centos_1448906370/c++/compilers/unix/../../src/algo/winmask/seq_masker_istat_factory.cpp”,第271行:错误:ncbi :: CSeqMaskerIstatFactory :: create() – 无法创build单位计数容器
我正在使用这个命令来创build爆炸本地数据库:
makeblastdb -in chr23.fa -parse_seqids -dbtype nucl
这是我执行blastn的命令:
blastn -task megablast -db HumanGenome/blastdb/chr23.fa -window_masker_taxid 9606 -query readBatch.txt -out blastOut.txt
任何帮助将非常感激..谢谢
出现您需要首先创建WindowMasker文件:
https://www.ncbi.nlm.nih.gov/books/NBK279687/
具体而言,链接内容提到第1步:
https://www.ncbi.nlm.nih.gov/books/NBK279681/#cookbook.Create_masking_information_usin_1
如果windowmasker不是绝对的要求,可以考虑RepeatMasker
有多个安装步骤,但是你很简单。 我发现他们的文档比windowmasker文档更好,他们声称更好的掩盖覆盖。
(我的经验特别是免费的HMMER / Dfam选项。)
(我第二个Vince的关于创建掩码数据库的要求…你有没有尝试运行相同的blastn没有windowmasker选项?)