Summary | Included libraries | Package variables | Synopsis | Description | General documentation | Methods |
WebCvs | Raw content |
use Bio::EnsEMBL::Analysis::Tools::GeneBuildUtils::ConfigDependent::TranscriptUtils qw( method );
or use Bio::EnsEMBL::Analysis::Tools::GeneBuildUtils::ConfigDependent::TranscriptUtils to get all methods
low_complexity_less_than_maximum | Description | Code |
low_complexity_less_than_maximum | code | next | Top |
Arg [1] : Bio::EnsEMBL::Transcript |
low_complexity_less_than_maximum | description | prev | next | Top |
my ($transcript, $complexity_threshold) = @_; my $peptide = $transcript->translate; my $seg = Bio::EnsEMBL::Analysis::Runnable::ProteinAnnotation::Seg->new ( -query => $peptide, -analysis => Bio::EnsEMBL::Analysis->new ( -logic_name => 'seg', -program_file => 'seg', ) ); $seg->run; my $low_complexity = $seg->get_low_complexity_length; logger_info(id($transcript)." has ".$low_complexity. " low complexity sequence"); #print_peptide($transcript);}
#print id($transcript)." has ".$low_complexity." low complexity sequence compared to".
# " ".$complexity_threshold."\n";
if($low_complexity >= $complexity_threshold){ warn(id($transcript)."'s low ". "complexity (".$low_complexity.") is above ". "the threshold ".$complexity_threshold. "\n"); return 0; } return 1; } 1;
CONTACT | Top |