Hi folks,
I was trying to recover a pdb of active site residues and close waters
(using ptraj from AmberTools1.2) with the following command:
strip "!((:1,22,144,145) | (.2478 <:4.2 & :WAT))"
This works fine in ambmask when checked without the negating !() to
keep the size down. I.e, ambmask should give the identify of the
residues left after the ptraj strip command, and ambmask does identify
the desired residues. However, ptraj drops core and gives a nonsense
pdb file (includes a bunch of seemingly random solute and water
atoms). Checkmask in rdparm just drops core with a segmentation fault
after printing a whole bunch of atom strings either with or without
the negating !().
Is anything obviously wrong with this? Has this behavior been
documented?
I tried the search engine on the archive, but may have failed due to
lack of search string foo.
Bud Dodson
PS,
strip "!(:1,22,144,145)"
works fine
strip "!(.2478 <:4.2 & :WAT)"
behaves the same as the whole strip mask given above.
--
M. L. Dodson
Business Email: activesitedynamics-at-comcast-dot-net
Personal Email: mldodson-at-comcast-dot-net
Phone: eight_three_two-56_three-386_one
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Fri Dec 05 2008 - 11:06:31 PST
