[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