[Biojava-l] question about BLAST output parsing
Louise Ott
anjolou at hotmail.com
Tue Mar 16 09:23:37 UTC 2010
Sorry i forgot : there is an example of using blast parser in here :
http://biojava.org/wiki/BioJava:CookBook:Blast:Parser
It should be enough for what you want to do.
> Date: Mon, 15 Mar 2010 19:09:50 -0400
> From: xuejiachen at gmail.com
> To: biojava-l at lists.open-bio.org
> Subject: [Biojava-l] question about BLAST output parsing
>
> Hi,
>
> Thanks advance for help.
>
> For the following piece of text appearing in a blast output. How can I get
> the fields of "Identities", "Positives", "Gaps" as well as the alignment
> information, such as "DK V L+D + G +S + +++ +E GA+K+ L + AAPE" and
> subject string?
>
> >sp|B9KAQ6.1|UPP_THENN RecName: Full=Uracil phosphoribosyltransferase;
> AltName: Full=UMP
> pyrophosphorylase; AltName: Full=UPRTase
> Length = 209
>
> Score = 32.0 bits (71), Expect = 9.7, Method: Compositional matrix
> adjust.
> Identities = 16/42 (38%), Positives = 27/42 (64%), Gaps = 2/42 (4%)
>
> Query: 360 DKNVLLVDDSIVRGTTSEQIIEMAREAGAKKVYLAS--AAPE 399
> DK V L+D + G +S + +++ +E GA+K+ L + AAPE
> Sbjct: 124 DKEVFLLDPMLATGVSSVKALDILKENGARKITLVALIAAPE 165
> _______________________________________________
> Biojava-l mailing list - Biojava-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/biojava-l
_________________________________________________________________
Consultez gratuitement vos emails Orange, Gmail, Free, ... directement dans HOTMAIL !
http://www.windowslive.fr/hotmail/agregation/
More information about the Biojava-l
mailing list