Skip to content

Instantly share code, notes, and snippets.

@wizi
Created June 29, 2010 15:02
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save wizi/457329 to your computer and use it in GitHub Desktop.
Save wizi/457329 to your computer and use it in GitHub Desktop.
#!/stawberry/perl/bin/perl -w
use warnings;
use Bio::Tools::Run::Alignment::Blat;
use Bio::SeqIO;
use Bio::DB::Fasta;
my $seqio_object = Bio::SeqIO -> new (-file => 'C:/Users/MINH/Desktop/yeastquery.fasta', -format => 'fasta');
my $database = 'C:/Users/MINH/Desktop/yeast.nt';
my $db = Bio::DB::Fasta->new($database);
my $factory = Bio::Tools::Run::Alignment::Blat->new();
while (my $blat_object = $seqio_object -> next_seq) {
my $results = $factory -> run($blat_object,$db);
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment