1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27
| #!/usr/bin/perl
use strict;
use warnings;
use Bio::SeqIO;
# fichier d'entrée
my $in = Bio::SeqIO->new(-file => "test.txt", '-format' => 'Fasta');
# les fichiers de sortie
my $outA= Bio::SeqIO->new(-file => ">sequencesA.fsa", '-format' => 'Fasta');
my $outB = Bio::SeqIO->new(-file => ">sequencesB.fsa", '-format' => 'Fasta');
my $outC = Bio::SeqIO->new(-file => ">sequencesC.fsa", '-format' => 'Fasta');
# récupération des séquences
while ( my $seq = $in->next_seq()){
if($seq->primary_id =~ /A/){
$outA->write_seq($seq);
}
elsif($seq->primary_id =~ /B/){
$outB->write_seq($seq);
}
elsif($seq->primary_id =~ /C/){
$outC->write_seq($seq);
}
} |