[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[modeller_usage] regarding chain identifiers ...



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