[Bioperl-l] Bio::SeqIO HOWTO

[chen li]

Hi folks,

 Here is one script copied from the Bio::SeqIO HOWTO:

     use Bio::SeqIO;
     my $usage = "getaccs.pl file format\n";
     my $file = shift or die $usage;
     my $format = shift or die $usage;

     my $inseq = Bio::SeqIO->new('-file'  => "<$file",
              '-format' => $format );
     while (my $seq = $inseq->next_seq) {
           print $seq->accession_number,"\n";
     }
     exit;


I have a small file called mouse.fasta kept in the
same directory. My question is that  how does the
script know to read in mouse.fasta? Where should I
make a small modification in the script?

Thanks,

Li
 


    


		
__________________________________ 
Yahoo! FareChase: Search multiple travel sites in one click.
http://farechase.yahoo.com