File: script1.pl

package info (click to toggle)
libtfbs-perl 0.5.svn.20100421-1
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 920 kB
  • sloc: perl: 6,073; ansic: 604; makefile: 30
file content (8 lines) | stat: -rw-r--r-- 322 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
8
 
#!/usr/bin/env perl  -w 
use Bio::DB::GenBank; 
use TFBS::DB::TRANSFAC; 
my $seq = Bio::DB::GenBank->new()->get_Seq_by_acc('AF100993'); 
my $db = TFBS::DB::TRANSFAC->connect(); 
my $pwm = $db->get_Matrix_by_ID('V$CEBPA_01','PWM'); 
my $siteset = $pwm->search_seq(-seqobj=>$seq, -threshold=>"80%"); 
print $siteset->GFF();