added get_seq

[biopieces.git] / bp_bin / format_genome

diff --git a/bp_bin/format_genome b/bp_bin/format_genome
index 35aa8b3af7354de14d65ef00aa444ebed767e149..2907aa4a61419aa31573d882cd0d2af192da76a0 100755 (executable)
--- a/bp_bin/format_genome
+++ b/bp_bin/format_genome
@@ -109,11 +109,13 @@ while ( $record = Maasha::Biopieces::get_record( $in ) )
     Maasha::Biopieces::put_record( $record, $out ) if not $options->{ "no_stream" };
 }
 
+close $fh_out if $fh_out;
+
 foreach $format ( @{ $options->{ 'formats' } } )
 {
     print STDERR qq(Creating format: $format for $genome ... ) if $options->{ 'verbose' };
 
-    if    ( $format =~ /^fasta$/i )     { Maasha::Fasta::fasta_index( "$fasta_dir/$genome.fna", "$dir/genomes/$genome/fasta/$genome.index" ) }
+    if    ( $format =~ /^fasta$/i )     { Maasha::Fasta::fasta_index( "$fasta_dir/$genome.fna", $fasta_dir, "$genome.index" ) }
     elsif ( $format =~ /^blast$/i )     { Maasha::NCBI::blast_index( "$genome.fna", $fasta_dir, "$dir/genomes/$genome/blast", "dna", $genome ) }
     elsif ( $format =~ /^blat$/i )      { warn "BLAT FORMAT NOT IMPLEMENTED" }
     elsif ( $format =~ /^vmatch$/i )    { Maasha::Match::vmatch_index( "$genome.fna", $fasta_dir, "$dir/genomes/$genome/vmatch", $tmp_dir ) }
@@ -124,8 +126,6 @@ foreach $format ( @{ $options->{ 'formats' } } )
     print STDERR qq(done.\n) if $options->{ 'verbose' };
 }
 
-close $fh_out if $fh_out;
-
 Maasha::Biopieces::close_stream( $in );
 Maasha::Biopieces::close_stream( $out );