Hi there! I have small problem.
Within BioPerl I use SearchIO to parse my BLAST-output. Within this, the output is a result, which contain Hits, which contain HSPs.
How can I compare two HSPs with each other?
The two HSPs i want to compare are the current one and the ones before. If the current one is longer than the one before and is on the same place on the Query as the HSP before, i want to keep the current HSP and replace the old one (e.g. in an array).
I want to do that for all hits for one query entry, so that i have a maximum of one hit on one range of the query.
Thanks for any help in advance!
For reference: http://www.bioperl.org/wiki/HOWTO:SearchIO
Within BioPerl I use SearchIO to parse my BLAST-output. Within this, the output is a result, which contain Hits, which contain HSPs.
How can I compare two HSPs with each other?
The two HSPs i want to compare are the current one and the ones before. If the current one is longer than the one before and is on the same place on the Query as the HSP before, i want to keep the current HSP and replace the old one (e.g. in an array).
I want to do that for all hits for one query entry, so that i have a maximum of one hit on one range of the query.
Thanks for any help in advance!
For reference: http://www.bioperl.org/wiki/HOWTO:SearchIO
Comment