[Bioperl-l] .qual are like fasta files

Chad Matsalla matsallac@EM.AGR.CA
Fri, 03 Aug 2001 19:05:42 -0400


Hi all,

I need to parse some files containing quality values. They are each named
<whatever>.qual and look very similar to multiple fasta files:

>ML4861F PHD_FILE: ML4861F.phd.1
6 6 9 9 9 11 10 9 12 11 9 16 14 14 13 13 16 23 17 22
20 20 22 33 44 33 22 13 13 16 25 10 11 11 24 24 24
10 10 10 13 17 19 17 11 11 10 11 9 11 10 10 23 11 11
>ML4861F PHD_FILE: ML4861R.phd.1
6 6 9 9 9 11 10 9 12 11 9 16 14 14 13 13 16 23 17 22
20 20 22 33 44 33 22 13 13 16 25 10 11 11 24 24 24
10 10 10 13 17 19 17 11 11 10 11 9 11 10 10 23 11 11
>ML4861F PHD_FILE: ML4869F.phd.1
6 6 9 9 9 11 10 9 12 11 9 16 14 14 13 13 16 23 17 22
20 20 22 33 44 33 22 13 13 16 25 10 11 11 24 24 24
10 10 10 13 17 19 17 11 11 10 11 9 11 10 10 23 11 11

Is there a SeqIO::Fasta-like mechanism that I can use to go through these
files? Ideally I would like it to behave just like SeqIO::Fasta.

If not, I have a _lot_ to go through and I am willing to write something
if it could be used by somebody. This seems like something that
could nicely be used by my Phred.pm module.

Anyone? ... Anyone? ... Anyone?

Chad Matsalla
Agriculture & Agri-Food Canada
Saskatoon, Saskatchewan, Canada