我在用自己的蛋白与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
。。。。。。
后面想可能是blast版本的问题
后来换用blast-2.2.26的blastall 运行(formatdb还是用2.2.20)
就有结果了。。。。
https://blog.sciencenet.cn/blog-786123-637119.html