13 Apr
2004
13 Apr
'04
12:09 p.m.
sorry to bother u with this simple question ...
I use following top script to extract sequences from pdb file '1PFZ'. works fine. But it reads all the different chain (A,B,C ..etc). what should i do if i want to get only sequence of chain A.
read_model FILE = '1PFZ', sequence_to_ali atom_files = '1PFZ', align_codes = 'gowtham' write_alignment file='seq.pir'
I tried the following and did not work read_model FILE = '1PFZ',MODEL_SEGMENT = '@:A' 'X:A'
#this works but writes only a chain break and * (/*) as a third line in output file 'seq.pir' many thanks in advance .... i could not get answer from list archives. i am sorry to bother you guys.
--
Ra. Gowthaman,
Graduate Student,
Bioinformatics Lab,
Malaria Research Group,
ICGEB , New Delhi.
INDIA
Phone: 91-9811261804
91-11-26173184; 91-11-26189360 #extn 314