玩蛇网提供最新Python编程技术信息以及Python资源下载!

下载序列并保存到文本中的方法

想要知道如何下载一些序列并保存到文本中的python方法,可以参考以下这段代码。下面这段代码示例是从NCBI下载基因序列,读取文件(每行一个ID)中的GeneID,下载其序列并保存到文件中的应用。

关注本篇文章的人同时还关注了:python open()函数文件操作方法

#!/usr/bin/env python

# -*- coding:utf-8 -*-

__date__    = '2008-10-30'


from Bio import Entrez

def read_id(file_name):
    '''从文件中读取GeneID'''
    id_array = []
    fh = open(file_name, 'r')
    lines = fh.readlines()

    for line in lines:
        id = line.strip()
        id_array.append(id)

    fh.close()

    id_array = ','.join(id_array)
    return id_array

def download_seq (id_array):
    '''根据GeneID下载相应格式的序列'''

    result_handle = Entrez.efetch(db="nucleotide", rettype="genbank",  id=id_array)
    result = result_handle.read()

    return result

def write_to_file(file_out_name, content):
    '''将序列写入文件中 '''
    fh = open(file_out_name, 'w')
    fh.write(content)
    fh.close()

def main():
    '''主控制程序'''
    file_name = 'id_list.txt'
    file_out_name = 'sequences.txt'
    id_array = read_id(file_name)
    result = download_seq(id_array)
    write_to_file(file_out_name, result)

#www.iplaypy.com

if __name__ == '__main__':
    main()

玩蛇网文章,转载请注明出处和文章网址:http://www.iplaypy.com/code/text-file/tf2612.html [复制]



我要小额赞助,鼓励作者写出更好的教程↓↓↓

玩蛇网Python QQ群,欢迎加入: ① 279974227 玩蛇网Python新手群
修订日期:2016年04月01日 - 09时00分26秒 发布自玩蛇网

我要分享到:
评论列表(网友评论仅供网友表达个人看法,并不表明本站同意其观点或证实其描述)
相关文章推荐
别人正在看
特别推荐
站长推荐:
去顶部去底部