[Biojava-l] question about BLAST output parsing

Jiachen Xue xuejiachen at gmail.com
Mon Mar 15 23:09:50 UTC 2010


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



More information about the Biojava-l mailing list