[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