aw1:10
qs2:20
dd3:30
de4:10
hg5:30
dfd6:20
gf4:20
hgh5:30
hgy3:10
and so on...
####
>aw1
ATGCTAGATGCTAGCTAGCTAGCACTGAT
CGATGCTAGCGTAGTCAGCTGATGCTGTA
CGATGCTAGTCGTACG
>qs2
CGAGCTAGTCGTAGTCGTGATGCTGATTA
CGATGCTAGTCGTAGCTAGCTGATGCTGC
CGATGCTAGTCGTAGTC
>dd3
CGTAGTCGTAGTCGTAGTCGATGCTGATG
GCTAGTCGATGCTAGCTAGTCGATGCTGG
CGATGCTGAT
>de4
CGTAGTCGTAGTCGTACGTAGTCGTGAGT
CGATTATTTAGGAGGGACAAGGATAGTA
>hg5
CGTAGTCGTAGTCTAGTCGTGATGCTAGA
>dfd6
CGATGCTACGTACGTAGTCAGTCGTGATG
AATTAGAGCAGATAGAGGGGGAAAGGGTT
AAACCCC
>gf4
CGTAGTCAGTCTAGCTGATGTCGATGCTG
>hgh5
CATGCTAGTCGTAGTCGTAGTCGATGCTT
TTTTAAGGGAACCCCC
>hgy3
CCCCGGGTTTGGGAAAAGGGGGGGGATAG
####
open (FILE1, ">10.txt");
open (FILE2, ">20.txt");
open (FILE3, ">30.txt");
my @files = glob('Data/*.fa');
my %hash = &readFile();
foreach my $f(@files){
open FILE, $f or die "Cannot open : $!\n";
while (my $line = ){
chomp $line;
#print "$line\n";
if($line =~ /^>/){
my @head = split (/ /, $line);
my $name = substr($line, 1);
foreach my $key(keys %hash){
if($key eq $name){
if($hash{$key} == 10){
select FILE1;
}
if($hash{$key} == 20){
select FILE2;
}
if($hash{$key} == 30){
select FILE3;
}
}
}
}
}
}