Bonjour,

Comment récupérer la ligne ENTIERE de mon fichier fasta?

exemple:
>eb|xxx|xxx; xxxx xxxxx xxxx xxxxx xxxxx
catgtacgtgtagttagtagtacgtcagtacgtca
catgacgtcatgatcctgatcggatacggta

A l'aide de mon script, je récupère seulement: eb|xxx|xxx;
(je pense que cela est du à l'espace après le point virgule!)
Je voudrais récupérer: eb|xxx|xxx; xxxx xxxxx xxxx xxxxx xxxxx

Merci par avance

Mon script:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
my $file = 'xxx.fa ';
my $in  = Bio::SeqIO->new(-file => $file , '-format' => 'fasta');
my %hash;
my $nb_seq++;
 
while ( my $seq = $in->next_seq() ){
    $nb_seq++;          
    $hash{$seq->seq} = $seq->primary_id ;
}
foreach my $seq (keys %hash){   
    print ">$hash{$seq}\n$seq\n";
}