Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Retrieve only top HSP per hit with blast+? erg55 Bioinformatics 2 06-24-2014 02:49 PM
how to identify blast hits with only 1 HSP (not limit the number of HSP) Kennels Bioinformatics 9 11-17-2013 10:17 PM
bioperl blast parsing problem Mark Bioinformatics 9 04-27-2013 04:22 PM
Blast table format hsp alignment !! empyrean Bioinformatics 0 03-17-2011 08:06 AM
Parsing BLAST results using BioPerl Ben Saville Bioinformatics 8 08-24-2010 07:43 AM

Thread Tools
Old 08-24-2015, 09:27 AM   #1
Junior Member
Location: Germany

Join Date: Aug 2015
Posts: 1
Thumbs up Blast HSP tiling BioPerl Module Query



I have one query regarding HSP tiling. I want to parse tblastn output to merge all overlapping HSPs and also non-overlaping HSP with some distance cutoff and would like to get every possible statics which blast tabular format provides.

What i have written untill now. but i dont have much experience with Bioperl blast tiling. Could you please help me with the same?

use Bio::SearchIO;

use Bio::Search::Tiling::MapTiling;
use Bio::Search::Tiling::TilingI;
open (FH,">trial_output.txt");
my $in = new Bio::SearchIO(-format => 'blast',
-file => 'output.tblastn');
while( my $result = $in->next_result ) {
## $result is a Bio::Search::Result::ResultI compliant object
while( my $hit = $result->next_hit ) {
## $hit is a Bio::Search::Hit::HitI compliant object

$tiling = Bio::Search::Tiling::MapTiling->new($hit);
my $ident = $tiling->identities('query','exact');
$query_length = $tiling->length('query');
#print $query_length."\t";
#print $ident."\t";
#print $result-> query_name,"\t";

my @contexts = $tiling->contexts('query');
my $best_context = pop @contexts;
my $max_frac = $tiling->frac_conserved(-type=>'query',

my $max_identical = $tiling->frac_identical(-type=>'query',

my $range = $tiling->range('query',$best_context);

#@alns = $tiling->get_tiled_alns('query', $best_context);

my @map = $tiling->coverage_map_as_text('query', $best_context);

foreach $m (@map)
print $m."\t";
print $max_frac."\t".$max_identical."\t".$range."\n";

swttalyan is offline   Reply With Quote

bioinformatics, bioperl, blast, hsp tiling

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 09:05 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO