[Biopython] NCBIXML: hit start and end
Wibowo Arindrarto
w.arindrarto at gmail.com
Wed Apr 24 03:04:02 EDT 2013
Hi Mic,
> How do I get hsp->start('hit') and hsp->end('hit') from the bioperl code in
> Biopython?
With NCBIXML, they should be hsp.sbjct_start and hsp.sbjct_end respectively.
> Why does blast_record.query appears immediately in sequence and not after
> the other two for loops has finished?
It may be because the first three queries in your BLAST XML results
(XA10_v3.0-snap.{1..3}) do not have any hits and hsps. Check with your
XML results to be sure.
Hope that helps :),
Bow
More information about the Biopython
mailing list