Skip to content

Commit

Permalink
Fix silly min vs max bug for genome chunking. D'oh!
Browse files Browse the repository at this point in the history
  • Loading branch information
tseemann committed Jun 22, 2015
1 parent 64bae17 commit daebebf
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion bin/snippy
Original file line number Diff line number Diff line change
Expand Up @@ -148,7 +148,7 @@ msg("Wrote $nfeat features to ref.gff");

my $refsize = -s "$refdir/ref.fa"; # rough size in bases
my $num_chunks = $cpus * 4; # oversample by 4 for run-time variation
my $chunk_size = min( $MIN_FREEBAYES_CHUNK_SIZE, int( $refsize / $num_chunks ) ); # bases per chunk
my $chunk_size = max( $MIN_FREEBAYES_CHUNK_SIZE, int( $refsize / $num_chunks ) ); # bases per chunk
msg("Freebayes will process $num_chunks chunks of $chunk_size bp, $cpus chunks at a time.");

# . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Expand Down

0 comments on commit daebebf

Please sign in to comment.