Hi I tried the code....need debuging could you pl help..thank you.
#!/usr/bin/perl -w
$sasa_file= 'File1';
$pdbfile='File2';
# Remove the newline from the DNA filename
open (MYXFILE,$sasa_file);
open (PDBFILE, $pdbfile);
while (<MYXFILE>){
chomp;
@sa = split(/\t/,$_);
#print "$sa[0] $sa[1] $sa[2]\n";
#$res = substr $_, 0, 4;
#$res_no = substr $_, 4, 5;
#$atom =substr $_, 9, 3;
while (<PDBFILE>){
chomp;
$res = substr $_, 17, 3;
$res_no = substr $_, 22, 4;
$atom =substr $_, 13, 3;
#print " $res $res_no $atom \n" ; # works fine
if (( $sa[0]eq$res)&&($sa[1]==$res_no)&&($sa[2] eq $atom))
{
print "$_"; # does not print anything.
}
}
}
close MYXFILE;
close PDBFILE;