[Bioperl-l] Alignment Blast parser
brian_osborne at cognia.com
Thu Jan 22 14:23:57 EST 2004
I think he wants to print out the exact same alignment as you see in the
BLAST report, which you can't do using AlignIO.
Nathan, do you mean "exactly", with that "Query: 2" and "Subject: 2"
business as well?
From: bioperl-l-bounces at portal.open-bio.org
[mailto:bioperl-l-bounces at portal.open-bio.org]On Behalf Of Jason Stajich
Sent: Thursday, January 22, 2004 2:14 PM
To: Agrin, Nathan
Cc: bioperl-l at portal.open-bio.org
Subject: Re: [Bioperl-l] Alignment Blast parser
You want to retain it meaning what?
Bio::SearchIO(-format => 'blast')
parses blast for you.
On Thu, 22 Jan 2004, Agrin, Nathan wrote:
> I need a way to retain the EXACT alignment found in a blast file.
> Query: 1 MAGRGK-GKTSGKKAVSRSAKAGLQFPVGRIARYLKKGKYAERIGAGAPVYLAAVLEYLT
> M+GRGK G + KA SRS++AGLQFPVGR+ R L+KG YAER+GAGAPVYLAAVLEYLT
> Sbjct: 1 MSGRGKQGGKARAKAKSRSSRAGLQFPVGRVHRLLRKGNYAERVGAGAPVYLAAVLEYLT
> Query: 60 AEVLELAGNAARDNKKNRIVPRHIQLAIRNDEELGKLLGEVTIASGGVLPNIHAVLLP
> AE+LELAGNAARDNKK RI+PRH+QLAIRNDEEL KLLG VTIA GGVLPNI AVLLP
> Sbjct: 61 AEILELAGNAARDNKKTRIIPRHLQLAIRNDEELNKLLGRVTIAQGGVLPNIQAVLLP
> 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.
> Nathan Agrin
> Research Associate
> UMass Medical Center
> 55 Lake Ave. N.
> Worcester MA, 01655
> nathan.agrin at umassmed.edu
jason at cgt.mc.duke.edu
Bioperl-l mailing list
Bioperl-l at portal.open-bio.org
More information about the Bioperl-l