[Bioperl-l] Directions for qual modules.
Wed, 08 Aug 2001 10:52:05 +0100
Chad Matsalla wrote:
> I know about the different formats that phred can write but at the moment
> I really only care about phd files:
> a 6 1
> c 6 20
> t 6 17
> and fasta-style files containing quality values only. These are what can
> be found in phrap'ed consed project directories.
> Should there be a Bio::SeqIO::phd and Bio::SeqIO::xbap and so on, one for
> each type of file? This sounds good to me and doesn't break the fact that
> at this time my Bio::SeqIO::qual only parses fasta-style quality files. It
> also seems more intuitive to do things this way rather then to pass in
> some flag when then SeqIO object is constructed.
> my $in_qual = Bio::SeqIO->new(-file => "<t/qualfile.qual" , '-format' =>
> my $in_qual = Bio::SeqIO->new(-file => "<t/phredfile.phd" , '-format' =>
> my $in_qual = Bio::SeqIO->new(-file => "<t/quality.xbap" , '-format' =>
> Heikki, how does this tie into your idea of a Bio::Seq::QualIO? How would
> this account for .phd files with both quality and sequence?
> Bio::Seq::QualIO might only be intuitively useful for files (like my
> fasta-quality files) with quality values only. I am really interested in
> phd files (like above) that have both quality and sequence.
I did not realise that the the sequence and quality are in the same
This just shows how long time it is I did any sequencing. ;-)
The above code example is fine. Go for it.
______ _/ _/_____________________________________________________
_/ _/ http://www.ebi.ac.uk/mutations/
_/ _/ _/ Heikki Lehvaslaiho email@example.com
_/_/_/_/_/ EMBL Outstation, European Bioinformatics Institute
_/ _/ _/ Wellcome Trust Genome Campus, Hinxton
_/ _/ _/ Cambs. CB10 1SD, United Kingdom
_/ Phone: +44 (0)1223 494 644 FAX: +44 (0)1223 494 468