²é¿´: 562  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[»ù½ðÉêÇë] »¶Ó­·¢À´filecodeµÄMz6ºóµÄ´úÂëÑéÖ¤Æä¹æÂÉ +31 ҽѧÀÏÄк¢ 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[»ù½ðÉêÇë] ÓÐʱºò£¬×ÔÈ»»ù½ðÕæµÄ²»ÄÜÌ«ÈÏÕæ £¨ÎÒµÄÉ걨¾­Ñ飩 +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[»ù½ðÉêÇë] ʱ¼ä´ÁÓÖ±äÁË8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[»ù½ðÉêÇë] filecode=ºóÃæµÚÒ»¸öÊÇ´óд×Öĸ +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[»ù½ðÉêÇë] ¹ØÓÚFilecode·ÖÎö·½·¨ +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[»ù½ðÉêÇë] ¸÷λµÀÓÑ£¬ÎÒҪȥÀ¥Ã÷Íæ¼¸Ì죬»ØÀ´¼û¡£ +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[»ù½ðÉêÇë] ÔÛÃÇÒ»ÆðÓÃÌúÖ¤·ÖÎö2026¹ú¼ÒÉç¿Æ»ù½ðÖбêÓë·ñ +7 ÆôÃȿƼ¼ 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[»ù½ðÉêÇë] ÊÇÕâÖܳö½á¹û»¹ÊÇÏÂÖܳö½á¹û£¿ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[»ù½ðÉêÇë] Ææ¹Ö£¬Á½¸öÈ˵Äfilecode¹Ì¶¨¶Î´ÓÍ·µ½Î²Ò»Ä£Ò»Ñù +8 ²¼²¼ºÍÒ»¶þ 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[»ù½ðÉêÇë] Ó¦¸ÃÊÇ93bebmhtakǰºóʮһ¸ö×Ö·û±È½Ï¹Ø¼ü +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[˶²©¼ÒÔ°] ¶Á²©µÄºÃ´¦ +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[»ù½ðÉêÇë] ÎҵĹú»ùÌáǰ֪µÀÖÐÁË£¬¿ÉÊÇͬʵIJÙ×÷ÈÃÎÒʵÔÚ½ÓÊܲ»ÁË£¬Ôõô»áÓÐÕâÑùµÄÈË +10 ¼ÒÓëÔ¶·½ 2026-08-10 15/750 2026-08-14 02:08 by ÃàÑò¸ç¸ç
[ÎÄѧ·¼²ÝÔ°] °¢ÒÌ +4 ÍôÍô¹ø 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[»ù½ðÉêÇë] ¾²µÈ»ù½ð½á¹û +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by ÇÒÌý»¢Ð¥
[»ù½ðÉêÇë] ÖØÒªÀ´Ô´£º±¾ÖÜÄ©³ö½á¹û +10 ˲ϢÓîÖæ 2026-08-12 10/500 2026-08-13 15:46 by likettle
[»ù½ðÉêÇë] ²»Ó¦¸Ã¿´fileCode +7 ÇÒÌý»¢Ð¥ 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[»ù½ðÉêÇë] ×ÛÊöÂÛÎÄ×÷Ϊ´ú±í×÷»á²»»áÓ°ÏìÆÀÉóר¼ÒµÄÓ¡Ïó·Ö£¿ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[»ù½ðÉêÇë] °ïæ¿´¿´fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by ±ùÐÄÓñºøÇç
[»ù½ðÉêÇë] ÎªÊ²Ã´ÍøÉϺܶàÈË˵±¾ÖÜ 12ºÅ³ö½á¹û +6 ˲ϢÓîÖæ 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[»ù½ðÉêÇë] È·¶¨ÁË£¬¹ú×ÔÈ»21ºÅ·Å°ñ +6 ²¼²¼ºÍÒ»¶þ 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û