<?xml version='1.0'?><rss version="2.0" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:georss="http://www.georss.org/georss" xmlns:atom="http://www.w3.org/2005/Atom" >
<channel>
	<title><![CDATA[BOL: Python script to convert fastq to fasta]]></title>
	<link>https://bioinformaticsonline.com/snippets/view/44493/python-script-to-convert-fastq-to-fasta?</link>
	<atom:link href="https://bioinformaticsonline.com/snippets/view/44493/python-script-to-convert-fastq-to-fasta?" rel="self" type="application/rss+xml" />
	<description><![CDATA[]]></description>
	
	<item>
	<guid isPermaLink="true">https://bioinformaticsonline.com/snippets/view/44493/python-script-to-convert-fastq-to-fasta</guid>
	<pubDate>Wed, 27 Mar 2024 20:30:32 -0500</pubDate>
	<link>https://bioinformaticsonline.com/snippets/view/44493/python-script-to-convert-fastq-to-fasta</link>
	<title><![CDATA[Python script to convert fastq to fasta]]></title>
	<description><![CDATA[<code>def fastq_to_fasta(fastq_file, fasta_file):
    with open(fastq_file, &#039;r&#039;) as fq:
        with open(fasta_file, &#039;w&#039;) as fa:
            while True:
                # Read four lines from the FASTQ file
                header = fq.readline().strip()
                sequence = fq.readline().strip()
                fq.readline()  # Skip the &#039;+&#039; line
                quality = fq.readline().strip()
                
                # Check for EOF
                if not header:
                    break
                
                # Write to the FASTA file
                fa.write(&#039;&gt;&#039; + header[1:] + &#039;\n&#039;)
                fa.write(sequence + &#039;\n&#039;)

# Usage example
fastq_to_fasta(&#039;input.fastq&#039;, &#039;output.fasta&#039;)</code>]]></description>
	<dc:creator>LEGE</dc:creator>
</item>

</channel>
</rss>