R
Rita
Here is my code i want to do this work its related to Bioinformatics i
hope there are some bioinformatics mamber too:
my $result=$in->next_result;
my $hit =$result->next_hit;
my $hsp1 =$hit->next_hsp;
my$start1=$hsp1->query->start;
my $end1=$hsp1->query->end;
print "start First Hsp ",$start1,"\n";
print "End First Hsp ",$end1,"\n\n";
#Print Starting and Ending Point of 2 Hsps-
my $hsp2 =$hit->next_hsp;
my$start2=$hsp2->query->start;
my $end2=$hsp2->query->end;
print "start Second Hsp ",$start2,"\n";
print "End Second Hsp ",$end2,"\n\n";
#Print Gep Sequence
print "Gap Sequence- ",$seq1->subseq($end1+1,$start2-1);
you can see that i did same work twice bcoz of to get value of $end1
and $start2
how i can put in loopthis thing.i tried while loop but it's not working
..
thanks
hope there are some bioinformatics mamber too:
my $result=$in->next_result;
my $hit =$result->next_hit;
my $hsp1 =$hit->next_hsp;
my$start1=$hsp1->query->start;
my $end1=$hsp1->query->end;
print "start First Hsp ",$start1,"\n";
print "End First Hsp ",$end1,"\n\n";
#Print Starting and Ending Point of 2 Hsps-
my $hsp2 =$hit->next_hsp;
my$start2=$hsp2->query->start;
my $end2=$hsp2->query->end;
print "start Second Hsp ",$start2,"\n";
print "End Second Hsp ",$end2,"\n\n";
#Print Gep Sequence
print "Gap Sequence- ",$seq1->subseq($end1+1,$start2-1);
you can see that i did same work twice bcoz of to get value of $end1
and $start2
how i can put in loopthis thing.i tried while loop but it's not working
..
thanks