C*********m
发帖数: 213 | 1 不知道版上各位有没有碰到过:
blastall -p blastp -d wwpdb -e 0.01 -F F -m 8 -i a.fasta | awk '$3 > 90'
wwpdb 是pdb 里所有蛋白序列，ftp://ftp.wwpdb.org/pub/pdb/derived_data/pdb_
seqres.txt,a.fasta 是 NP_113671.3 序列 (862 a.a)，blast 2.2.15给出8个结果，
blast2.2.25只给出 7个,缺失的那个2p4r_T 只有25 a.a, 好像是 -F F 在 2.2.25版
本里不起作用． |
|
