[Bioperl-l] Bio::index::Fasta- where's the indexed file?
    zhihua li 
    lzhtom at hotmail.com
       
    Fri Feb  2 03:27:39 UTC 2007
    
    
  
Sorry guys, the former empty mail was sent out by mistake.
I'm using Bio::index::Fasta to index a file containing lots of sequences in 
fasta format. All is fine except one thing.
According to the bioperl tutorial and the documents, the following code 
will make a indexed file:
my $inx = Bio::Index::Fasta->new(-filename => "test.fasta.idx",
                                     -write_flag => 1);
    $inx->make_index("test.fasta");
And in another script I can access the indexed file by sayinig
$ENV{BIOPERL_INDEX} = "."; # find index in current directory
 my $inx = Bio::Index::Fasta->new(-filename => "test.fasta.idx");
my $seq=$inx->fetch("ent1001");        #fetch the sequence named ent1001
However, after running the first script, I cannot find a new file 
test.fasta.idx in my current directory. And not surprisingly, when I ran 
the second script, perl told me it couldn't find "test.fasta.idx".
What's going on here?
Thanks a lot!
_________________________________________________________________
与联机的朋友进行交流,请使用 MSN Messenger:  http://messenger.msn.com/cn  
    
    
More information about the Bioperl-l
mailing list