[Bioperl-l] Alignment Blast parser

Agrin, Nathan Nathan.Agrin at umassmed.edu
Thu Jan 22 12:47:26 EST 2004


I need a way to retain the EXACT alignment found in a blast file.
Ex-
 
 
Query: 1   MAGRGK-GKTSGKKAVSRSAKAGLQFPVGRIARYLKKGKYAERIGAGAPVYLAAVLEYLT
59
           M+GRGK G  +  KA SRS++AGLQFPVGR+ R L+KG YAER+GAGAPVYLAAVLEYLT
Sbjct: 1   MSGRGKQGGKARAKAKSRSSRAGLQFPVGRVHRLLRKGNYAERVGAGAPVYLAAVLEYLT
60
 
 
Query: 60  AEVLELAGNAARDNKKNRIVPRHIQLAIRNDEELGKLLGEVTIASGGVLPNIHAVLLP
117
           AE+LELAGNAARDNKK RI+PRH+QLAIRNDEEL KLLG VTIA GGVLPNI AVLLP
Sbjct: 61  AEILELAGNAARDNKKTRIIPRHLQLAIRNDEELNKLLGRVTIAQGGVLPNIQAVLLP
118
 
I have tried everything to no avail.  AlignIO almost works out, but
instead of doing a typical blast alignment, the best it can do is a
clustalw like alignment.
 
Any help is much appreciated.
 
-Nate
 
Nathan Agrin
Research Associate
UMass Medical Center
55 Lake Ave. N.
Worcester MA, 01655
(508)-856-6018
nathan.agrin at umassmed.edu
 


More information about the Bioperl-l mailing list