#!usr/bin/perl -w
use strict;
use warnings;
my @input_files=<*.seq>;
my $local_count=0;
my %hash;
foreach my $input_file(@input_files)
{
unless (open(INPUT, $input_file))
{
print "Cannot open file \"$input_file\"\n\n";
exit;
}
my $sequence='ABCD';
my @headings=('Tags', 'Frequency');
my $headings=join("\t",@headings);
while (my $line=)
{
if ($local_count==0){
my $outfile=$input_file;
$outfile=~s/.seq/.tag.txt/gi;
unless (open (OUTPUT, ">$outfile"))
{
print "Cannot open file \"$outfile\"\n\n";
exit;
}
}
chomp $line;
foreach($line=~m/$sequence/i){
if ($line=~m/$sequence(.{11})(.{11})$sequence/){
print OUTPUT "\n",$headings,"\n",$1,"\n",$2;
}
$local_count++;
}
}
}
####
Tags Frequency
CDDDDDDDDDD
BCDDEDDDDDR
Tags Frequency
CDEDEDDDESE
CEEESEEDESE
Tags Frequency