博文
关于Blastall 时出现Selenocysteine (U) 却没最终结果的说明
|||
我在用自己的蛋白与swissprot数据库比对进行蛋白功能注释时，
出现如下警告，最终程序跑完，却没有最终比对结果：
swissprot用 blast-2.2.20 formatdb 格式化，
采用blast-2.2.20的blastall： 命令如下：
”blastall -p blastp -i protein.fasta -d uniprot_sprot.fasta -v 10 -b 10 -e 0.00001 -a 10 -m 8 -o protein_vs_swissprot.txt >log 2>errlog &“
报警结果：
Selenocysteine (U) at position 681 replaced by X
Selenocysteine (U) at position 696 replaced by X
Selenocysteine (U) at position 47 replaced by X
Selenocysteine (U) at position 47 replaced by X
Selenocysteine (U) at position 637 replaced by X
Selenocysteine (U) at position 690 replaced by X
Selenocysteine (U) at position 637 replaced by X
Selenocysteine (U) at position 690 replaced by X
Selenocysteine (U) at position 637 replaced by X
。。。。。。
查了下这个警告是正常的，其实也不算是警告！NCBI（http://www.ncbi.nlm.nih.gov/blast/blastcgihelp.shtml）上有说明，swissprot数据库里面有U。
后面想可能是blast版本的问题
后来换用blast-2.2.26的blastall 运行（formatdb还是用2.2.20）
就有结果了。。。。
https://blog.sciencenet.cn/blog-786123-637119.html
扫一扫，分享此博文