²é¿´: 515  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©»¯¹¤´óѧ 070300 ѧ˶ 336·Ö Çóµ÷¼Á +6 vvÃÔ 2026-03-22 6/300 2026-03-23 22:00 by fuyu_
[¿¼ÑÐ] 298-Ò»Ö¾Ô¸Öйúũҵ´óѧ-Çóµ÷¼Á +10 ÊÖ»úÓû§ 2026-03-17 11/550 2026-03-23 16:30 by lingjue
[¿¼ÑÐ] 352Çóµ÷¼Á +3 ´óÃ×·¹£¡ 2026-03-22 3/150 2026-03-22 23:28 by king123£¡
[¿¼ÑÐ] 293Çóµ÷¼Á +3 ÌÎÌÎWjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[¿¼ÑÐ] 311Çóµ÷¼Á +6 ¶¬Ê®Èý 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[¿¼ÑÐ] Çóµ÷¼Á +6 Ê®Èý¼ÓÓÍ 2026-03-21 6/300 2026-03-22 17:00 by i_cooler
[¿¼ÑÐ] 308Çóµ÷¼Á +3 īīĮ 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[¿¼ÑÐ] 0856²ÄÁÏר˶353Çóµ÷¼Á +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÖпƼ¼´óѧ071000£¬Çóµ÷¼Á +4 ÑØ°¶Óб´¿Ç6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[¿¼ÑÐ] 286·ÖÈ˹¤ÖÇÄÜרҵÇëÇóµ÷¼ÁÔ¸Òâ¿ç¿¼£¡ +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[¿¼ÑÐ] »¯Ñ§µ÷¼Á +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[¿¼ÑÐ] ¡¾¿¼Ñе÷¼Á¡¿»¯Ñ§×¨Òµ 281·Ö£¬Ò»Ö¾Ô¸ËÄ´¨´óѧ£¬³ÏÐÄÇóµ÷¼Á +11 ³Ô³Ô³Ô²ÅÓÐÒâÒå 2026-03-19 11/550 2026-03-21 18:23 by ѧԱ8dgXkO
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖØÇì´óѧ085700×ÊÔ´Óë»·¾³×Ü·Ö308Çóµ÷¼Á +7 īīĮ 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[»ù½ðÉêÇë] ѧУÒѾ­Ìá½»µ½NSFC£¬»¹ÄÜÐÞ¸ÄÂ𣿠40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[¿¼ÑÐ] 265Çóµ÷¼Á +12 ÁºÁºÐ£Ð£ 2026-03-19 14/700 2026-03-21 13:38 by lature00
[¿¼ÑÐ] 317Çóµ÷¼Á +5 Éê×ÓÉêÉê 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷ÄϽ»Í¨ ר˶ ²ÄÁÏ355 ±¾¿ÆË«·Ç Çóµ÷¼Á +5 Î÷ÄϽ»Í¨×¨²Ä355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ ÄϾ©º½¿Õº½Ìì´óѧ´óѧ £¬080500²ÄÁÏ¿ÆÑ§Ó빤³Ìѧ˶ +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[¿¼ÑÐ] ÕÐÊÕµ÷¼Á˶ʿ +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[ÂÛÎÄͶ¸å] ÓÐûÓдóÀз¢Ð¡ÂÛÎÄÄÜ´øÎÒ¸ö¶þ×÷ +3 ÔöÈñ©ÈË 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û