²é¿´: 543  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] *O_o^ +3 jk31rb6vyz 2026-06-12 3/150 2026-06-14 09:08 by 8bmv3h4hvh
[ÂÛÎÄͶ¸å] ¡Ñ¨Œ¡Ñ +3 m584r2z4j3 2026-06-11 3/150 2026-06-14 08:46 by 8bmv3h4hvh
[˶²©¼ÒÔ°] T^T +7 m584r2z4j3 2026-06-11 8/400 2026-06-14 08:34 by 8bmv3h4hvh
[¹«Åɳö¹ú] T^T +6 q1zg6x4lf3 2026-06-12 6/300 2026-06-14 08:28 by 8bmv3h4hvh
[¿¼²©] ¡ú_¡ú +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-14 08:28 by 8bmv3h4hvh
[¿¼²©] T^T +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-14 08:26 by 8bmv3h4hvh
[˶²©¼ÒÔ°] O(¡É_¡É)O +6 jk31rb6vyz 2026-06-12 8/400 2026-06-14 08:10 by 8bmv3h4hvh
[¹«Åɳö¹ú] ¦Ð_¦Ð +5 m584r2z4j3 2026-06-11 5/250 2026-06-14 08:08 by 8bmv3h4hvh
[¿¼ÑÐ] ¨q(¨s¦Å¨t)¨r +5 m584r2z4j3 2026-06-11 5/250 2026-06-14 08:03 by 8bmv3h4hvh
[»ù½ðÉêÇë] º¯ÆÀ +3 ²¼²¼ºÍÒ»¶þ 2026-06-13 3/150 2026-06-13 23:35 by ÓóÔ²ÏãÆ®Æ®
[ÕÒ¹¤×÷] ÊÛSCIÎÄÕ£¬ÎÒ:8.O.5.51O.5.4,¿ÆÄ¿ÆëÈ«,¿ÉÊ®¼± +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[¹«Åɳö¹ú] ‡Ö +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 19:05 by 1js5n6upxm
[¿¼ÑÐ] T^T +6 m584r2z4j3 2026-06-11 6/300 2026-06-13 18:29 by lkmi35jdzt
[ÎÄѧ·¼²ÝÔ°] ÏÌ²Ë +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[¿¼²©] ¦Ð_¦Ð +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[¿¼ÑÐ] o>_<o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[ÕÒ¹¤×÷] ÊÛSCIÎÄÕ£¬ÎÒ:8.O.5.51O.5.4,¿ÆÄ¿ÆëÈ«,¿ÉÊ®¼± +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[½Ìʦ֮¼Ò] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[¿¼²©] »¹ÓпÎÌâ×éÓв©Ê¿Ãû¶îÂð +5 ¶³ÓêÌì 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
[¿¼²©] ²©Ê¿ÉêÇë +3 ÕųÇì³ 2026-06-07 3/150 2026-06-08 09:15 by A´óÅóͬѧ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û