[Bioperl-l] blat module
Jason Stajich
jason at cgt.duhs.duke.edu
Wed Sep 17 15:28:08 EDT 2003
You want:
my $blat_parser = new Bio::Tools::Blat(-file => "./temp");
Otherwise, in absence of any arguments that Root::IO knows about, it will
try and read from STDIN.
The new - more experimental Bio::SearchIO::blast is in CVS on the main
trunk if you want to give that whirl too. It *might* require you to have
stripped out the header line ahead of time, can't remember if I added in
a check to remove that or not.
-jason
On Wed, 17 Sep 2003, Jack Chen wrote:
> Hi All,
>
> I am trying to use the Blat.pm module. I have a simple script like this:
>
> #!/usr/bin/perl -w
> use strict;
> use Bio::Tools::Blat;
>
> my @bllat_feat;
>
> my $blat_parser = new Bio::Tools::Blat(-filename => "./temp");
>
> while( my $blat_feat = $blat_parser->next_result() ) {
> warn "here...\n";
> push @bllat_feat, $blat_feat;
> }
>
>
> In the script, temp is a blat result file, sitting in the current
> directory. The content of the file is (sorry that the lines in the file
> are wrapped). The probelm I have is that the program never runs. It gets
> stuck right after I execute the command. Any suggestion?
>
> Jack
>
> psLayout version 3
>
> match mis- rep. N's Q gap Q gap T gap T gap strand Q
> Q Q Q T T T T blockSizes
> qStarts tStarts
> match match count bases count bases
> name size start end name size start
> end
> ---------------------------------------------------------------------------------------------------------------------------------------------------------------
> 80 0 0 0 0 0 0 0 +
> origi 80 0 80 CHROMOSOME_III 13783268 13703494
> 13703574 1 80, 0, 13703494,
>
>
>
>
> ++++++++++++++++++++++++++++++++++++++++++++
> o-o Jack Chen, Stein Laboratory
> o---o Cold Spring Harbor Laboratory
> o----o #5 Williams, 1 Bungtown Road
> O----O Cold Spring Harbor, NY, 11724
> 0--o Tel: 1 516 367 8394
> O e-mail: chenn at cshl.org
> o-o Website: http://www.wormbase.org
> +++++++++++++++++++++++++++++++++++++++++++++
>
>
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
--
Jason Stajich
Duke University
jason at cgt.mc.duke.edu
More information about the Bioperl-l
mailing list