²é¿´: 516  |  »Ø¸´: 1

houqingzhen

гæ (³õÈëÎÄ̳)

[ÇóÖú] python-×Ô¶¯·ÃÎÊÍøÂ磬½«htm½á¹û´æÅ̵½±¾µØ

ÐÅϢѧÐÂÊÖÇóÖú£º

ÎÒÏÖÔÚÓиöÎĵµ£¬ÏëÖð¸ö¶ÁÈ¡ÎĵµÀïµÄÄÚÈÝ£¬È»ºóÔÚÍøÒ³Éϵã»÷£¬×îºó½«·µ»ØµÄhtml½á¹û´æÅÌ¡£
ÍøÒ³ÍøÖ·ÊÇ
http://www.ncbi.nlm.nih.gov/COG/old/xognitor.html£¬
ÐèҪѡÔñµÄÑ¡ÏîÊÇBeTs to 17 clades£¬µã»÷µÄÑ¡ÏîÊÇcompare to COGs£¬
×îºó½«·µ»ØµÄhtml ½á¹û´æÅÌ¡£
ÍøÉÏÓиö½Å±¾£¬µ«ÊÇÎÒ²»ÓÃNCBIWWWÕâ¸öÄ£¿é£¬Ï£ÍûÄÜÓÃurl


ÔÚÍøÉϲéÕÒµÄÏà¹Ø½Å±¾ÊÇ£º
from Bio import Fasta
from Bio.Blast import NCBIWWW

filename_for_blast = 'seq/all.fasta' # Òª×öblastµÄfastaÎļþ

parser = Fasta.RecordParser()

file_for_blast = file(filename_for_blast, 'r')
iterator = Fasta.Iterator(file_for_blast, parser)
while 1:
cur_record = iterator.next()
if cur_record is None:
break
print u'´ò¿ªÊä³öÎļþ%s' % cur_record.title+'.htm'
file_for_blast_result = file(cur_record.title+'.htm', 'w')
print 'Now BLAST >', cur_record.title
''' NCBIWWW.blastº¯ÊýĬÈϵÄblast²ÎÊýÊÇ:
expect='10', alignments='50', timeout=20
'''
try:
blast_result = NCBIWWW.blast('blastn', 'nr', cur_record, \
expect='0.001', alignments='1', timeout=1)
file_for_blast_result.write(blast_result.read())
print u'OK£¬¸ã¶¨Ò»¸ö'
except:
print u'³¬Ê±ÁË£¬×öÏÂÒ»¸ö'
file_for_blast_result.close()
file_for_blast.close()

print u'È«²¿¸ã¶¨!'
raw_input()
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

jumper158

½ð³æ (ÕýʽдÊÖ)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

Õâ¸ö½Å±¾ÓÐÎÊÌ⣬ÄãµÄprint Óï·¨´òµÄ¾ÍÓÐÎÊÌ⣬Õâ¸öÍøÖ·ÊÇ×öʲôÓõģ¿ÄÜÓÐÏêϸµÄÔ­°æÃ´£¿Ëµ²»¶¨ÎÒ¿ÉÒÔ°ïÉÏÄã
ÓÎ×ßÓÚÅ£AÓëÅ£CÖ®¼ä£¬Ë­µÄÇà´º²»¿à±ï
2Â¥2012-06-10 20:49:32
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ houqingzhen µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Çóµ÷¼Á +7 Ê®Èý¼ÓÓÍ 2026-03-21 7/350 2026-03-23 23:48 by ÈÈÇéɳĮ
[¿¼ÑÐ] Ò»Ö¾Ô¸¼ª´ó»¯Ñ§322Çóµ÷¼Á +3 17501029541 2026-03-23 4/200 2026-03-23 23:47 by Txy@872106
[¿¼ÑÐ] »·¾³Ñ§Ë¶288Çóµ÷¼Á +8 ƤƤƤ123456 2026-03-22 8/400 2026-03-23 23:47 by ÈÈÇéɳĮ
[¿¼ÑÐ] 0854 ¿¼Ñе÷¼Á ÕÐÉúÁË£¡AI ·½Ïò +4 pk3725069 2026-03-19 16/800 2026-03-23 23:09 by Íô£¡£¿£¡
[¿¼ÑÐ] ÉϺ£µçÁ¦´óѧ²ÄÁÏ·À»¤ÓëвÄÁÏÖØµãʵÑéÊÒÕÐÊÕµ÷¼ÁÑо¿Éú£¨²ÄÁÏ¡¢»¯Ñ§¡¢µç»¯Ñ§£¬»·¾³£© +3 ÎÒ°®Ñ§µç³Ø 2026-03-23 3/150 2026-03-23 17:16 by AZMK
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤¿¼Ñе÷¼Á +4 ‹üÈA 2026-03-22 4/200 2026-03-23 16:13 by Ò»ÐݸçFU
[¿¼ÑÐ] 328Çóµ÷¼Á +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[¿¼ÑÐ] ÇóÀÏʦÊÕÎÒ +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[¿¼ÑÐ] Ò»Ö¾Ô¸ Î÷±±´óѧ £¬070300»¯Ñ§Ñ§Ë¶£¬×Ü·Ö287£¬Ë«·ÇÒ»±¾£¬Çóµ÷¼Á¡£ +3 ³¿»èÏßÓëÐǺ£ 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[¿¼ÑÐ] 260Çóµ÷¼Á +3 ÖìÜÆÁÕ 2026-03-20 4/200 2026-03-22 15:12 by ÖìÜÆÁÕ
[¿¼ÑÐ] ¿¼Ñе÷¼Á +4 À´ºÃÔËÀ´À´À´ 2026-03-21 4/200 2026-03-22 12:15 by ÐÇ¿ÕÐÇÔÂ
[¿¼ÑÐ] 354Çóµ÷¼Á +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by ÈËÀ´Ê¢
[¿¼ÑÐ] 0805²ÄÁÏ320Çóµ÷¼Á +3 ÉÎïÓï 2026-03-20 3/150 2026-03-21 15:46 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] 085601µ÷¼Á 358·Ö +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[¿¼ÑÐ] 274Çóµ÷¼Á +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤ 322Çóµ÷¼Á +4 È»11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[¿¼ÑÐ] ²ÄÁÏѧ˶297ÒѹýËÄÁù¼¶Çóµ÷¼ÁÍÆ¼ö +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[ÂÛÎÄͶ¸å] ÉêÇë»Ø¸åÑÓÆÚÒ»¸öÔ£¬±à¼­Í¬ÒâÁË¡£µ«ÏµÍ³ÉϵÄʱ¼äû±ä£¬¸ø±à¼­ÓÖдÓʼþÁË£¬Ã»»Ø¸´ 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤Çóµ÷¼Á +6 Ð÷ÐÒÓë×Ó 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[¿¼²©] 26²©Ê¿ÉêÇë +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by ÇáËɲ»ÉÙËæ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û