Skip to content

Instantly share code, notes, and snippets.

@arupgsh
Created March 7, 2015 07:20
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
Star You must be signed in to star a gist
Save arupgsh/7c7be76a80a3ff55b8aa to your computer and use it in GitHub Desktop.
print "Enter PDB File Name:";
$pdbn=<STDIN>;
print "Enter Output PDB File Name:";
$pdbo=<STDIN>;
print "Enter Chain Para 1:";
$pdbc1=<STDIN>;
print "Enter Chain Para 2:";
$pdbc2=<STDIN>;
print "Enter RES Para 1:";
$pdbr1=<STDIN>;
print "Enter RES Para 2:";
$pdbr2=<STDIN>;
print "Enter Chain Name:";
chomp($pdb1=<STDIN>);
print "Amino Acid Name:";
chomp($pdb2=<STDIN>);
open(in,"../$pdbn");
open(out,">$pdbo");
@pdb=<in>;
for($i=0;$i<=$#pdb;$i++)
{
chomp $pdb[$i];
if($pdb[$i]=~/^ATOM/)
{
$chain=substr($pdb[$i],$pdbc1,$pdbc2);
$residue=substr($pdb[$i],$pdbr1,$pdbr2);
if($chain eq "$pdb1" && $residue eq "$pdb2")
{
print out "$pdb[$i] \n";
print "$pdb[$i] \n";
}
}
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment