***  For details on how to be removed from this list visit the  ***
***          CCP4 home page http://www.ccp4.ac.uk         ***


This perl snippet tells you where the gaps are. (It's a bit rough at the edges and ignores chain-id).

#!/usr/bin/perl
$n=0;
while(<>){
    if(($_=~/^ATOM/)||($_=~/^HETATM/)){
        $m=substr($_,22,4);
        unless(($n==$m)||($n==$m-1)){
            print "$n .. $m\n";
            };
        $n=$m;
        }
};
print "$m ..\n";

Cheers,
Charlie
--
Charlie Bond
Professorial Fellow
University of Western Australia
School of Biomedical, Biomolecular and Chemical Sciences
M310
35 Stirling Highway
Crawley WA 6009
Australia
[EMAIL PROTECTED]
+61 8 6488 4406

Reply via email to