Dear Martin et al:
This, along with some scripting in coot and my favorite shell (zsh),
did the trick.
Many thanks to all who replied, and apologizes for what in retrospect
must have been a very poorly-worded question.
Maybe this functionality can be incorporated in the next release of
perlmole.
Bill Scott
On Jul 31, 2008, at 9:05 PM, Martin.Laurberg wrote:
> Bill,
> Almost there then (I wrote the line without checking)..
>
> I moved the 'P' one position to the left and now it works, provided
> first
> atom record of every atom record is the phosphor (/ P /) atom :
>
>
> perl -ne 'BEGIN{$i=0}if (/^ATOM/) \
> {if (substr($_,12,4) =~ / P /){$i++}; \
> substr($_,22,4,sprintf("%4.0f",$i))}{print}' file.pdb
>
>
>
>
>
> On Thu, 31 Jul 2008, William G. Scott wrote:
>
>> turned them all to zeros
>>
>>
>> William G. Scott
>>
>> Contact info:
>> http://chemistry.ucsc.edu/~wgscott/
>>
>>
>> On Jul 31, 2008, at 5:03 PM, Martin.Laurberg wrote:
>>
>>> Bill, if each residue starts at the same atom name, something
>>> close to
>>> this might work:
>>>
>>> perl -ne 'BEGIN{$i=0}if (/^ATOM/){if (substr($_,12,4) =~ / P /){$i+
>>> +};
>>> substr($_,22,4,sprintf("%4.0f",$i))}{print}' file.pdb
>>>
>>>
>>> /Martin
>>>
>>>
>>> --
>>> Martin Laurberg, PhD
>>> Noller Laboratory
>>> 225 Sinsheimer Laboratories
>>> University of California at Santa Cruz
>>> CA-95064 Santa Cruz
>>> USA
>>>
>>> Tel +1 (831) 459 35 84
>>> Fax +1 (831) 459 37 37
>>>
>>>
>>> On Thu, 31 Jul 2008, William Scott wrote:
>>>
>>>> Hi folks:
>>>>
>>>> I am hoping there is a simple answer I have overlooked to the
>>>> following question. I have a pdb file in it that has multiple
>>>> residues that have the same number and chain ID, and I want to
>>>> force
>>>> them to be renumbered sequentially. Is there a simple way to do so
>>>> (eg, pdbset) or am I doomed to fixing it manually?
>>>>
>>>> Thanks.
>>>>
>>>> Bill Scott
>>>>
>>>>
>>>> William G. Scott
>>>>
>>>> Contact info:
>>>> http://chemistry.ucsc.edu/~wgscott/
>>>>
>>
>>
|