Hello,
I would like to parse the file fastq in order to obtain a hash table with the key ID and value for the next 3 lines.
I use the Bio::SeqIO::fastq module but how to obtain the score line?
--input--
@HCDPQ1D0501
GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT.
+HCDPQ1D0501
!''*((((***+))%%%++)(%%%%).1***-+*''))**55CCF>>>>>>CCCCCCC65.....
I would like to parse the file fastq in order to obtain a hash table with the key ID and value for the next 3 lines.
I use the Bio::SeqIO::fastq module but how to obtain the score line?
Code:
my $seqio = Bio::SeqIO::fastq->new('-format'=>'fastq' , '-file'=>'test.fastq'); while((my $rec = $seqio->next_seq())) { my $name = $rec->id(); my $seq= $rec->seq ; }
@HCDPQ1D0501
GATTTGGGGTTCAAAGCAGTATCGATCAAATAGTAAATCCATTTGTTCAACTCACAGTTT.
+HCDPQ1D0501
!''*((((***+))%%%++)(%%%%).1***-+*''))**55CCF>>>>>>CCCCCCC65.....
Comment