玖叶教程网

情景描述，很多时候，我们拿到的往往是fastq的数据，但是我们需要转换成fasta数据进行后续操作，此时，我们需要转换文件格式，显然手动不可能完成，因此，需要用到编程技术。

#!/usr/bin/python
#-*- coding:utf-8 -*-

from Bio import SeqIO

def fq2fa(my_file,name):

    with open(my_file) as handle:

        record=SeqIO.parse(handle, "fastq")
        filename=name
        SeqIO.write(record, filename+".fasta", "fasta")

    handle.close()  

fq2fa(fastqfile,filename)

biopython包提供的SeqIO.parse提供了转换方案。