Module 11

FASTQ

Reading quality data with SeqIO.

Quality data

SeqIO can also read FASTQ files. Each record includes the quality scores.

from Bio import SeqIO
from io import StringIO
fastq = "@seq1\nATGC\n+\nIIII"
handle = StringIO(fastq)
for record in SeqIO.parse(handle, "fastq"):
    print(record.id)
Your task

Read FASTQ

  • Parse a FASTQ string and print the record id.