Mercurial > repos > mvdbeek > damidseq_findpeaks

comparison find_peaks @ 12:7b0ca503bb95 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/damidseq_findpeaks commit 5bf4f4dbd82988d1612d951d69125bc5f5b439be-dirty
author mvdbeek
date Thu, 14 Jun 2018 11:14:52 -0400
parents be830200e987
children
comparison
equal deleted inserted replaced
11:c851dab0d2d9 12:7b0ca503bb95
408 ($chr, $start, $end, $score) = @line[0,3,4,5]; 408 ($chr, $start, $end, $score) = @line[0,3,4,5];
409 } 409 }
410 410
411 next unless $start; 411 next unless $start;
412 412
413 $score = 0 if $score eq "NA"; 413 $score = 0 if $score eq "NA" or $score eq 0.0;
414 414
415 push (@$in_ref, [$chr, $start, $end, $score]); 415 push (@$in_ref, [$chr, $start, $end, $score]);
416 } 416 }
417 417
418 close (IN); 418 close (IN);
434 my $total_coverage; 434 my $total_coverage;
435 435
436 foreach my $l (@$in_ref) { 436 foreach my $l (@$in_ref) {
437 my ($chr, $start, $end, $score) = @$l; 437 my ($chr, $start, $end, $score) = @$l;
438 next unless $score; 438 next unless $score;
439 $score = 0 if $score eq "NA";
440 $total_coverage += $end-$start; 439 $total_coverage += $end-$start;
441 push @frags, $score; 440 push @frags, $score;
442 } 441 }
443 442
444 @frags = sort {$a <=> $b} @frags; 443 @frags = sort {$a <=> $b} @frags;