my @l= split /\t/, $_;
next unless length $l[2];
$l[3]='1001' if $l[3] eq '>1000';
- my $k= sprintf "%-23s\t%-31s", $l[0], $l[1];
- $s{$k}[$ab*2]= $l[2];
- $s{$k}[$ab*2+1]= $l[3];
+ my $k= sprintf "%5d\t%-23s", $l[2], $l[0];
+ $s{$k} ||= [ 0,0 ];
+ $s{$k}[$ab]= $l[3];
}
close F or die $!;
}
my @r= @{ $s{$k} };
my $whynot=
!defined($r[0]) ? 'ADD' :
- !defined($r[2]) ? 'DEL' :
- $r[0] != $r[2] ? 'pri' :
- $r[1] != $r[3] ? (sprintf "qt %4d",($r[3]-$r[1])) :
+ !defined($r[1]) ? 'DEL' :
+ $r[0] != $r[1] ? (sprintf "qt %4d",($r[1]-$r[0])) :
undef;
next unless $whynot;
print $whynot,"\t",$k;