[Biopython-dev] [Bug 2767] Bio.SeqIO support for FASTQ and QUAL files

bugzilla-daemon at portal.open-bio.org bugzilla-daemon at portal.open-bio.org
Thu Feb 19 12:42:08 EST 2009


http://bugzilla.open-bio.org/show_bug.cgi?id=2767





------- Comment #1 from biopython-bugzilla at maubp.freeserve.co.uk  2009-02-19 12:42 EST -------
Created an attachment (id=1244)
 --> (http://bugzilla.open-bio.org/attachment.cgi?id=1244&action=view)
Read/write support for FASTQ and QUAL files, using the annotation dict

This patch stores the PHRED qualities as a list of integers in the SeqRecord's
annotations dictionary.  Changing this to use say a property, or a separate
per-letter-annotation dictionary should be trivial.

For QUAL files, the SeqRecord's seq is set to None.  This requires a few
changes to test_SeqIO.py which does not expect this.  We could also consider
introducing an UnknownSeq object (giving it a character like "?", "N", or "X",
an alphabet, and a length).  This would have a __str__ output like "?"*length.


-- 
Configure bugmail: http://bugzilla.open-bio.org/userprefs.cgi?tab=email
------- You are receiving this mail because: -------
You are the assignee for the bug, or are watching the assignee.


More information about the Biopython-dev mailing list