Entering edit mode
6.8 years ago
darya_boris
•
0
Hello!
I don't know, how I should do these tasks, please, help me!
Tasks:
#Per sequence, save information in a hash
#Go over hash, and print information of sequences one by one
Here is my script:
#!/usr/bin/perl
use strict;
use Bio::SeqIO;
my $sequence = "sequence.fa";
my $multifasta = Bio::SeqIO ->new (-file => "<$sequence",-format=> "fasta");
while (my $seq= $multifasta->next_seq()) {
my $na = $seq->display_id;
my $des = $seq->description;
my $ss = $seq->seq;
print "$na \n $des \n $ss \n_\n";
}
...and...what I should write next?
Thank you