biopython文件处理:fastq文件转换为fasta文件

#!/usr/bin/python

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



from Bio import SeqIO

def fq2fa(my_file):

    with open(my_file) as handle:

        record=SeqIO.parse(handle, "fastq")

        SeqIO.write(record, "./new.fasta", "fasta")

    handle.close()    

 

posted @ 2017-12-11 14:35  遗世独立的愚公  阅读(2507)  评论(0编辑  收藏  举报