#!/usr/bin/perl
# p e r l
# r u l e
# s b i g
# t i m e
use strict;
use warnings;
my %fields;
my ($base1, $base2);
my $linenr = 0;
my $pos = 14;
while(<>) {
chomp;
my $line = $_;
$base1 = substr $line,$pos,1;
if ( $base1 eq "t" ) {
pos($line) = 14;
substr($line,pos($line),1,"a");
}
if ( $base1 eq "a" ) {
pos($line) = 14;
substr($line,pos($line),1,"t");
}
pos($line) = 0;
$linenr++;
compare_line($line,$linenr);
}
sub compare_line {
my $line = shift;
my $linenr = shift;
open (DATA, "test_data.dat")
or die "Kan het bestand niet openen\n";
while(<DATA>) {
chomp;
$base2 = substr $line,$pos,1;
$base2 = "a" if $base2 eq "t";
$base2 = "t" if $base2 eq "a";
if ( $line eq $_ ) {
$fields{$linenr} = {
base1 => $base1,
base2 => $base2,
};
}
}
}
print "SPNL\tBASE1\tBASE2\n";
$pos += 1;
foreach my $linenr (sort keys %fields) {
print "$pos\t$fields{$linenr}{'base1'}\t$fields{$linenr}{'base2'}\n";
}