[Bioperl-l] Re ad in a Qual File like a Fasta?

Robert Buels rmb32 at cornell.edu
Fri Jul 2 00:18:58 UTC 2010


You probably want to do something like:

use Bio::SeqIO;

my $qual_in = Bio::SeqIO->new(
    -format => 'qual',
    -file => $ARGV[0],
   );

while( my $q = $qual_in->next_seq ) {
    print "qual for ".$q->id." is: ";
    print "$_ " for @{ $q->qual };
    print "\n";
}


twaddlac wrote:
> Hi,
> 
>      I am having a bit of trouble reading in a qual file; I would like to
> input the qual file into a hash but I can't seem to find a method that
> returns the qual sequence. Any ideas? Thank you in advance!





More information about the Bioperl-l mailing list