1.

Can I Get Domain Number From Hmmpfam Or Hmmsearch Output?

Answer»

For EXAMPLE:

 SH2_5: domain 2 of 2, from 349 to 432: SCORE 104.4, E = 1.9e-26
Not directly but you can COMPUTE it since the domains are numbered by their order on the protein:
my @domains = $hit->domains;
my $domainnum = 1;
my $total = scalar @domains;
foreach my $domain ( SORT { $a->START <=> $b->start } $hit->domains ) {
print "domain $domainnum of $total,n";
$domainnum++;
}

For example:

 SH2_5: domain 2 of 2, from 349 to 432: score 104.4, E = 1.9e-26
Not directly but you can compute it since the domains are numbered by their order on the protein:
my @domains = $hit->domains;
my $domainnum = 1;
my $total = scalar @domains;
foreach my $domain ( sort { $a->start <=> $b->start } $hit->domains ) {
print "domain $domainnum of $total,n";
$domainnum++;
}


Discussion

No Comment Found

Related InterviewSolutions