Skip to content
Snippets Groups Projects
Commit 0858b9cf authored by Andreas Kusalananda Kähäri's avatar Andreas Kusalananda Kähäri
Browse files

Avoid division by zero in generate_mapping_stats().

parent 23c5faec
No related branches found
No related tags found
No related merge requests found
......@@ -691,13 +691,35 @@ sub generate_mapping_stats {
my $novel_total = $stats->{'mapped_novel'} + $stats->{'lost_novel'};
# no split into known and novel for exons
unless ($type eq 'exon') {
$result .= sprintf($fmt2, 'known', $stats->{'mapped_known'},
$stats->{'lost_known'}, $stats->{'mapped_known'}/$known_total*100);
$result .= sprintf($fmt2, 'novel', $stats->{'mapped_novel'},
$stats->{'lost_novel'}, $stats->{'mapped_novel'}/$novel_total*100);
}
unless ( $type eq 'exon' ) {
if ( $known_total > 0 ) {
$result .= sprintf( $fmt2,
'known',
$stats->{'mapped_known'},
$stats->{'lost_known'},
$stats->{'mapped_known'}/$known_total*100 );
} else { # Avoid division by zero
$result .= sprintf( $fmt2,
'known',
$stats->{'mapped_known'},
$stats->{'lost_known'}, 0 );
}
if ( $novel_total > 0 ) {
$result .= sprintf( $fmt2,
'novel',
$stats->{'mapped_novel'},
$stats->{'lost_novel'},
$stats->{'mapped_novel'}/$novel_total*100 );
} else { # Avoid division by zero
$result .= sprintf( $fmt2,
'novel',
$stats->{'mapped_novel'},
$stats->{'lost_novel'}, 0 );
}
} ## end unless ( $type eq 'exon' )
$result .= sprintf($fmt2, 'total', $mapped_total, $lost_total,
$mapped_total/($known_total + $novel_total)*100);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment