[Bioperl-l] help parsing msf file or clustalW file reports
    Paola Bisignano 
    paola_bisignano at yahoo.it
       
    Tue Sep  1 12:20:25 UTC 2009
    
    
  
Hi, 
I'm trying to parse fasta files, where I have couple of alignments....I need to identify my residue in my alignment......I have separate lists that derived from ligplot parsing files.. so I have to manipulate string...but I don't now how to start..it seems complicated..
I used Bio::AlignIO to parse the fasta file, so I can have a parsed file in msf or clustalW forma
here an example:
CLUSTAL W(1.81) multiple sequence alignment
Sequence/9-273         DKWEMERTDITMKHKLGGGQYGEVYEGVWKKYSLTVAVKTLKEDTMEVEEFLKEAAVMKE
2pl0:A/6-268           DEWEVPRETLKLVERLGAGQFGEVWMGYYNGHT-KVAVKSLKQGSMSPDAFLAEANLMKQ
                       *:**: *  :.: .:**.**:***: * :: :: .****:**:.:*. : ** ** :**:
Sequence/9-273         IKHPNLVQLLGVCTREPPFYIITEFMTYGNLLDYLRECNRQEVSAVVLLYMATQISSAME
2pl0:A/6-268           LQHQRLVRLYAVVTQEP-IYIITEYMENGSLVDFLKTPSGIKLTINKLLDMAAQIAEGMA
                       ::* .**:* .* *:** :*****:*  *.*:*:*:  .  :::   ** **:**:..* 
I  choose two residue for example...how can I extract them...starting from their position in the pdb file?
I need to walk...to my sequence 
I don't know if it is clear because I cannot explain the question correctly in english...are there any Italians?
could anyone help me?
      
    
    
More information about the Bioperl-l
mailing list