[BioPython] target sequence length in blast parsing

Ann Loraine aloraine at gmail.com
Sun Jan 7 06:23:19 UTC 2007


Dear all,

I have a question about blast parsing in biopython - any tips would be
much appreciated.

How can I access the length of the target sequence (e.g., 669 in the
following text) from alignment (or other?) objects retrieved from a
blast report parse?

** example **

>gi|34908012|ref|NP_915353.1| putative carboxypeptidase D [Oryza sativa (japo
nica
           cultivar-group)]
          Length = 669

 Score =  247 bits (625), Expect(2) = 3e-71
 Identities = 108/167 (64%), Positives = 132/167 (78%)
 Frame = +2

Yours,

Ann

-- 
Ann Loraine
Assistant Professor
Departments of Genetics, Biostatistics, and
Section on Statistical Genetics
University of Alabama at Birmingham
http://www.ssg.uab.edu
http://www.transvar.org



More information about the Biopython mailing list