[Bioperl-l] HSP tiling problem
biology0046 at hotmail.com
Mon Sep 15 16:00:02 EDT 2008
I try to use bioperl's SearhIO module to parse the blast output,
the output of blast contain many HSPs from the same hit, each HSP come out with an count of identity residues and alignment length, but what i want to get is the global identiy and alignment length of the query sequence and hit sequence, which may be done by tiling the HSPs and build a HSP contig,
the bioperl BLASTutil module contain "hsp_tiling and length_aln, frac_identical" method fullfill my purpose, but i read through the mailing list, some user said the method used by bioperl were not precise in many cases, they sugguest to use WUBLAST.
but i don't know which WUBLAST opinion should i use, i use "links" opinion, but the output still give me lots of HSPs, are these HSPs created from some small HSPs?best regards!wenkaichinese academy of sciences
More information about the Bioperl-l