[Bioperl-l] problem parsing msf file
Paola Bisignano
paola.bisignano at gmail.com
Fri Sep 4 12:28:03 UTC 2009
I have a problem with the parsing of msf file...I can't find the exact
object of Bio::SimpleAlign for my case...
I have to identify residues (from a list) in aligned sequences...but
when I parse the alignment from fasta file, I save as msf file, where
I have to identify my residue (from the list, numbering as the pdb
file) and the residue aligned in the aligned sequences...
this is a piece of the file...
NoName MSF: 2 Type: P Wed Aug 26 10:32:50 2009 Check: 00 ..
Name: Sequence/23-178 Len: 156 Check: 8937 Weight: 1.00
Name: 2zhz:A/1-148 Len: 156 Check: 9006 Weight: 1.00
//
1 50
Sequence/23-178 NDPRVAAYGE VDELNSWVGY TKSLINSHTQ VLSNELEEIQ QLLFDCGHDL
2zhz:A/1-148 DDARIAAIGD VDELNSQIGV L--LAEPLPD DVRAALSAIQ HDLFDLGGEL
51 100
Sequence/23-178 ATPADDERHS FKFKQEQPTV WLEEKIDNYT QVVPAVKKHI LPGGTQLASA
2zhz:A/1-148 CIPGHAAITD AHLARLDG-- WLA----HYN GQLPPLEEFI LPGGARGAAL
101 150
Sequence/23-178 LHVARTITRR AERQIVQLMR EEQINQDVLI FINRLSDYFF AAARYANYLE
2zhz:A/1-148 AHVCRTVCRR AERSIVALGA SEPLNAAPRR YVNRLSDLLF VLARVLNRAA
151 200
Sequence/23-178 QQPDML
2zhz:A/1-148 GGADVL
for example in this I have to identify the residue that is in front of
Val 28 (that is in Sequen) in 2zhz:A (that manually conting is Ile
5)....
Tyr4-> has no residue in front of it because the alignment starts from
N23 of Sequence...
how can I find the way to enter the residue of my sequen, and extract
the residue from the other????
I wish you all dear friends..and I'm actually in atrouble with this..
Thanks for suggestions
Paola
More information about the Bioperl-l
mailing list