Mercurial > repos > mvdbeek > damidseq_findpeaks
changeset 0:be830200e987 draft
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/damidseq_findpeaks commit 87c99a97cb2ce55640afdd2e55c8b3ae5ad99324
author | mvdbeek |
---|---|
date | Fri, 20 Apr 2018 04:34:17 -0400 |
parents | |
children | e726f7546561 |
files | damidseq_findpeaks.xml find_peaks test-data/hp1.bed test-data/peaks.gff |
diffstat | 4 files changed, 2744 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/damidseq_findpeaks.xml Fri Apr 20 04:34:17 2018 -0400 @@ -0,0 +1,68 @@ +<tool id="damidseq_find_peaks" name="damidseq find peaks" version="0.1.3"> + <description>Simple FDR random permutation peak caller</description> + <requirements> + <requirement type="package" version="5.26">perl</requirement> + </requirements> + <version_command><![CDATA['$__tool_directory__/find_peaks' --help 2>&1| grep find_peaks]]></version_command> + <command detect_errors="aggressive"><![CDATA[ +'$__tool_directory__/find_peaks' + --fdr=$fdr + --frac=$fdr + --min_count=$min_count + --min_quant=$min_quant + --n=$n + --step=$step + --unified_peaks=$unified_peaks + '$input_file' && +mv peak_analysis*/*.gff peaks.gff && +mv peak_analysis*/*data log.txt + ]]></command> + <inputs> + <param name="input_file" type="data" format="gff,bed,bedgraph" label="Select dam-fusion/dam ratio files" help="You can use damidseq_core to produce this file."/> + <param argument="--fdr" type="float" min="0" max="1" value="0.01" label="Set the False Discovery Rate (FDR)"/> + <param argument="--frac" type="float" min="0" max="1" value="0.01" label="Number of random fragments to consider per iteration"/> + <param argument="--min_count" type="integer" min="1" value="2" label="Minimum number of fragments to consider as a peak"/> + <param argument="--min_quant" type="float" min="0" max="1" value="0.95" label="Minimum quantile for considering peaks"/> + <param argument="--n" type="integer" min="1" value="100" label="Number of iterations"/> + <param argument="--step" type="float" min="0" max="1" value="0.01" label="Stepping for quantiles"/> + <param argument="--unified_peaks" type="select" label="Set the False Discovery Rate (FDR)"> + <option value="max">call maximum overlap as peak</option> + <option value="min">call minimum overlap as peak</option> + </param> + </inputs> + <outputs> + <data name="peaks" format="gff3" from_work_dir="peaks.gff" label="${tool.name} on ${on_string}"/> + <data name="logs" format="txt" hidden="true" from_work_dir="log.txt" label="${tool.name} logs on ${on_string}"/> + </outputs> + <tests> + <test> + <param name="input_file" value="hp1.bed" ftype="bed"/> + <param name="min_count" value="1"/> + <param name="fdr" value="0.1"/> + <output name="peaks" value="peaks.gff"/> + </test> + </tests> + <help><![CDATA[ + +Options: + --fdr False discovery rate value + [Current value: 0.01] + --frac Number of random fragments to consider per iteration + --min_count Minimum number of fragments to consider as a peak + [Current value: 2] + --min_quant Minimum quantile for considering peaks + [Current value: 0.95] + --n Number of iterations + [Current value: 100] + --step Stepping for quantiles + [Current value: 0.01] + --unified_peaks Method for calling peak overlaps (two options): + 'min': call minimum overlapping peak area + 'max': call maximum overlap as peak + [Current value: max] + + ]]></help> + <citations> + <citation type="doi">10.1093/bioinformatics/btv386</citation> + </citations> +</tool>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/find_peaks Fri Apr 20 04:34:17 2018 -0400 @@ -0,0 +1,633 @@ +#!/usr/bin/perl -w + +# Copyright © 2012-2015, Owen Marshall + +# This program is free software; you can redistribute it and/or modify +# it under the terms of the GNU General Public License as published by +# the Free Software Foundation; either version 2 of the License, or (at +# your option) any later version. +# +# This program is distributed in the hope that it will be useful, but +# WITHOUT ANY WARRANTY; without even the implied warranty of +# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +# General Public License for more details. +# +# You should have received a copy of the GNU General Public License +# along with this program; if not, write to the Free Software +# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 +# USA + +use strict; +use File::Basename; +use 5.010; +$|++; + +my $version = "1.0.1"; + +print STDERR "\nfind_peaks v$version\nCopyright © 2012-15, Owen Marshall\n\n"; + +my %vars = ( + 'fdr' => 0.01, + 'min_count' => 2, + 'n' => 100, + 'frac' => 0, + 'min_quant' => 0.95, + 'step' => 0.01, + 'unified_peaks' => 'max', +); + +my %vars_details = ( + 'fdr' => 'False discovery rate value', + 'min_count' => 'Minimum number of fragments to consider as a peak', + 'n' => 'Number of iterations', + 'frac' => 'Number of random fragments to consider per iteration', + 'min_quant' => 'Minimum quantile for considering peaks', + 'step' => 'Stepping for quantiles', + 'unified_peaks' => "Method for calling peak overlaps (two options):\n\r'min': call minimum overlapping peak area\n\r'max': call maximum overlap as peak", +); + + +my @in_files; +process_cli(); + +# Time and date +my ($sec,$min,$hour,$mday,$mon,$year) = localtime(); +my $date = sprintf("%04d-%02d-%02d.%02d-%02d-%02d",$year+1900,$mon+1,$mday,$hour,$min,$sec); + +help() unless @in_files; + +foreach my $fn (@in_files) { + my @in; + my @unified_peaks; + my @sig_peaks; + my @peakmins; + + my %peaks; + my %peak_count; + my %peak_count_real; + my %log_scores; + my %regression; + my %peak_fdr_cutoff; + my %fdr; + + # Output file names + my ($name,$dir,$ext) = fileparse($fn, qr/\.[^.]*/); + + # filenames + my $fn_base_date = "peak_analysis.".$name.".$date"; + my $base_dir = "$fn_base_date/"; + my $out = "$base_dir"."$name-FDR$vars{'fdr'}"; + my $out_peak_unified_track = $out.".peaks.gff"; + my $out_peaks = $out."_FDR-data"; + + # Load gff data files + load_gff(FILE=>$fn, IN_REF=>\@in); + + my $probes = @in; + + find_quants(IN_REF=>\@in, PEAKMINS_REF=>\@peakmins); + find_randomised_peaks(IN_REF=>\@in, PEAKMINS_REF=>\@peakmins, PEAK_COUNT=>\%peak_count, PEAKS=>\%peaks); + + # Make directory + mkdir($base_dir); + + # Open peaks file for writing + open(OUTP, ">$out_peaks")|| die "Cannot open peak file for writing: $!\n"; + print OUTP "FDR peak call v$version\n\n"; + print OUTP "Input file: $fn\n"; + + calculate_regressions(IN_REF=>\@in, PEAKMINS_REF=>\@peakmins, PEAK_COUNT=>\%peak_count, PEAKS=>\%peaks, LOG_SCORES=>\%log_scores, REGRESSION=>\%regression); + + call_peaks_unified_redux(ITER=>1, REAL=>1, AREF=>\@in, PEAKMINS_REF=>\@peakmins, PEAK_COUNT=>\%peak_count, PEAK_COUNT_REAL=>\%peak_count_real, PEAKS=>\%peaks); + + calculate_fdr(IN_REF=>\@in, PEAKMINS_REF=>\@peakmins, PEAK_COUNT=>\%peak_count, PEAK_COUNT_REAL=>\%peak_count_real, PEAK_FDR_CUTOFF=>\%peak_fdr_cutoff, FDR=>\%fdr, LOG_SCORES=>\%log_scores, REGRESSION=>\%regression); + + find_significant_peaks(PEAKMINS=>\@peakmins, SIG_PEAKS=>\@sig_peaks, PEAKS=>\%peaks, PEAK_FDR_CUTOFF=>\%peak_fdr_cutoff, FDR=>\%fdr); + + make_unified_peaks(SIG_PEAKS=>\@sig_peaks, UNIFIED_PEAKS=>\@unified_peaks, OUT=>$out_peak_unified_track, TYPE=>$vars{'unified_peaks'}); + + print STDERR "$#unified_peaks peaks found.\n\n"; + + close OUTP; +} + +print STDERR "All done.\n\n"; +exit 0; + + +###################### +# Subroutines start here +# + +sub find_significant_peaks { + my (%opts) = @_; + + my $peakmins = $opts{PEAKMINS}; + my $peaks = $opts{PEAKS}; + my $sig_peaks = $opts{SIG_PEAKS}; + my $fdr = $opts{FDR}; + my $peak_fdr_cutoff = $opts{PEAK_FDR_CUTOFF}; + + # Generate significant peaks and unify peaks + print STDERR "Selecting significant peaks ... \n"; + + foreach my $pm (@$peakmins) { + for my $i (0 .. $#{$$peaks{$pm}}) { + my ($chr, $pstart, $pend, $mean_pscore, $pscore, $count, $size) = @{ $$peaks{$pm}[$i] }; + if ($count >= $$peak_fdr_cutoff{$pm}) { + push (@$sig_peaks, [ @{$$peaks{$pm}[$i]}, $$fdr{$pm}[$count] ]); + } + } + } + + print OUTP "\nNumber of peaks: $#$sig_peaks\n"; +} + +sub calculate_fdr { + my (%opts) = @_; + my $in_ref = $opts{IN_REF}; + my $peakmins = $opts{PEAKMINS_REF}; + my $peak_count = $opts{PEAK_COUNT}; + my $log_scores = $opts{LOG_SCORES}; + my $regression = $opts{REGRESSION}; + my $fdr = $opts{FDR}; + my $peak_fdr_cutoff = $opts{PEAK_FDR_CUTOFF}; + my $peak_count_real = $opts{PEAK_COUNT_REAL}; + + foreach my $pm (@$peakmins) { + # get regression variables + my ($m,$b) = @{$$regression{$pm}} if $$regression{$pm}[0]; + + for my $i (0 .. $#{$$peak_count_real{$pm}}) { + next unless $$peak_count_real{$pm}[$i]; + my $expect = 10**($m*$i + $b); + + my $real_count = $$peak_count_real{$pm}[$i]; + my $fdr_conservative = $expect/$real_count; + $$fdr{$pm}[$i]= $fdr_conservative; + } + } + + # print FDR rates + print OUTP "\n"; + + foreach my $pm (@$peakmins) { + print OUTP "Peak min = $pm\n"; + for my $c (0 .. $#{$$fdr{$pm}}) { + next unless defined($$fdr{$pm}[$c]); + print OUTP "Peak size: $c\tCount: $$peak_count_real{$pm}[$c]\tFDR: $$fdr{$pm}[$c]\n"; + $$peak_fdr_cutoff{$pm} = $c if (($$fdr{$pm}[$c]<$vars{'fdr'}) && (!$$peak_fdr_cutoff{$pm})); + } + $$peak_fdr_cutoff{$pm}||= 10**10; # clumsy hack to prevent errors + print OUTP "\n"; + } + + foreach my $pm (@$peakmins) { + print OUTP "Peak min $pm: peak cutoff size for alpha = $vars{'fdr'} was $$peak_fdr_cutoff{$pm}\n\n"; + } +} + +sub calculate_regressions { + my (%opts) = @_; + my $in_ref = $opts{IN_REF}; + my $peakmins = $opts{PEAKMINS_REF}; + my $peaks = $opts{PEAKS}; + my $peak_count = $opts{PEAK_COUNT}; + my $log_scores = $opts{LOG_SCORES}; + my $regression = $opts{REGRESSION}; + + my $in_num = @$in_ref; + + foreach my $pm (@$peakmins) { + print OUTP "Peak min = $pm\n"; + for my $c (0 .. $#{$$peak_count{$pm}}) { + my $peak_count_avg = $$peak_count{$pm}[$c]/$vars{'n'} if $$peak_count{$pm}[$c]; + next unless $peak_count_avg; + + if ($vars{'frac'}) { + $peak_count_avg = $peak_count_avg * $in_num/$vars{'frac'}; + } + $$log_scores{$pm}[$c] = log($peak_count_avg)/log(10); + print OUTP "Peak size: $c\tCount:$peak_count_avg\n"; + } + + # calculate exponential decay rates + # y= a+bx for log(y) + my ($sumx, $sumy, $sumxsq, $sumxy); + my $n=0; + for my $i (0 .. $#{$$peak_count{$pm}}) { + next unless $$peak_count{$pm}[$i]; + $n++; + $sumx += $i; + $sumy += $$log_scores{$pm}[$i]; + $sumxsq += $i ** 2; + $sumxy += $i * $$log_scores{$pm}[$i]; + } + + next unless $n > 1; + + my $mean_x = $sumx/$n; + my $mean_y = $sumy/$n; + my $mean_xsq = $sumxsq/$n; + my $mean_xy = $sumxy/$n; + + my $b = ($mean_xy - ($mean_x * $mean_y)) / ($mean_xsq - ($mean_x **2)); + my $a = $mean_y - ($b * $mean_x); + + # store values + $$regression{$pm}=[$b, $a]; + print OUTP "regression: log(y) = $b(x) + $a\n"; + + for my $i (0 .. $#{$$peak_count{$pm}}) { + next unless $$peak_count{$pm}[$i]; + my ($b,$a) = @{$$regression{$pm}}; + my $logval = ($b*$i + $a); + my $val = 10**$logval; + print OUTP "lin regress: $i\t$$log_scores{$pm}[$i]\t$logval\t$val\n"; + } + + print OUTP "\n"; + } +} + +sub find_randomised_peaks { + my (%opts) = @_; + my $in_ref = $opts{IN_REF}; + my $peakmins_ref = $opts{PEAKMINS_REF}; + my $peaks = $opts{PEAKS}; + my $peak_count = $opts{PEAK_COUNT}; + + print STDERR "Duplicating ... \n"; + my @inr = @$in_ref; + + # Call peaks on input file + print STDERR "Calling peaks on input file ...\n"; + + for my $iter (1 .. $vars{'n'}) { + #print STDERR "Iteration $iter ... \r"; + print STDERR "Iteration $iter: [shuffling] \r"; + + if ($vars{'frac'}) { + # only use a fraction of the array per rep + my @a = inside_out(\@inr, $vars{'frac'}); + call_peaks_unified_redux(ITER=>$iter, AREF=>\@a, PEAKMINS_REF=>$peakmins_ref, PEAK_COUNT=>$peak_count); + } else { + shuffle(\@inr); + call_peaks_unified_redux(ITER=>$iter, AREF=>\@inr, PEAKMINS_REF=>$peakmins_ref, PEAK_COUNT=>$peak_count); + } + } + +} + +sub call_peaks_unified_redux { + my (%opts) = @_; + my $iter = $opts{ITER}; + my $real = $opts{REAL}; + my $a = $opts{AREF}; + my $peakmins_ref = $opts{PEAKMINS_REF}; + my $peak_count_real = $opts{PEAK_COUNT_REAL}; + my $peaks = $opts{PEAKS}; + my $peak_count = $opts{PEAK_COUNT}; + + my ($pstart, $pend, $inpeak, $pscore, $count); + $pstart=$pend=$pscore=$inpeak=$count=0; + + my @tmp_peak; + my $total = $#$a; + + if ($real) { + print STDERR "Calling real peaks ... \r"; + } else { + print STDERR "Iteration $iter: [processing ...] \r"; + } + + my $old_chr=""; + foreach my $pm (@$peakmins_ref) { + for my $i (0 .. $total) { + my ($chr, $start, $end, $score) = @{ @$a[$i] }; + next unless $score; + + if ($real) { + unless ($chr eq $old_chr) { + # Next chromosome + # (Peaks can't carry over chromosomes, but we don't use this shortcut when randomly shuffling) + $pstart=$pend=$pscore=$inpeak=$count=0; + @tmp_peak = () if $real; + } + } + $old_chr = $chr if $real; + + unless ($inpeak) { + next unless $score >= $pm; + # record new peak + $pstart = $start; + $pend = $end; + $pscore = $score * ($end-$start)/1000; + $count++; + push @tmp_peak, $score if $real; + $inpeak = 1; + } else { + if ($score >= $pm) { + # still in peak + $count++; + + # Fragment score to deal with scoring peaks made from uneven sized fragments + my $fragment_score = $score * ($end-$start)/1000; + + push @tmp_peak, $score if $real; + $pscore += $fragment_score; + $pend = $end; + } else { + # out of a peak + if ($count >= $vars{'min_count'}) { + # record peak + if ($real) { + $$peak_count_real{$pm}[$count]++; + my $mean_pscore = sprintf('%0.2f',($pscore/(($pend-$pstart)/1000))); + push (@{$$peaks{$pm}},[($chr, $pstart, $pend, $mean_pscore, $pscore, $count, ($pend-$pstart))]); + } else { + $$peak_count{$pm}[$count]++; + } + } + + # reset + $pstart=$pend=$pscore=$inpeak=$count=0; + @tmp_peak = () if $real; + } + } + } + } +} + + +sub shuffle { + # Fisher-Yates shuffle (Knuth shuffle) + my ($array) = @_; + my $i = @$array; + while ( --$i ) { + my $j = int rand( $i+1 ); + @$array[$i,$j] = @$array[$j,$i]; + } +} + +sub inside_out { + my ($array, $frac) = @_; + my @a; + for my $i (0 .. $frac-1) { + my $j = int rand( $i+1 ); + if ($j != $i) { + $a[$i] = $a[$j] + } + $a[$j] = @$array[$i] + } + return @a; +} + +sub load_gff { + my (%opts) = @_; + my $fn = $opts{FILE}; + my $in_ref = $opts{IN_REF}; + + print STDERR "Reading input file: $fn ...\n"; + open (IN, "<$fn") || die "Unable to read $fn: $!\n"; + + my $i; + while (<IN>) { + $i++; + print STDERR "Read $i lines ...\r" if $i%10000 == 0; + chomp; + + my @line = split('\t'); + + my ($chr, $start, $end, $score); + if ($#line == 3) { + # bedgraph + ($chr, $start, $end, $score) = @line; + } else { + # GFF + ($chr, $start, $end, $score) = @line[0,3,4,5]; + } + + next unless $start; + + $score = 0 if $score eq "NA"; + + push (@$in_ref, [$chr, $start, $end, $score]); + } + + close (IN); + + print STDERR "Sorting ... \n"; + @$in_ref = sort { $a->[1] <=> $b->[1] } @$in_ref; + @$in_ref = sort { $a->[0] cmp $b->[0] } @$in_ref; +} + +sub find_quants { + my (%opts) = @_; + + my $in_ref = $opts{IN_REF}; + my $peakmins_ref = $opts{PEAKMINS_REF}; + + my %seg; + my @frags; + + my $total_coverage; + + foreach my $l (@$in_ref) { + my ($chr, $start, $end, $score) = @$l; + next unless $score; + $score = 0 if $score eq "NA"; + $total_coverage += $end-$start; + push @frags, $score; + } + + @frags = sort {$a <=> $b} @frags; + + print STDERR "Total coverage was $total_coverage bp\n"; + + my @quants; + for (my $q=0;$q<=1;$q+=$vars{'step'}) { + push @quants, [$q, int($q * @frags)] if $q > $vars{'min_quant'}; + } + + foreach (@quants) { + my $cut_off = @{$_}[0]; + my $score = $frags[@{$_}[1]]; + printf(" Quantile %0.2f: %0.2f\n",$cut_off,$score); + $seg{$cut_off} = $score; + } + + foreach my $c (sort {$a <=> $b} keys %seg) { + push (@$peakmins_ref, $seg{$c}); + } +} + +sub make_unified_peaks { + my (%opts) = @_; + my $ref = $opts{SIG_PEAKS}; + my $out_peak_unified_track = $opts{OUT}; + my $unified_peaks = $opts{UNIFIED_PEAKS}; + my $type = $opts{TYPE}; + my $total = @$ref; + + # Unify overlapping peaks, and make significant peaks file + my $skipped_peaks; + print STDERR "Combining significant peaks ...\n"; + + # unroll chromosomes for speed + foreach my $chr (uniq( map @{$_}[0], @$ref )) { + my @c = grep {@{$_}[0] eq $chr} @$ref; + + my @unified_peaks_chr; + foreach my $ar (@c) { + if (state $i++ % 100 == 0) { + my $pc = sprintf("%0.2f",($i*100)/$total); + print STDERR "$pc\% processed ...\r"; + } + + my ($chra, $start, $end, $score, $total_score, $count, $peaklen, $fdr) = @$ar; + + # next if @unified_peaks_chr already overlaps + next if grep { + @{$_}[3] < $end + && @{$_}[4] > $start + } @unified_peaks_chr; + + # Grab all elements that overlap + my @test = grep { + @{$_}[1] < $end + && @{$_}[2] > $start + } @c; + + for my $j (0 .. $#test) { + my ($chr1, $start1, $end1, $score1, $total_score1, $count1, $peaklen1, $fdr1) = @{$test[$j]}; + + next unless $start1 < $end; + next unless $end1 > $start; + + if ($type eq 'min') { + $start = max($start, $start1); + $end = min($end, $end1); + } else { + $start = min($start, $start1); + $end = max($end, $end1); + } + + $score = max($score, $score1); + $fdr = min($fdr, $fdr1); + } + + push @unified_peaks_chr, [($chr, '.', '.', $start, $end, $score, '.', '.', "FDR=$fdr")]; + } + + @$unified_peaks = (@$unified_peaks, @unified_peaks_chr); + } + + $total = $#$unified_peaks; + + print STDERR "Sorting unified peaks ...\n"; + @$unified_peaks = sort { $a->[3] <=> $b->[3] } @$unified_peaks; + @$unified_peaks = sort { $a->[0] cmp $b->[0] } @$unified_peaks; + + print STDERR "Writing unified peaks file ...\n"; + open(PEAKOUTUNI, ">$out_peak_unified_track") || die "Unable to open peak output track for writing: $!\n\n"; + for my $j (0 .. $#$unified_peaks) { + print PEAKOUTUNI join("\t", @{$$unified_peaks[$j]}), "\n"; + } +} + +sub max { + my ($max, @vars) = @_; + my $index=0; + $max||=0; + for my $i (0..$#vars) { + ($max, $index) = ($vars[$i], $i+1) if $vars[$i] > $max; + } + return $max; +} + +sub min { + my ($min, @vars) = @_; + my $index=0; + $min||=0; + for my $i (0..$#vars) { + ($min, $index) = ($vars[$i],$i+1) if $vars[$i] < $min; + } + return $min; +} + +sub uniq { + my %seen; + return grep { !$seen{$_}++ } @_; +} + + + +sub process_cli { + foreach (@ARGV) { + if (/--(.*)=(.*)/) { + unless (defined($vars{$1})) { + print STDERR "Did not understand $_ ...\n"; + help(); + } + my ($v, $opt) = ($1,$2); + $vars{$v} = $opt; + next; + } elsif (/--h[elp]*/) { + help(); + } elsif (/--(.*)/) { + print STDERR "Please add a parameter to $_ ...\n\n"; + exit 1; + } + push @in_files, $_; + } +} + + +sub help { + print STDOUT <<EOT; +Simple FDR random permutation peak caller +Usage: [options] [files in bedgraph or GFF format] + +Options: +EOT + + my $opt_len = 0; + foreach (keys %vars) { + my $l = length($_); + $opt_len = $l if $l > $opt_len; + } + + $opt_len+=2; + + my $cols= `tput cols` || 80; + + my ($v, $val, $def, $def_format); + my $help_format = "format STDOUT =\n" + .' '.'^'.'<'x$opt_len . ' '. '^' . '<'x($cols-$opt_len-4) . "\n" + .'$v, $def_format'."\n" + .' '.'^'.'<'x$opt_len . ' '. '^' . '<'x($cols-$opt_len-6) . "~~\n" + .'$v, $def_format'."\n" + .".\n"; + + eval $help_format; + die $@ if $@; + + foreach my $k (sort (keys %vars)) { + ($v, $val, $def) = ($k, $vars{$k}, $vars_details{$k}); + $def||=""; + $def_format = $val ? "$def\n\r[Current value: $val]" : $def; + $v = "--$v"; +# format = +# ^<<<<<<<<<<<<<<<<<<<< ^<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< +#$v, $def_format +# ^<<<<<<<<<<<<<<<<<<<< ^<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< ~~ +#$v, $def_format +#. + + write(); + + } + print STDOUT "\n"; + exit 1; +}
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/hp1.bed Fri Apr 20 04:34:17 2018 -0400 @@ -0,0 +1,2000 @@ +3L 17415961 17416091 -2.4153674931314 +3L 17416091 17416196 -2.55492586585599 +3L 17416196 17416656 2.05932465978737 +3L 17416656 17417546 -0.599249371648761 +3L 17417546 17418033 -4.07452348956241 +3L 17418033 17418337 -6.20594954430726 +3L 17418337 17418447 -5.26249039904586 +3L 17418447 17418719 -2.77645084089879 +3L 17418719 17420249 -1.51483717395 +3L 17420249 17420459 -3.807732051728 +3L 17420459 17420660 -4.17530933603838 +3L 17420660 17421213 -4.17530933603838 +3L 17421213 17421223 -2.90724258636772 +3L 17421223 17421228 -3.42663311503462 +3L 17421228 17421615 -3.56515643745284 +3L 17421615 17421629 0 +3L 17421629 17421762 -1.05910489381442 +3L 17421762 17421772 -2.0877734257891 +3L 17421772 17421946 -2.0877734257891 +3L 17421946 17421980 -1.39256885351249 +3L 17421980 17422107 -1.39256885351249 +3L 17422107 17422172 -2.90724258636772 +3L 17422172 17422355 -3.05563971869915 +3L 17422355 17422490 -2.0877734257891 +3L 17422490 17422792 -1.20186343010074 +3L 17422792 17423430 -2.84332836994737 +3L 17423430 17423460 -2.0877734257891 +3L 17423460 17423535 -1.39256885351249 +3L 17423535 17423600 0 +3L 17423600 17423686 0 +3L 17423686 17423739 -1.39256885351249 +3L 17423739 17423973 -2.34017592319915 +3L 17423973 17424090 -2.4153674931314 +3L 17424090 17424215 -3.42663311503462 +3L 17424215 17424809 -4.35793833134167 +3L 17424809 17424921 -3.56515643745284 +3L 17424921 17425440 -3.42663311503462 +3L 17425440 17425825 -2.26084893649928 +3L 17425825 17425970 -2.0877734257891 +3L 17425970 17426273 -3.42663311503462 +3L 17426273 17426428 -0.575168107609865 +3L 17426428 17426746 0.00482907107359622 +3L 17426746 17426777 -0.00467530052002015 +3L 17426777 17426885 -2.55492586585599 +3L 17426885 17427159 -2.77645084089879 +3L 17427159 17427895 -3.38234266122887 +3L 17427895 17427930 -2.55492586585599 +3L 17427930 17427961 0 +3L 17427961 17427970 0 +3L 17427970 17428057 0 +3L 17428057 17428207 0 +3L 17428207 17428349 -2.0877734257891 +3L 17428349 17428494 1.19717008618447 +3L 17428494 17428510 2.97408106030375 +3L 17428510 17428639 5.06185448609285 +3L 17428639 17428708 4.10441267255218 +3L 17428708 17428781 0 +3L 17428781 17428864 0 +3L 17428864 17428938 0 +3L 17428938 17429057 0 +3L 17429057 17429071 0 +3L 17429071 17430048 0 +3L 17430048 17430148 0 +3L 17430148 17430342 -1.78164771235138 +3L 17430342 17430402 -2.55492586585599 +3L 17430402 17431777 -2.53376819312325 +3L 17431777 17431963 -3.72147453126565 +3L 17431963 17432187 -5.53955621377484 +3L 17432187 17432272 -4.57024464575376 +3L 17432272 17432335 0 +3L 17432335 17432373 0 +3L 17432373 17432404 0 +3L 17432404 17432455 0 +3L 17432455 17432587 0 +3L 17432587 17432785 0 +3L 17432785 17432846 0 +3L 17432846 17432868 -1.39256885351249 +3L 17432868 17432876 -2.0877734257891 +3L 17432876 17433362 -1.73208775926791 +3L 17433362 17434012 -1.56098342977366 +3L 17434012 17434745 -2.0877734257891 +3L 17434745 17434875 -1.39256885351249 +3L 17434875 17435361 2.35624268381577 +3L 17435361 17435492 -1.05910489381442 +3L 17435492 17435510 0 +3L 17435510 17435577 0 +3L 17435577 17435637 -1.39256885351249 +3L 17435637 17436191 1.28498697177989 +3L 17436191 17436240 1.54948680669619 +3L 17436240 17436340 -4.10890705137129 +3L 17436340 17436797 -5.06648320422561 +3L 17436797 17436977 -3.1901861920131 +3L 17436977 17437344 -2.0877734257891 +3L 17437344 17437377 -2.0877734257891 +3L 17437377 17437461 -2.0877734257891 +3L 17437461 17437783 -2.0877734257891 +3L 17437783 17443646 -0.318351482828144 +3L 17443646 17443808 -0.858163070035469 +3L 17443808 17444686 -0.345861399909242 +3L 17444686 17444691 -2.0877734257891 +3L 17444691 17444866 -2.74181239868136 +3L 17444866 17445180 -4.21373152611049 +3L 17445180 17445580 -5.74298869455855 +3L 17445580 17445763 -4.2995935275974 +3L 17445763 17446580 -1.05910489381442 +3L 17446580 17446681 0 +3L 17446681 17446688 0 +3L 17446688 17446769 0 +3L 17446769 17446858 0 +3L 17446858 17447621 0 +3L 17447621 17448967 -0.424601364136886 +3L 17448967 17449441 -3.25302896243619 +3L 17449441 17450191 -2.18474136668502 +3L 17450191 17450215 -3.42663311503462 +3L 17450215 17450221 -3.42663311503462 +3L 17450221 17450517 -3.02715636847465 +3L 17450517 17450614 -2.0877734257891 +3L 17450614 17450779 -0.858163070035469 +3L 17450779 17451018 -2.0877734257891 +3L 17451018 17451388 -3.91526259934797 +3L 17451388 17451450 -4.10890705137129 +3L 17451450 17451785 0 +3L 17451785 17451802 0 +3L 17451802 17451823 0 +3L 17451823 17452104 0 +3L 17452104 17452122 0 +3L 17452122 17452220 0 +3L 17452220 17452312 0 +3L 17452312 17452430 -4.43223310796149 +3L 17452430 17452793 -6.04479410220437 +3L 17452793 17453018 -5.03106084870693 +3L 17453018 17453081 -3.42663311503462 +3L 17453081 17453141 0 +3L 17453141 17453572 0 +3L 17453572 17453728 0 +3L 17453728 17453785 -5.20006743060371 +3L 17453785 17454084 -6.15420576828151 +3L 17454084 17454230 -4.57024464575376 +3L 17454230 17454392 -2.68216688920208 +3L 17454392 17455057 1.47052080903416 +3L 17455057 17455156 -2.55492586585599 +3L 17455156 17455177 -2.90724258636772 +3L 17455177 17455349 -2.90724258636772 +3L 17455349 17458191 1.31727823744687 +3L 17458191 17458535 -6.50687455466741 +3L 17458535 17459678 -3.48014999994755 +3L 17459678 17459795 -3.10189685114339 +3L 17459795 17459885 -2.0877734257891 +3L 17459885 17459891 0 +3L 17459891 17462120 0 +3L 17462120 17462127 0 +3L 17462127 17462135 0 +3L 17462135 17463404 0 +3L 17463404 17463583 1.14542132688302 +3L 17463583 17464217 -0.445953270420312 +3L 17464217 17465445 -2.27041110948774 +3L 17465445 17466650 -1.6339116082142 +3L 17466650 17466680 0 +3L 17466680 17466760 0 +3L 17466760 17466884 0 +3L 17466884 17466977 0 +3L 17466977 17467092 -1.39256885351249 +3L 17467092 17467211 -4.43223310796149 +3L 17467211 17467224 -5.81418010758948 +3L 17467224 17467427 -5.56454704262344 +3L 17467427 17467484 -4.75527859770636 +3L 17467484 17468158 -3.42663311503462 +3L 17468158 17468165 -2.0877734257891 +3L 17468165 17468184 -0.561313377016609 +3L 17468184 17468225 -0.477701727681987 +3L 17468225 17468235 -0.477701727681987 +3L 17468235 17468843 -1.01864478926391 +3L 17468843 17469130 -0.768562895851621 +3L 17469130 17469147 -5.32232407845369 +3L 17469147 17469263 -0.454542538543331 +3L 17469263 17469274 2.15461189972513 +3L 17469274 17469286 2.15461189972513 +3L 17469286 17469357 1.63522137105823 +3L 17469357 17469804 0.386152176625108 +3L 17469804 17469992 -3.42663311503462 +3L 17469992 17470111 -5.06648320422561 +3L 17470111 17470411 -5.63702259370271 +3L 17470411 17470494 -5.06648320422561 +3L 17470494 17470551 -4.46799333170228 +3L 17470551 17470740 -5.13482089478576 +3L 17470740 17470749 -5.32232407845369 +3L 17470749 17471171 -4.01533135855775 +3L 17471171 17471545 -5.9719432181294 +3L 17471545 17472158 -4.50976745435994 +3L 17472158 17472208 -4.57024464575376 +3L 17472208 17472528 -3.62972994085222 +3L 17472528 17472911 -4.57024464575376 +3L 17472911 17473380 -1.51483717395 +3L 17473380 17473972 -3.2733826212464 +3L 17473972 17474085 -2.90724258636772 +3L 17474085 17474212 0 +3L 17474212 17474283 2.71184381903969 +3L 17474283 17474436 2.64648699296144 +3L 17474436 17474901 -1.65438943074425 +3L 17474901 17475220 -4.78959902445488 +3L 17475220 17475340 -3.10189685114339 +3L 17475340 17475435 -1.39256885351249 +3L 17475435 17475531 -2.0877734257891 +3L 17475531 17475671 -2.90724258636772 +3L 17475671 17475678 -0.735183064892071 +3L 17475678 17476292 -0.576591791462247 +3L 17476292 17476480 0.0660948934056327 +3L 17476480 17476781 -0.0323534318009542 +3L 17476781 17477316 -2.8952683936586 +3L 17477316 17477588 -5.27465730379902 +3L 17477588 17477672 -4.35793833134167 +3L 17477672 17477826 -2.90724258636772 +3L 17477826 17477948 -1.66324290730034 +3L 17477948 17478026 0 +3L 17478026 17478144 0 +3L 17478144 17478181 -3.42663311503462 +3L 17478181 17478185 -4.35793833134167 +3L 17478185 17478499 -4.5302095028321 +3L 17478499 17478658 -3.10189685114339 +3L 17478658 17478792 -1.05910489381442 +3L 17478792 17478922 1.47301539422532 +3L 17478922 17479040 4.2397784396702 +3L 17479040 17479161 3.35792826791411 +3L 17479161 17479470 1.87731795947831 +3L 17479470 17479738 1.29186720687598 +3L 17479738 17479813 1.22281083040574 +3L 17479813 17479896 0 +3L 17479896 17479966 0 +3L 17479966 17480025 0 +3L 17480025 17480327 -1.72053005185583 +3L 17480327 17480437 -0.447537621337495 +3L 17480437 17480608 5.06185448609285 +3L 17480608 17480672 4.10441267255218 +3L 17480672 17481129 -1.73208775926791 +3L 17481129 17481261 -4.35793833134167 +3L 17481261 17481637 -6.35281248517157 +3L 17481637 17482130 -0.418070913410033 +3L 17482130 17484190 -0.463127980294898 +3L 17484190 17484263 1.54948680669619 +3L 17484263 17484320 0.593861154390571 +3L 17484320 17485061 -0.78488328593786 +3L 17485061 17485093 -5.20006743060371 +3L 17485093 17485171 -4.57024464575376 +3L 17485171 17486914 -1.66324290730034 +3L 17486914 17486923 -2.0877734257891 +3L 17486923 17487308 -1.58598905941738 +3L 17487308 17487901 -0.0037354757900331 +3L 17487901 17488016 -4.10890705137129 +3L 17488016 17488613 -5.00289597475257 +3L 17488613 17489570 -3.00149997803535 +3L 17489570 17490118 -1.56247537501058 +3L 17490118 17490580 -5.16314576788544 +3L 17490580 17491778 -1.87861712142024 +3L 17491778 17493333 0.759971479660135 +3L 17493333 17493632 2.10425581506495 +3L 17493632 17493999 -2.37330132707517 +3L 17493999 17494030 -6.18030962448475 +3L 17494030 17494714 -3.73914355418382 +3L 17494714 17494916 -2.0877734257891 +3L 17494916 17495035 -2.0877734257891 +3L 17495035 17495050 -2.90724258636772 +3L 17495050 17495349 -2.90724258636772 +3L 17495349 17495368 -2.90724258636772 +3L 17495368 17496292 -1.81084455985473 +3L 17496292 17496815 -3.06108643753096 +3L 17496815 17497107 -7.60551312373146 +3L 17497107 17497281 -6.39057025028289 +3L 17497281 17498771 -2.63260885305845 +3L 17498771 17498843 -2.0877734257891 +3L 17498843 17499122 -1.84829023815112 +3L 17499122 17500261 -1.69376753678793 +3L 17500261 17500303 -3.1901861920131 +3L 17500303 17501284 -4.10890705137129 +3L 17501284 17502403 -2.45154299114717 +3L 17502403 17502668 0 +3L 17502668 17503445 -2.90724258636772 +3L 17503445 17504943 -0.38913829674294 +3L 17504943 17505702 -3.06840310089978 +3L 17505702 17507677 -0.67989333316847 +3L 17507677 17507808 0 +3L 17507808 17507832 0 +3L 17507832 17507852 -4.23878936063925 +3L 17507852 17507856 -5.20006743060371 +3L 17507856 17508156 -5.53955621377484 +3L 17508156 17508488 -5.98681066577781 +3L 17508488 17508515 -6.36234497216546 +3L 17508515 17508887 -6.0209168975093 +3L 17508887 17509040 -3.69153660144476 +3L 17509040 17509128 -4.19678119846456 +3L 17509128 17509386 -5.75026801676681 +3L 17509386 17509604 -0.989370954409652 +3L 17509604 17509706 2.15461189972513 +3L 17509706 17509742 1.19717008618447 +3L 17509742 17509902 -4.27960886543453 +3L 17509902 17510154 -5.48823712025589 +3L 17510154 17510453 -4.78959902445488 +3L 17510453 17510953 -4.38624911229275 +3L 17510953 17511001 -3.1901861920131 +3L 17511001 17511294 0 +3L 17511294 17511451 -1.39256885351249 +3L 17511451 17512057 -3.47430077637544 +3L 17512057 17512265 -4.50288856075644 +3L 17512265 17512513 -3.42663311503462 +3L 17512513 17512639 -3.62972994085222 +3L 17512639 17513119 -0.774562500606804 +3L 17513119 17513206 -0.465831973201576 +3L 17513206 17513372 -2.0877734257891 +3L 17513372 17513473 -1.39256885351249 +3L 17513473 17513724 0.597550404454718 +3L 17513724 17513978 0.612784511019419 +3L 17513978 17514100 0.605070183043923 +3L 17514100 17514207 -0.993730880144887 +3L 17514207 17514502 -1.02768019764921 +3L 17514502 17514646 -1.10799571635366 +3L 17514646 17514775 -1.66324290730034 +3L 17514775 17514957 0 +3L 17514957 17514981 0 +3L 17514981 17515158 0 +3L 17515158 17515438 0 +3L 17515438 17515480 0 +3L 17515480 17515513 0 +3L 17515513 17515697 -2.68216688920208 +3L 17515697 17515933 -4.79805306077022 +3L 17515933 17516353 -5.58213625649124 +3L 17516353 17516774 -5.43502497407227 +3L 17516774 17516853 -2.0877734257891 +3L 17516853 17516859 -2.90724258636772 +3L 17516859 17517387 -2.0877734257891 +3L 17517387 17518193 0 +3L 17518193 17518719 0 +3L 17518719 17518938 0 +3L 17518938 17518968 0 +3L 17518968 17520878 0 +3L 17520878 17521329 0 +3L 17521329 17522145 0 +3L 17522145 17522303 0 +3L 17522303 17522312 0 +3L 17522312 17522535 0 +3L 17522535 17522608 0 +3L 17522608 17522657 0 +3L 17522657 17522669 0 +3L 17522669 17522852 0 +3L 17522852 17524407 0 +3L 17524407 17524567 -2.0877734257891 +3L 17524567 17524770 -3.42663311503462 +3L 17524770 17524930 -2.0877734257891 +3L 17524930 17525095 -2.90724258636772 +3L 17525095 17525736 -1.68371795492887 +3L 17525736 17525756 -1.41022228832274 +3L 17525756 17525819 -0.946973963512023 +3L 17525819 17526336 -1.46591272074667 +3L 17526336 17526821 -6.93872222643839 +3L 17526821 17527360 -1.19901669582865 +3L 17527360 17527404 -5.20006743060371 +3L 17527404 17527807 -3.42663311503462 +3L 17527807 17527835 -2.90724258636772 +3L 17527835 17527967 -2.90724258636772 +3L 17527967 17528159 -3.2733826212464 +3L 17528159 17528167 -3.42663311503462 +3L 17528167 17529070 -3.70081461656489 +3L 17529070 17529247 -2.90724258636772 +3L 17529247 17529271 -2.0877734257891 +3L 17529271 17529362 -1.39256885351249 +3L 17529362 17529386 0 +3L 17529386 17529443 0 +3L 17529443 17529454 -1.39256885351249 +3L 17529454 17529515 -2.0877734257891 +3L 17529515 17531020 -1.51483717395 +3L 17531020 17531727 0.97620445390234 +3L 17531727 17533027 0.242037294667923 +3L 17533027 17534833 -1.28119532926073 +3L 17534833 17536554 -1.05910489381442 +3L 17536554 17537886 -3.09228099773779 +3L 17537886 17539148 -3.37733650571732 +3L 17539148 17539876 -1.2478082245722 +3L 17539876 17540663 -4.27227252726769 +3L 17540663 17541422 -0.297000344899299 +3L 17541422 17541576 -3.69153660144476 +3L 17541576 17541767 -3.53175247535673 +3L 17541767 17541890 -2.55492586585599 +3L 17541890 17542025 0 +3L 17542025 17542096 -1.39256885351249 +3L 17542096 17542357 -2.0877734257891 +3L 17542357 17542363 -2.0877734257891 +3L 17542363 17542427 -2.0877734257891 +3L 17542427 17542901 -1.92086439808696 +3L 17542901 17543489 -5.03503988948319 +3L 17543489 17543518 -5.53955621377484 +3L 17543518 17543670 -4.01533135855775 +3L 17543670 17543750 0 +3L 17543750 17543843 0 +3L 17543843 17543996 0 +3L 17543996 17544084 -2.4153674931314 +3L 17544084 17544844 -4.39893625820657 +3L 17544844 17545104 -1.05285097901053 +3L 17545104 17545566 -0.210459594582568 +3L 17545566 17545727 -2.55492586585599 +3L 17545727 17545821 -2.55492586585599 +3L 17545821 17545955 -2.90724258636772 +3L 17545955 17546526 -3.62972994085222 +3L 17546526 17546663 -4.69619873626645 +3L 17546663 17546890 -5.77188817111398 +3L 17546890 17546957 -4.83959573744425 +3L 17546957 17547289 -2.29310500926643 +3L 17547289 17547433 -2.4153674931314 +3L 17547433 17547495 -1.39256885351249 +3L 17547495 17548267 -1.30899491083192 +3L 17548267 17548293 0 +3L 17548293 17548748 -1.73208775926791 +3L 17548748 17548916 -1.05910489381442 +3L 17548916 17549057 1.89754591271408 +3L 17549057 17549666 0.4779733961301 +3L 17549666 17549697 -4.57024464575376 +3L 17549697 17549872 -3.807732051728 +3L 17549872 17550014 -2.0877734257891 +3L 17550014 17550128 0 +3L 17550128 17550234 -1.05910489381442 +3L 17550234 17550437 -2.0877734257891 +3L 17550437 17550527 -1.05910489381442 +3L 17550527 17550532 0 +3L 17550532 17550614 0 +3L 17550614 17550680 -1.39256885351249 +3L 17550680 17551338 -3.03988511098966 +3L 17551338 17551598 -5.43502497407227 +3L 17551598 17551921 -5.25022001190748 +3L 17551921 17551957 -4.35793833134167 +3L 17551957 17552635 -3.84084102162073 +3L 17552635 17553047 -4.01533135855775 +3L 17553047 17553153 -2.0877734257891 +3L 17553153 17553261 0 +3L 17553261 17553380 -1.05910489381442 +3L 17553380 17553737 -2.47282073788985 +3L 17553737 17553906 -3.1901861920131 +3L 17553906 17554009 -2.90724258636772 +3L 17554009 17554447 -4.06954396796188 +3L 17554447 17554877 -4.01526047650326 +3L 17554877 17555107 -2.25406869409716 +3L 17555107 17555137 -2.25032852800284 +3L 17555137 17555149 -2.25032852800284 +3L 17555149 17555225 -2.24802544331815 +3L 17555225 17555288 -2.0877734257891 +3L 17555288 17555297 -2.0877734257891 +3L 17555297 17555438 -4.27960886543453 +3L 17555438 17555694 -5.77188817111398 +3L 17555694 17556220 -5.66039399675287 +3L 17556220 17556327 -3.69153660144476 +3L 17556327 17556409 0 +3L 17556409 17556539 -1.05910489381442 +3L 17556539 17556923 -1.89106755736324 +3L 17556923 17556942 0 +3L 17556942 17557111 -5.5055471133017 +3L 17557111 17557358 -7.06918138875176 +3L 17557358 17557793 -4.81203395618529 +3L 17557793 17558334 -0.784441637777609 +3L 17558334 17558933 -1.48854786879403 +3L 17558933 17559785 -4.27960886543453 +3L 17559785 17559822 -4.91925594869484 +3L 17559822 17559911 -3.56515643745284 +3L 17559911 17560555 -2.97508771941039 +3L 17560555 17560651 -4.10890705137129 +3L 17560651 17561015 -5.32232407845369 +3L 17561015 17561024 -5.32232407845369 +3L 17561024 17561169 -3.1901861920131 +3L 17561169 17561434 -5.32232407845369 +3L 17561434 17561705 -5.32232407845369 +3L 17561705 17561716 -5.53955621377484 +3L 17561716 17561723 -5.53955621377484 +3L 17561723 17562025 -4.68408678906606 +3L 17562025 17562069 -4.35793833134167 +3L 17562069 17562172 -3.10189685114339 +3L 17562172 17562216 -2.0877734257891 +3L 17562216 17563368 -1.95175728612821 +3L 17563368 17563375 -2.90724258636772 +3L 17563375 17564258 -3.24829090327753 +3L 17564258 17565016 -1.30899491083192 +3L 17565016 17565308 -1.84829023815112 +3L 17565308 17567904 -0.128139338442937 +3L 17567904 17569246 -2.61659853787386 +3L 17569246 17570037 -0.452591759356602 +3L 17570037 17570580 0.215290702927938 +3L 17570580 17574745 -0.301246061756575 +3L 17574745 17574829 0 +3L 17574829 17574992 -1.05910489381442 +3L 17574992 17575203 -0.606784390052116 +3L 17575203 17575481 -0.483868520600485 +3L 17575481 17575535 -0.666464442155458 +3L 17575535 17576070 -3.84900088806128 +3L 17576070 17577396 -1.68899058844381 +3L 17577396 17577473 -0.288428787059227 +3L 17577473 17577700 -0.101270964454808 +3L 17577700 17578751 1.09016723029279 +3L 17578751 17579362 -5.4824206591574 +3L 17579362 17579399 -4.10890705137129 +3L 17579399 17580033 -3.32629566508554 +3L 17580033 17580177 -1.05910489381442 +3L 17580177 17580430 0 +3L 17580430 17581273 -2.26084893649928 +3L 17581273 17582441 -2.02136769931885 +3L 17582441 17582792 -3.42663311503462 +3L 17582792 17582989 -5.81418010758948 +3L 17582989 17583125 -4.27960886543453 +3L 17583125 17583265 0 +3L 17583265 17583283 0 +3L 17583283 17583303 0.138248241579792 +3L 17583303 17583308 0.142598537398012 +3L 17583308 17583313 0.142598537398012 +3L 17583313 17583318 0.142598537398012 +3L 17583318 17583758 -0.624577987581715 +3L 17583758 17583900 -3.42663311503462 +3L 17583900 17585543 0.281560551003223 +3L 17585543 17585904 -0.285607350758448 +3L 17585904 17586249 -5.58911233312691 +3L 17586249 17586266 -5.20006743060371 +3L 17586266 17586766 -4.81897392259593 +3L 17586766 17586803 -6.92224372353247 +3L 17586803 17587287 -2.35013771517629 +3L 17587287 17587448 -0.648089094246175 +3L 17587448 17587653 -3.35071059435935 +3L 17587653 17587716 -6.11414243765132 +3L 17587716 17588129 -1.01623854719582 +3L 17588129 17589090 -1.03788438086145 +3L 17589090 17589600 -0.841423097807887 +3L 17589600 17589629 -2.0877734257891 +3L 17589629 17589690 0.474682731699959 +3L 17589690 17590148 -0.587120068470401 +3L 17590148 17590242 -2.90724258636772 +3L 17590242 17590296 -4.83959573744425 +3L 17590296 17590527 -5.91478453479568 +3L 17590527 17591027 -4.2995935275974 +3L 17591027 17591188 -0.701739429236591 +3L 17591188 17591498 -0.71853044186611 +3L 17591498 17592725 -0.813683386614058 +3L 17592725 17592975 -4.35793833134167 +3L 17592975 17593126 -3.42663311503462 +3L 17593126 17593148 0 +3L 17593148 17594218 -2.68216688920208 +3L 17594218 17594513 -6.37370143374959 +3L 17594513 17594948 -4.45277592863354 +3L 17594948 17595187 -5.13482089478576 +3L 17595187 17595214 -5.20006743060371 +3L 17595214 17595361 -4.01533135855775 +3L 17595361 17595423 -2.55492586585599 +3L 17595423 17597356 -0.164019374992697 +3L 17597356 17597717 0 +3L 17597717 17597893 -1.05910489381442 +3L 17597893 17598236 -1.89106755736324 +3L 17598236 17598366 0 +3L 17598366 17598374 0 +3L 17598374 17598425 -4.35793833134167 +3L 17598425 17599160 -3.64685059781814 +3L 17599160 17599290 -3.42663311503462 +3L 17599290 17599797 -0.136652434608378 +3L 17599797 17599802 0.138248241579792 +3L 17599802 17600074 -3.31324822703235 +3L 17600074 17600811 -3.50383471530534 +3L 17600811 17600978 -1.05910489381442 +3L 17600978 17602232 0 +3L 17602232 17602246 0 +3L 17602246 17602312 -2.55492586585599 +3L 17602312 17602707 -4.57024464575376 +3L 17602707 17603585 -4.00038075543875 +3L 17603585 17603694 -3.62972994085222 +3L 17603694 17604342 -2.29310500926643 +3L 17604342 17604789 -3.42663311503462 +3L 17604789 17605967 -2.57821119676241 +3L 17605967 17606105 0.878621930812341 +3L 17606105 17607076 0.285527518055443 +3L 17607076 17608158 -1.07009921718573 +3L 17608158 17608541 -1.60165217054032 +3L 17608541 17608967 -1.60424617912047 +3L 17608967 17609177 -1.09695858295787 +3L 17609177 17609196 -0.842950586461244 +3L 17609196 17609607 -1.86038923743962 +3L 17609607 17609624 -6.18030962448475 +3L 17609624 17610218 0.139956089353332 +3L 17610218 17610426 -2.74181239868136 +3L 17610426 17610811 -3.35204179418636 +3L 17610811 17612576 -1.42786134193199 +3L 17612576 17613187 -1.32508429577872 +3L 17613187 17613423 -1.79798622879982 +3L 17613423 17614645 -0.338520660620992 +3L 17614645 17616072 -2.99809932148944 +3L 17616072 17616535 -2.55492586585599 +3L 17616535 17617137 -3.1613526114503 +3L 17617137 17617545 -2.21504823233841 +3L 17617545 17617705 0.782245620658319 +3L 17617705 17618337 -1.94084088676422 +3L 17618337 17618418 -3.62972994085222 +3L 17618418 17618605 -6.38122304655397 +3L 17618605 17619027 -8.10087927855092 +3L 17619027 17619885 -4.13138244802135 +3L 17619885 17620733 -2.80770304525312 +3L 17620733 17620748 -3.807732051728 +3L 17620748 17621218 -3.807732051728 +3L 17621218 17623011 -0.157223115524371 +3L 17623011 17623380 -1.27170196313972 +3L 17623380 17623588 -2.18514083349501 +3L 17623588 17624133 -0.294812617799881 +3L 17624133 17624767 -2.42766436630248 +3L 17624767 17625011 -6.91448886541995 +3L 17625011 17625904 -1.50243303568788 +3L 17625904 17626130 -1.81673475484425 +3L 17626130 17626196 -1.81673475484425 +3L 17626196 17626295 -1.30049048607261 +3L 17626295 17626600 0.914699491996825 +3L 17626600 17626951 -0.00462871813275598 +3L 17626951 17626968 -2.0877734257891 +3L 17626968 17627060 -2.0877734257891 +3L 17627060 17627492 -3.54616326734426 +3L 17627492 17628436 -3.24829090327753 +3L 17628436 17628481 -1.39256885351249 +3L 17628481 17628611 0 +3L 17628611 17628715 0 +3L 17628715 17628763 -2.55492586585599 +3L 17628763 17628774 -3.42663311503462 +3L 17628774 17628949 -3.56515643745284 +3L 17628949 17629032 0.474682731699959 +3L 17629032 17629240 1.43212454524063 +3L 17629240 17629676 -2.04892675949615 +3L 17629676 17629865 -3.10189685114339 +3L 17629865 17630617 -4.69070590795494 +3L 17630617 17632169 0.371004838048949 +3L 17632169 17633308 -4.05701944763538 +3L 17633308 17633492 -2.74181239868136 +3L 17633492 17633532 -3.807732051728 +3L 17633532 17634514 -2.54509341053212 +3L 17634514 17635093 0.0380992621605252 +3L 17635093 17635776 -1.50720774125216 +3L 17635776 17636157 -4.23878936063925 +3L 17636157 17636741 -4.83046685785534 +3L 17636741 17636753 -4.91925594869484 +3L 17636753 17637011 -3.93583934205124 +3L 17637011 17637090 -3.1901861920131 +3L 17637090 17637236 -2.0877734257891 +3L 17637236 17637443 -2.90724258636772 +3L 17637443 17637540 -1.66324290730034 +3L 17637540 17637719 3.56078882001442 +3L 17637719 17637772 5.06185448609285 +3L 17637772 17637966 -0.13074778143034 +3L 17637966 17638396 -3.6405238550957 +3L 17638396 17638405 -1.93572795170745 +3L 17638405 17638934 -0.435118901742244 +3L 17638934 17639067 0.140389670025455 +3L 17639067 17639208 -2.90724258636772 +3L 17639208 17639274 1.54948680669619 +3L 17639274 17639844 1.15643888068263 +3L 17639844 17639865 -2.90724258636772 +3L 17639865 17639883 0 +3L 17639883 17640066 0 +3L 17640066 17640073 0 +3L 17640073 17641106 -1.85135270932415 +3L 17641106 17641395 -0.389717565869056 +3L 17641395 17641451 0.369856894136833 +3L 17641451 17641558 0.296680620824183 +3L 17641558 17641818 -4.82312193437512 +3L 17641818 17641823 -5.06648320422561 +3L 17641823 17641869 -5.06648320422561 +3L 17641869 17642588 -4.35793833134167 +3L 17642588 17642956 -5.13482089478576 +3L 17642956 17643012 -3.96616443097239 +3L 17643012 17644161 -3.807732051728 +3L 17644161 17645005 -0.499882197163522 +3L 17645005 17645127 0.134155704979806 +3L 17645127 17645324 -4.35793833134167 +3L 17645324 17645397 0 +3L 17645397 17645510 0.458891968871037 +3L 17645510 17645833 0.562393588032954 +3L 17645833 17646037 0.380371358807581 +3L 17646037 17646406 0.401664203993259 +3L 17646406 17647830 -0.726151156521355 +3L 17647830 17648359 -4.54535317265227 +3L 17648359 17648494 -1.66324290730034 +3L 17648494 17648535 0 +3L 17648535 17648556 0 +3L 17648556 17648613 0 +3L 17648613 17648709 0 +3L 17648709 17648801 -2.55492586585599 +3L 17648801 17648845 -3.42663311503462 +3L 17648845 17649078 -4.48903158717769 +3L 17649078 17649089 -4.35793833134167 +3L 17649089 17649209 -0.363580659150097 +3L 17649209 17649420 -0.477701727681987 +3L 17649420 17649725 -3.74762766109031 +3L 17649725 17650288 -5.9719432181294 +3L 17650288 17650357 -4.91925594869484 +3L 17650357 17650417 0 +3L 17650417 17650518 0.878621930812341 +3L 17650518 17650537 0.952947434721555 +3L 17650537 17651040 -0.256258206449979 +3L 17651040 17651082 -4.35793833134167 +3L 17651082 17651117 -3.42663311503462 +3L 17651117 17651244 0 +3L 17651244 17655021 0.947030237728085 +3L 17655021 17655247 -6.9975824378003 +3L 17655247 17655378 -2.04447102960565 +3L 17655378 17655400 0.952947434721555 +3L 17655400 17655591 0.491609840339092 +3L 17655591 17655794 -1.43283672364324 +3L 17655794 17655832 1.16435396148042 +3L 17655832 17656189 1.84329380023654 +3L 17656189 17656200 1.46985051136244 +3L 17656200 17656285 0 +3L 17656285 17656357 4.10441267255218 +3L 17656357 17656401 5.06185448609285 +3L 17656401 17656506 -0.246943231713576 +3L 17656506 17656597 -0.621538328926804 +3L 17656597 17657304 -1.38801684858639 +3L 17657304 17657938 0.124815224248221 +3L 17657938 17657947 1.18270602828767 +3L 17657947 17658425 -0.0625631866871179 +3L 17658425 17658831 -1.09961511870992 +3L 17658831 17659195 -2.10395238631383 +3L 17659195 17659257 -2.3938408066937 +3L 17659257 17659358 -0.763461267090944 +3L 17659358 17659366 -0.746685502460338 +3L 17659366 17659687 -0.70422180099761 +3L 17659687 17659791 -0.479488176977641 +3L 17659791 17660124 -0.389214158926654 +3L 17660124 17660569 -0.304799229560287 +3L 17660569 17660897 1.00036031747807 +3L 17660897 17661115 0.390135544601459 +3L 17661115 17661230 0.34641397688408 +3L 17661230 17662153 -1.29181411672887 +3L 17662153 17662265 -1.39256885351249 +3L 17662265 17662272 0 +3L 17662272 17662331 -1.39256885351249 +3L 17662331 17662337 -2.0877734257891 +3L 17662337 17662353 -4.46799333170228 +3L 17662353 17662790 -5.43502497407227 +3L 17662790 17662867 -1.57898014246747 +3L 17662867 17663294 -1.7808984540691 +3L 17663294 17663594 0.0240014907903035 +3L 17663594 17663729 0.782245620658319 +3L 17663729 17663765 2.15461189972513 +3L 17663765 17664031 1.32271093190903 +3L 17664031 17664260 -1.80880845102141 +3L 17664260 17664483 -3.05563971869915 +3L 17664483 17664537 0 +3L 17664537 17664588 0 +3L 17664588 17664661 0 +3L 17664661 17666051 0 +3L 17666051 17668139 0 +3L 17668139 17673548 -0.341642528914632 +3L 17673548 17674019 -1.51483717395 +3L 17674019 17676972 0 +3L 17676972 17677537 -3.32629566508554 +3L 17677537 17677973 -2.99390012510116 +3L 17677973 17678133 -2.90724258636772 +3L 17678133 17678179 -3.42663311503462 +3L 17678179 17678559 -4.0628777571469 +3L 17678559 17680923 -2.52976635235712 +3L 17680923 17681135 0 +3L 17681135 17681586 -0.118002065303967 +3L 17681586 17681806 -0.169762019400443 +3L 17681806 17682078 -0.608224247767604 +3L 17682078 17682652 -2.65139317599471 +3L 17682652 17682698 -2.0877734257891 +3L 17682698 17682852 -2.4153674931314 +3L 17682852 17683002 -2.4153674931314 +3L 17683002 17683077 -3.1901861920131 +3L 17683077 17683287 -5.20006743060371 +3L 17683287 17683629 -6.66004130597075 +3L 17683629 17684625 -4.18455072110617 +3L 17684625 17685243 -2.68216688920208 +3L 17685243 17687853 -0.239528056952404 +3L 17687853 17689597 0 +3L 17689597 17692474 -0.502619804033915 +3L 17692474 17692511 -1.39256885351249 +3L 17692511 17692794 -2.0877734257891 +3L 17692794 17692799 -2.0877734257891 +3L 17692799 17693606 0 +3L 17693606 17693649 -1.39256885351249 +3L 17693649 17697988 -0.351064340770172 +3L 17697988 17698787 0 +3L 17698787 17700433 -1.09136340279836 +3L 17700433 17701701 -1.57853798825356 +3L 17701701 17707052 -0.0628581865992481 +3L 17707052 17708279 -0.968018095839878 +3L 17708279 17709771 0 +3L 17709771 17710460 0 +3L 17710460 17712562 -0.642010754900314 +3L 17712562 17712918 -4.50288856075644 +3L 17712918 17714019 -0.722626030916987 +3L 17714019 17714054 0 +3L 17714054 17714263 -0.858163070035469 +3L 17714263 17715005 -1.98381532262656 +3L 17715005 17716159 -0.858163070035469 +3L 17716159 17718347 -0.148431542395496 +3L 17718347 17719504 -0.968018095839878 +3L 17719504 17719606 -4.10890705137129 +3L 17719606 17720224 -4.95329601781634 +3L 17720224 17720503 -3.66713075241527 +3L 17720503 17720905 -2.26084893649928 +3L 17720905 17723215 0 +3L 17723215 17724421 -0.654152189111727 +3L 17724421 17724757 -4.31930513412542 +3L 17724757 17726368 -1.61765775642534 +3L 17726368 17727106 -0.37355103300116 +3L 17727106 17727184 -4.83959573744425 +3L 17727184 17727427 -5.83487006914932 +3L 17727427 17727575 -4.57024464575376 +3L 17727575 17727744 -2.4153674931314 +3L 17727744 17727903 -2.0877734257891 +3L 17727903 17728394 -4.29106254778434 +3L 17728394 17728701 -5.78044632164322 +3L 17728701 17729028 -3.62972994085222 +3L 17729028 17729114 -3.42663311503462 +3L 17729114 17729534 -5.2419813887421 +3L 17729534 17731014 -1.92086439808696 +3L 17731014 17731036 -2.0877734257891 +3L 17731036 17731059 -1.39256885351249 +3L 17731059 17731206 0 +3L 17731206 17731651 0 +3L 17731651 17731754 0 +3L 17731754 17731763 0 +3L 17731763 17731914 0 +3L 17731914 17732357 0 +3L 17732357 17732412 0 +3L 17732412 17732421 0 +3L 17732421 17732634 0 +3L 17732634 17732649 0 +3L 17732649 17733801 0 +3L 17733801 17734448 0 +3L 17734448 17735514 0 +3L 17735514 17735741 0 +3L 17735741 17735966 0 +3L 17735966 17736192 0 +3L 17736192 17736656 0 +3L 17736656 17737687 0 +3L 17737687 17737844 0 +3L 17737844 17737866 0 +3L 17737866 17737889 0 +3L 17737889 17738036 0 +3L 17738036 17738670 0 +3L 17738670 17738838 0 +3L 17738838 17740882 -0.399901407109278 +3L 17740882 17740896 -0.666464442155458 +3L 17740896 17740931 -0.666464442155458 +3L 17740931 17741136 -0.688413530673959 +3L 17741136 17741171 -0.752325621496629 +3L 17741171 17741300 -0.871444287947066 +3L 17741300 17741366 -1.39256885351249 +3L 17741366 17741492 0 +3L 17741492 17741504 0 +3L 17741504 17741522 0 +3L 17741522 17741618 0 +3L 17741618 17741710 1.37031788464809 +3L 17741710 17741911 0.714955746730745 +3L 17741911 17742119 0.417623087339104 +3L 17742119 17742568 -0.994800581944547 +3L 17742568 17743507 -3.45745683701254 +3L 17743507 17744041 -4.26951175376548 +3L 17744041 17744192 0.653546233646707 +3L 17744192 17744329 0.933171563914282 +3L 17744329 17745345 -0.257646747899449 +3L 17745345 17745743 0.806791399922359 +3L 17745743 17745838 1.47301539422532 +3L 17745838 17745953 -2.0877734257891 +3L 17745953 17746036 -2.0877734257891 +3L 17746036 17746100 -2.90724258636772 +3L 17746100 17746231 -3.42663311503462 +3L 17746231 17746612 -3.62972994085222 +3L 17746612 17746637 -2.90724258636772 +3L 17746637 17746757 -3.807732051728 +3L 17746757 17747085 -5.25022001190748 +3L 17747085 17748063 -0.644127553411268 +3L 17748063 17748387 2.91636518721007 +3L 17748387 17748495 0 +3L 17748495 17748545 0 +3L 17748545 17748620 -2.55492586585599 +3L 17748620 17748636 -3.42663311503462 +3L 17748636 17748900 -3.42663311503462 +3L 17748900 17749088 0.806791399922359 +3L 17749088 17749321 0.973611952890587 +3L 17749321 17749436 0.583065402019151 +3L 17749436 17749630 -3.42663311503462 +3L 17749630 17749728 -2.4153674931314 +3L 17749728 17750168 -3.807732051728 +3L 17750168 17750233 -2.90724258636772 +3L 17750233 17750276 -1.39256885351249 +3L 17750276 17750302 -2.0877734257891 +3L 17750302 17750787 -2.98334834764705 +3L 17750787 17751414 -0.705545564574508 +3L 17751414 17751978 1.23457948792951 +3L 17751978 17752206 -1.00059170929627 +3L 17752206 17752460 -1.91511825818251 +3L 17752460 17752476 -5.06648320422561 +3L 17752476 17752495 -5.06648320422561 +3L 17752495 17752699 -5.13482089478576 +3L 17752699 17752910 -4.75527859770636 +3L 17752910 17753205 -5.57933635326641 +3L 17753205 17753457 -3.62972994085222 +3L 17753457 17753481 0 +3L 17753481 17753785 4.75068421171918 +3L 17753785 17753917 4.10441267255218 +3L 17753917 17754058 -1.05910489381442 +3L 17754058 17754185 -2.4153674931314 +3L 17754185 17754347 -2.4153674931314 +3L 17754347 17754479 -2.90724258636772 +3L 17754479 17755293 -3.91526259934797 +3L 17755293 17755368 -2.55492586585599 +3L 17755368 17755541 -2.4153674931314 +3L 17755541 17755703 -2.4153674931314 +3L 17755703 17755807 -0.313705954818164 +3L 17755807 17756029 -0.634645775943271 +3L 17756029 17756044 -1.86038923743962 +3L 17756044 17756308 -1.87002455654501 +3L 17756308 17756508 -3.9177526090449 +3L 17756508 17756773 -7.50323570073859 +3L 17756773 17756873 -5.81418010758948 +3L 17756873 17757959 1.62817589197509 +3L 17757959 17758155 -1.37563464349077 +3L 17758155 17758427 -1.40028315350858 +3L 17758427 17758744 -2.34181989623228 +3L 17758744 17758789 -7.73002496590584 +3L 17758789 17760033 -4.43223310796149 +3L 17760033 17760090 -0.253525658789485 +3L 17760090 17760117 -0.138212944510863 +3L 17760117 17760508 -0.244165820788913 +3L 17760508 17760643 0.140389670025455 +3L 17760643 17760761 -2.0877734257891 +3L 17760761 17760962 -3.42663311503462 +3L 17760962 17761041 -2.55492586585599 +3L 17761041 17761071 0 +3L 17761071 17761210 -2.5533536176369 +3L 17761210 17761227 -2.62329829623993 +3L 17761227 17761546 -2.6414149902627 +3L 17761546 17761683 -2.7630423648933 +3L 17761683 17762044 -3.42663311503462 +3L 17762044 17762095 -2.90724258636772 +3L 17762095 17762209 -3.2733826212464 +3L 17762209 17762992 -0.0727540107117133 +3L 17762992 17763141 -3.42663311503462 +3L 17763141 17764123 -3.0353520127332 +3L 17764123 17764164 -1.39256885351249 +3L 17764164 17764198 -2.0877734257891 +3L 17764198 17764977 -0.728547392956818 +3L 17764977 17765584 -1.48999701854234 +3L 17765584 17765913 -1.34538961827804 +3L 17765913 17766143 -2.6489316875206 +3L 17766143 17766805 -1.39676547041076 +3L 17766805 17766957 0.878621930812341 +3L 17766957 17767696 0 +3L 17767696 17768720 0 +3L 17768720 17769796 0 +3L 17769796 17771844 0 +3L 17771844 17772898 0 +3L 17772898 17773072 0 +3L 17773072 17773978 -1.32217178994419 +3L 17773978 17774506 -1.60018500752036 +3L 17774506 17775399 -0.585097891945795 +3L 17775399 17777046 -0.190733851261496 +3L 17777046 17779034 0 +3L 17779034 17781783 0 +3L 17781783 17784486 -0.525380204801075 +3L 17784486 17784848 -4.50288856075644 +3L 17784848 17785812 -3.39513642834552 +3L 17785812 17786884 -1.45091897165284 +3L 17786884 17787298 -1.89106755736324 +3L 17787298 17789285 -0.660399787137409 +3L 17789285 17789392 0 +3L 17789392 17790888 -0.827028284974454 +3L 17790888 17791625 0 +3L 17791625 17792135 -1.39256885351249 +3L 17792135 17792557 -3.65654337623127 +3L 17792557 17792708 -1.05910489381442 +3L 17792708 17792833 0 +3L 17792833 17792913 -2.0877734257891 +3L 17792913 17793109 -3.42663311503462 +3L 17793109 17793291 -4.01533135855775 +3L 17793291 17793348 -2.90724258636772 +3L 17793348 17793396 0 +3L 17793396 17793540 -4.10890705137129 +3L 17793540 17793859 -5.63702259370271 +3L 17793859 17793949 -4.83959573744425 +3L 17793949 17794147 -3.96616443097239 +3L 17794147 17794390 -5.13482089478576 +3L 17794390 17794435 -4.10890705137129 +3L 17794435 17794449 0 +3L 17794449 17794573 0 +3L 17794573 17795234 -2.84332836994737 +3L 17795234 17795277 -1.39256885351249 +3L 17795277 17795320 0 +3L 17795320 17795425 -2.68216688920208 +3L 17795425 17795983 -3.62972994085222 +3L 17795983 17796058 -1.39256885351249 +3L 17796058 17796073 0 +3L 17796073 17796087 -2.0877734257891 +3L 17796087 17796592 -2.4967599966257 +3L 17796592 17796710 -2.0877734257891 +3L 17796710 17796715 -3.42663311503462 +3L 17796715 17797122 -4.63459580306002 +3L 17797122 17797343 -2.24908659244355 +3L 17797343 17797478 0.195213662174823 +3L 17797478 17798180 -3.34207230947735 +3L 17798180 17798971 -3.42663311503462 +3L 17798971 17799079 -2.68216688920208 +3L 17799079 17799352 -2.90724258636772 +3L 17799352 17799359 -2.90724258636772 +3L 17799359 17799488 -1.66324290730034 +3L 17799488 17799542 -2.90724258636772 +3L 17799542 17799718 -3.91526259934797 +3L 17799718 17799725 -4.10890705137129 +3L 17799725 17800253 -2.21948744191938 +3L 17800253 17800362 0 +3L 17800362 17800387 0 +3L 17800387 17800446 0 +3L 17800446 17800776 0 +3L 17800776 17800817 0 +3L 17800817 17800941 0 +3L 17800941 17801230 0 +3L 17801230 17802045 1.84900954716158 +3L 17802045 17803003 -1.61607188218705 +3L 17803003 17803251 -2.20299160853309 +3L 17803251 17803296 -2.1697838140285 +3L 17803296 17803425 0 +3L 17803425 17803567 -4.83959573744425 +3L 17803567 17803787 -5.81418010758948 +3L 17803787 17804113 0.0813916759995955 +3L 17804113 17804178 0.677779557517565 +3L 17804178 17804719 -0.491970436125132 +3L 17804719 17805069 -2.0877734257891 +3L 17805069 17805436 -3.10798368064812 +3L 17805436 17805712 0.34173270480825 +3L 17805712 17805868 0.653546233646707 +3L 17805868 17806231 -0.819416713040243 +3L 17806231 17806931 1.1359676808816 +3L 17806931 17807063 0 +3L 17807063 17807255 0 +3L 17807255 17807416 -2.68216688920208 +3L 17807416 17807442 -4.10890705137129 +3L 17807442 17807688 -4.10890705137129 +3L 17807688 17807711 -3.1901861920131 +3L 17807711 17807766 0 +3L 17807766 17807864 0 +3L 17807864 17807915 0 +3L 17807915 17807994 0 +3L 17807994 17808044 0 +3L 17808044 17808115 0 +3L 17808115 17808258 -2.90724258636772 +3L 17808258 17808493 -4.75527859770636 +3L 17808493 17808814 -4.91925594869484 +3L 17808814 17808873 -3.62972994085222 +3L 17808873 17808914 0 +3L 17808914 17808940 0 +3L 17808940 17809104 -1.05910489381442 +3L 17809104 17809310 -2.0877734257891 +3L 17809310 17809337 -2.0877734257891 +3L 17809337 17809419 -1.39256885351249 +3L 17809419 17809431 -1.39256885351249 +3L 17809431 17809438 -2.0877734257891 +3L 17809438 17809668 -2.0877734257891 +3L 17809668 17809752 -1.39256885351249 +3L 17809752 17809803 -4.91925594869484 +3L 17809803 17810058 -5.8952170740068 +3L 17810058 17810145 -4.91925594869484 +3L 17810145 17810157 0 +3L 17810157 17810200 0 +3L 17810200 17810224 0 +3L 17810224 17810268 -2.0877734257891 +3L 17810268 17810862 -1.42576257804752 +3L 17810862 17810879 -0.477701727681987 +3L 17810879 17810897 -0.477701727681987 +3L 17810897 17810948 -0.465831973201576 +3L 17810948 17811356 -2.26084893649928 +3L 17811356 17812039 -4.10890705137129 +3L 17812039 17812148 -3.56515643745284 +3L 17812148 17812358 -4.75527859770636 +3L 17812358 17812499 -3.807732051728 +3L 17812499 17812882 -5.0896239740736 +3L 17812882 17813056 -4.69619873626645 +3L 17813056 17813068 -4.91925594869484 +3L 17813068 17813079 -4.91925594869484 +3L 17813079 17813376 -4.88791768737353 +3L 17813376 17813508 -3.42663311503462 +3L 17813508 17815750 0 +3L 17815750 17815876 -2.0877734257891 +3L 17815876 17816477 -3.56515643745284 +3L 17816477 17816517 -2.0877734257891 +3L 17816517 17816564 -1.39256885351249 +3L 17816564 17816645 -1.39256885351249 +3L 17816645 17816851 -1.78164771235138 +3L 17816851 17817000 0 +3L 17817000 17817229 -2.90724258636772 +3L 17817229 17817278 -2.90724258636772 +3L 17817278 17817298 -2.90724258636772 +3L 17817298 17817387 -1.66324290730034 +3L 17817387 17817507 0 +3L 17817507 17817604 0 +3L 17817604 17818316 0 +3L 17818316 17818442 0 +3L 17818442 17818572 0 +3L 17818572 17818827 -0.722626030916987 +3L 17818827 17819040 -2.0877734257891 +3L 17819040 17819113 -1.39256885351249 +3L 17819113 17819167 0 +3L 17819167 17819209 0 +3L 17819209 17819602 0 +3L 17819602 17819683 -1.39256885351249 +3L 17819683 17820560 -2.58530095496211 +3L 17820560 17820566 -2.90724258636772 +3L 17820566 17820638 -2.55492586585599 +3L 17820638 17820817 -2.0877734257891 +3L 17820817 17820901 -2.68216688920208 +3L 17820901 17820914 -3.807732051728 +3L 17820914 17820952 -3.807732051728 +3L 17820952 17821270 -3.807732051728 +3L 17821270 17821384 -2.4153674931314 +3L 17821384 17821485 0 +3L 17821485 17821512 0 +3L 17821512 17823646 0 +3L 17823646 17825165 0 +3L 17825165 17826011 0 +3L 17826011 17828139 -0.642010754900314 +3L 17828139 17828575 -0.550325869629912 +3L 17828575 17829505 -1.17015064642504 +3L 17829505 17829782 0 +3L 17829782 17831410 -1.43088984003895 +3L 17831410 17831745 -4.78959902445488 +3L 17831745 17832334 -5.11239941900137 +3L 17832334 17833362 -2.35553101928376 +3L 17833362 17833548 0 +3L 17833548 17833615 -2.0877734257891 +3L 17833615 17834528 -2.58530095496211 +3L 17834528 17834909 -1.89106755736324 +3L 17834909 17835095 0 +3L 17835095 17835162 0 +3L 17835162 17836015 0 +3L 17836015 17836031 0 +3L 17836031 17836037 0 +3L 17836037 17836749 -2.18474136668502 +3L 17836749 17836929 -3.42663311503462 +3L 17836929 17837458 -3.05563971869915 +3L 17837458 17837881 -3.42663311503462 +3L 17837881 17837898 -2.0877734257891 +3L 17837898 17837947 -1.39256885351249 +3L 17837947 17838062 0 +3L 17838062 17838165 0 +3L 17838165 17838555 -2.0877734257891 +3L 17838555 17838559 -2.90724258636772 +3L 17838559 17838573 -2.90724258636772 +3L 17838573 17838727 -2.68216688920208 +3L 17838727 17839378 -3.62972994085222 +3L 17839378 17839838 -4.87427577601333 +3L 17839838 17840659 1.11624373549678 +3L 17840659 17840913 -1.73331756706611 +3L 17840913 17840924 -4.91925594869484 +3L 17840924 17840961 0 +3L 17840961 17841131 0 +3L 17841131 17841735 2.75945610189416 +3L 17841735 17841749 2.97408106030375 +3L 17841749 17841800 0 +3L 17841800 17841911 0 +3L 17841911 17841918 0 +3L 17841918 17841938 0 +3L 17841938 17841942 0 +3L 17841942 17841993 2.71184381903969 +3L 17841993 17842004 2.97408106030375 +3L 17842004 17842240 2.97408106030375 +3L 17842240 17842320 2.01663924676308 +3L 17842320 17842479 -2.0877734257891 +3L 17842479 17842545 -2.0877734257891 +3L 17842545 17843716 -0.992229304112676 +3L 17843716 17844247 -0.421783257901303 +3L 17844247 17844352 -4.46799333170228 +3L 17844352 17844517 -5.60525984962007 +3L 17844517 17844746 -6.98837746081263 +3L 17844746 17845324 -1.43970343689831 +3L 17845324 17845556 -1.90502139127575 +3L 17845556 17845836 -3.111186905227 +3L 17845836 17845985 -4.69619873626645 +3L 17845985 17846055 0 +3L 17846055 17846128 0 +3L 17846128 17846941 -1.66324290730034 +3L 17846941 17847245 -3.5113324334873 +3L 17847245 17847355 -2.55492586585599 +3L 17847355 17847376 -1.39256885351249 +3L 17847376 17847806 -1.89106755736324 +3L 17847806 17847848 -1.39256885351249 +3L 17847848 17848138 -2.0877734257891 +3L 17848138 17848297 -1.05910489381442 +3L 17848297 17848462 0 +3L 17848462 17848561 0 +3L 17848561 17848900 -2.68216688920208 +3L 17848900 17850390 -0.38913829674294 +3L 17850390 17850452 -2.13916208570192 +3L 17850452 17850772 -2.24506652947103 +3L 17850772 17850932 -2.50749373447316 +3L 17850932 17851216 -4.68408678906606 +3L 17851216 17851261 -4.75527859770636 +3L 17851261 17851493 -3.31324822703235 +3L 17851493 17851570 -2.0877734257891 +3L 17851570 17852832 -1.89106755736324 +3L 17852832 17852874 0 +3L 17852874 17852892 0 +3L 17852892 17853050 0 +3L 17853050 17853195 0 +3L 17853195 17853269 0 +3L 17853269 17853323 0 +3L 17853323 17853459 -1.05910489381442 +3L 17853459 17853549 -2.0877734257891 +3L 17853549 17853741 1.89754591271408 +3L 17853741 17854148 1.59970175116322 +3L 17854148 17854229 -2.0877734257891 +3L 17854229 17854424 0 +3L 17854424 17854580 -3.56515643745284 +3L 17854580 17854602 -5.06648320422561 +3L 17854602 17854903 -5.06648320422561 +3L 17854903 17854988 -6.62245383238921 +3L 17854988 17855008 -7.35355047739808 +3L 17855008 17855328 -7.11501315791427 +3L 17855328 17855456 -2.0877734257891 +3L 17855456 17855764 -1.36031120560441 +3L 17855764 17856038 -1.24232725740425 +3L 17856038 17857266 -1.36821335855061 +3L 17857266 17857618 -3.1901861920131 +3L 17857618 17857696 0.138248241579792 +3L 17857696 17858185 -0.0847757186743665 +3L 17858185 17858323 -2.90724258636772 +3L 17858323 17858931 1.55052205260555 +3L 17858931 17858937 1.68215682577453 +3L 17858937 17858952 1.68215682577453 +3L 17858952 17859435 0.408296077472383 +3L 17859435 17859615 1.63522137105823 +3L 17859615 17859724 1.54948680669619 +3L 17859724 17859747 0 +3L 17859747 17859902 0 +3L 17859902 17860636 -4.05349300925994 +3L 17860636 17861050 -0.90207925846968 +3L 17861050 17861152 2.01663924676308 +3L 17861152 17861184 0 +3L 17861184 17861281 0 +3L 17861281 17862030 -2.78884325466564 +3L 17862030 17862796 -0.670002262006125 +3L 17862796 17863190 0 +3L 17863190 17863220 -2.0877734257891 +3L 17863220 17864153 -3.04475103245643 +3L 17864153 17864516 0 +3L 17864516 17864836 0 +3L 17864836 17865151 0 +3L 17865151 17866093 0 +3L 17866093 17866175 0 +3L 17866175 17866305 0 +3L 17866305 17866378 0 +3L 17866378 17866694 0 +3L 17866694 17867011 0 +3L 17867011 17868570 0 +3L 17868570 17869183 -2.68216688920208 +3L 17869183 17869198 -3.42663311503462 +3L 17869198 17869641 1.10596612513559 +3L 17869641 17869821 -3.88912196962424 +3L 17869821 17869930 -3.42663311503462 +3L 17869930 17869953 -1.39256885351249 +3L 17869953 17870107 2.5016839262 +3L 17870107 17870406 3.02456717647289 +3L 17870406 17870956 2.31862062585054 +3L 17870956 17871537 0.0606070532697907 +3L 17871537 17872304 -0.37355103300116 +3L 17872304 17872698 2.67814000004346 +3L 17872698 17872728 5.06185448609285 +3L 17872728 17873245 3.04619503201503 +3L 17873245 17873275 1.63522137105823 +3L 17873275 17873789 -4.10890705137129 +3L 17873789 17874126 -2.29310500926643 +3L 17874126 17874155 -3.42663311503462 +3L 17874155 17874464 -4.05349300925994 +3L 17874464 17874546 -2.90724258636772 +3L 17874546 17875152 2.75945610189416 +3L 17875152 17875158 2.97408106030375 +3L 17875158 17875173 2.97408106030375 +3L 17875173 17875836 3.94730805922951 +3L 17875836 17875945 0 +3L 17875945 17875968 0 +3L 17875968 17876123 0 +3L 17876123 17876387 -1.39256885351249 +3L 17876387 17877117 -2.44176656376487 +3L 17877117 17877226 -1.66324290730034 +3L 17877226 17877249 -2.0877734257891 +3L 17877249 17878063 -3.23238403846164 +3L 17878063 17878070 -3.807732051728 +3L 17878070 17878233 -5.39842757918527 +3L 17878233 17878596 -7.40057160428963 +3L 17878596 17879130 -5.46187636685508 +3L 17879130 17879164 -4.91925594869484 +3L 17879164 17879329 -3.96616443097239 +3L 17879329 17879711 -3.24310848589432 +3L 17879711 17879763 -0.982939616111517 +3L 17879763 17880004 0.260230321403625 +3L 17880004 17880293 0.693457228644697 +3L 17880293 17880323 0.593861154390571 +3L 17880323 17880431 0 +3L 17880431 17880556 0 +3L 17880556 17880609 0 +3L 17880609 17880747 0 +3L 17880747 17881100 -3.10189685114339 +3L 17881100 17881160 -3.1901861920131 +3L 17881160 17881167 -2.0877734257891 +3L 17881167 17881238 -1.39256885351249 +3L 17881238 17881344 -1.39256885351249 +3L 17881344 17881833 -3.48014999994755 +3L 17881833 17881894 -2.0877734257891 +3L 17881894 17881949 0 +3L 17881949 17882159 0 +3L 17882159 17882261 0 +3L 17882261 17882272 0 +3L 17882272 17882285 0 +3L 17882285 17882351 0 +3L 17882351 17882408 0 +3L 17882408 17882471 0 +3L 17882471 17882483 0 +3L 17882483 17882647 -2.0877734257891 +3L 17882647 17882850 -3.42663311503462 +3L 17882850 17883026 -3.42663311503462 +3L 17883026 17883721 -1.20186343010074 +3L 17883721 17884138 -3.07564608326844 +3L 17884138 17884887 -3.30248457954299 +3L 17884887 17885181 -5.06648320422561 +3L 17885181 17885409 -5.10105663981123 +3L 17885409 17885431 -4.91925594869484 +3L 17885431 17885679 -3.72147453126565 +3L 17885679 17885685 -3.807732051728 +3L 17885685 17886130 -3.75907525084849 +3L 17886130 17886319 -2.55492586585599 +3L 17886319 17886489 0.878621930812341 +3L 17886489 17886781 1.25412243436485 +3L 17886781 17886825 1.25412243436485 +3L 17886825 17886891 -2.90724258636772 +3L 17886891 17887459 -3.94939591375755 +3L 17887459 17887719 -2.34017592319915 +3L 17887719 17887729 0 +3L 17887729 17887877 0 +3L 17887877 17887906 0 +3L 17887906 17888370 -2.4967599966257 +3L 17888370 17888428 -2.0877734257891 +3L 17888428 17888433 0 +3L 17888433 17888499 0 +3L 17888499 17889411 -1.80764503055372 +3L 17889411 17889563 -1.05910489381442 +3L 17889563 17889884 -2.0877734257891 +3L 17889884 17890010 -3.2733826212464 +3L 17890010 17890208 -4.61877510527196 +3L 17890208 17890565 -2.90724258636772 +3L 17890565 17890687 0 +3L 17890687 17890849 0 +3L 17890849 17890885 0 +3L 17890885 17890966 0 +3L 17890966 17891080 -1.39256885351249 +3L 17891080 17891092 -2.0877734257891 +3L 17891092 17891526 -1.73208775926791 +3L 17891526 17891622 0 +3L 17891622 17891699 0 +3L 17891699 17892087 -2.0877734257891 +3L 17892087 17892097 -2.90724258636772 +3L 17892097 17892242 -2.4153674931314 +3L 17892242 17892281 -2.0877734257891 +3L 17892281 17892405 -2.90724258636772 +3L 17892405 17892610 -4.17530933603838 +3L 17892610 17892910 -2.90724258636772 +3L 17892910 17893175 -2.34017592319915 +3L 17893175 17893622 -4.06954396796188 +3L 17893622 17894152 -4.07452348956241 +3L 17894152 17894641 -0.389830889240997 +3L 17894641 17895390 -1.85761054842779 +3L 17895390 17895404 -5.9719432181294 +3L 17895404 17895414 -5.9719432181294 +3L 17895414 17895779 -6.9975824378003 +3L 17895779 17895877 -6.47207677441559 +3L 17895877 17896160 -6.56417012142421 +3L 17896160 17896398 -6.07988491452019 +3L 17896398 17896448 -4.99474683265014 +3L 17896448 17897310 -4.39556391357549 +3L 17897310 17897631 -4.75527859770636 +3L 17897631 17897924 -4.23878936063925 +3L 17897924 17898229 -0.523524929314923 +3L 17898229 17899003 1.18642892770275 +3L 17899003 17899324 -2.29310500926643 +3L 17899324 17899380 -1.39256885351249 +3L 17899380 17899437 -2.90724258636772 +3L 17899437 17899540 -3.96616443097239 +3L 17899540 17899703 -3.10189685114339 +3L 17899703 17899764 -2.90724258636772 +3L 17899764 17900186 -3.65654337623127 +3L 17900186 17900257 -2.0877734257891 +3L 17900257 17900271 0 +3L 17900271 17900467 0 +3L 17900467 17900576 0 +3L 17900576 17900585 0 +3L 17900585 17900639 0 +3L 17900639 17900651 0 +3L 17900651 17900720 -1.27536207441091 +3L 17900720 17900744 -1.30049048607261 +3L 17900744 17901015 -1.33429515485634 +3L 17901015 17901405 -2.86141621373138 +3L 17901405 17901598 -4.75527859770636 +3L 17901598 17902051 -2.62114699228392 +3L 17902051 17902281 0.104359381974235 +3L 17902281 17902463 -1.13540991625203 +3L 17902463 17903052 -1.26924057419307 +3L 17903052 17903296 -2.27593798522636 +3L 17903296 17903514 -3.62972994085222 +3L 17903514 17903522 -2.90724258636772 +3L 17903522 17903530 -2.90724258636772 +3L 17903530 17903547 -2.90724258636772 +3L 17903547 17904735 -2.71054276640109 +3L 17904735 17904850 -1.66324290730034 +3L 17904850 17904882 0 +3L 17904882 17904967 0 +3L 17904967 17906145 0.865167444331671 +3L 17906145 17908245 -0.291840495777527 +3L 17908245 17908896 -2.47282073788985 +3L 17908896 17908918 0 +3L 17908918 17908977 0 +3L 17908977 17909126 0 +3L 17909126 17909181 -2.90724258636772 +3L 17909181 17909391 -4.35793833134167 +3L 17909391 17909754 -4.94486100099495 +3L 17909754 17909977 -4.61877510527196 +3L 17909977 17909998 -4.75527859770636 +3L 17909998 17910191 -3.96616443097239 +3L 17910191 17910259 -3.1901861920131 +3L 17910259 17910607 -1.23712259930451 +3L 17910607 17910718 1.87166829407975 +3L 17910718 17910788 4.10441267255218 +3L 17910788 17910833 0 +3L 17910833 17910897 0 +3L 17910897 17910975 0 +3L 17910975 17911177 -1.05910489381442 +3L 17911177 17911743 -3.42663311503462 +3L 17911743 17912087 -4.75527859770636 +3L 17912087 17912121 -3.807732051728 +3L 17912121 17912195 -6.11414243765132 +3L 17912195 17912582 -3.17833938888357 +3L 17912582 17915892 -1.63292514427073 +3L 17915892 17916295 -3.00785017442921 +3L 17916295 17916340 -2.0877734257891 +3L 17916340 17916430 0 +3L 17916430 17918198 0.723969350973548 +3L 17918198 17918586 1.61724655256022 +3L 17918586 17918807 0 +3L 17918807 17918958 0 +3L 17918958 17919338 0 +3L 17919338 17919537 0 +3L 17919537 17920271 -1.39256885351249 +3L 17920271 17920549 0 +3L 17920549 17920815 0 +3L 17920815 17921036 -2.0877734257891 +3L 17921036 17921871 -4.08607596748718 +3L 17921871 17921977 -5.72831892824831 +3L 17921977 17922243 -6.42918031327152 +3L 17922243 17922258 -6.47207677441559 +3L 17922258 17922400 -1.76153525843927 +3L 17922400 17922596 0.427258683032826 +3L 17922596 17922894 0.430287447829839 +3L 17922894 17923286 -2.4967599966257 +3L 17923286 17923365 -2.0877734257891 +3L 17923365 17923992 -1.48854786879403 +3L 17923992 17924164 0 +3L 17924164 17924169 0 +3L 17924169 17924174 0 +3L 17924174 17924184 0 +3L 17924184 17924273 -1.39256885351249 +3L 17924273 17924315 -2.0877734257891 +3L 17924315 17924463 -3.69153660144476 +3L 17924463 17924757 -4.78959902445488 +3L 17924757 17924793 -4.10890705137129 +3L 17924793 17925098 -2.90724258636772 +3L 17925098 17925183 -2.0877734257891 +3L 17925183 17925265 -2.55492586585599 +3L 17925265 17925275 -3.42663311503462 +3L 17925275 17925675 -3.96616443097239 +3L 17925675 17926039 -4.10890705137129 +3L 17926039 17926471 -2.90724258636772 +3L 17926471 17926970 -1.39256885351249 +3L 17926970 17927035 -1.39256885351249 +3L 17927035 17931125 -0.437915121923503 +3L 17931125 17931895 -3.06840310089978 +3L 17931895 17933298 -1.30036494577117 +3L 17933298 17934457 -3.34753161773866 +3L 17934457 17934869 -0.227389498648522 +3L 17934869 17934991 0.690595949639123 +3L 17934991 17935008 0 +3L 17935008 17935133 -3.56515643745284 +3L 17935133 17935431 -5.34557959308659 +3L 17935431 17935657 -5.56454704262344 +3L 17935657 17936054 -3.807732051728 +3L 17936054 17936188 -1.66324290730034 +3L 17936188 17936209 -2.90724258636772 +3L 17936209 17936813 -4.01533135855775 +3L 17936813 17937334 1.85454197726282 +3L 17937334 17938048 -2.0877734257891 +3L 17938048 17938250 -3.1901861920131 +3L 17938250 17938780 -0.998196487798505 +3L 17938780 17939346 -0.984092657790277 +3L 17939346 17940503 -0.267491566674919 +3L 17940503 17940539 0 +3L 17940539 17940592 0 +3L 17940592 17942272 0 +3L 17942272 17942312 0 +3L 17942312 17942321 0 +3L 17942321 17943114 -0.13074778143034 +3L 17943114 17943122 -1.41022228832274 +3L 17943122 17943628 -4.15540558373432 +3L 17943628 17943699 -4.75527859770636 +3L 17943699 17943792 -2.99927831812144 +3L 17943792 17944017 -3.03229782580838 +3L 17944017 17946841 -2.14281966648736 +3L 17946841 17947981 0.492458533894907 +3L 17947981 17948012 -2.90724258636772 +3L 17948012 17948162 -1.66324290730034 +3L 17948162 17948255 -1.05910489381442 +3L 17948255 17949476 -1.47452610711612 +3L 17949476 17950037 -0.616540069924127 +3L 17950037 17950764 0.274668934561904 +3L 17950764 17950915 -2.2809124867484 +3L 17950915 17951582 -1.03171768656281 +3L 17951582 17952835 0.783581298676581 +3L 17952835 17952848 -2.90724258636772 +3L 17952848 17952972 -2.90724258636772 +3L 17952972 17953555 -1.2101962777544 +3L 17953555 17953940 -1.76721350523057 +3L 17953940 17954149 -5.29271740587966 +3L 17954149 17955217 -3.17295528665581 +3L 17955217 17955708 -1.94281499061958 +3L 17955708 17955750 0 +3L 17955750 17955814 0 +3L 17955814 17955843 0 +3L 17955843 17955970 0 +3L 17955970 17956039 0 +3L 17956039 17956141 -1.05910489381442 +3L 17956141 17957541 -2.30314310659826 +3L 17957541 17957732 -5.32232407845369 +3L 17957732 17958051 -5.86334640627401 +3L 17958051 17958141 -0.962070531673424 +3L 17958141 17958595 -0.346964418007989 +3L 17958595 17959034 -0.151330194529063 +3L 17959034 17960089 0.160997041792354 +3L 17960089 17960121 -3.62972994085222 +3L 17960121 17960991 -2.18850082387457 +3L 17960991 17961638 -1.1773410317231 +3L 17961638 17961647 -6.52396367771668 +3L 17961647 17962713 -3.02715636847465 +3L 17962713 17963847 -3.45363969489135 +3L 17963847 17964617 -3.23238403846164 +3L 17964617 17965995 -2.25223925265358 +3L 17965995 17966794 -4.13138244802135 +3L 17966794 17966884 -2.55492586585599 +3L 17966884 17969476 -1.28974662020957 +3L 17969476 17970337 -2.61994870380323 +3L 17970337 17971015 -3.28946129483649 +3L 17971015 17971107 -2.0877734257891 +3L 17971107 17971611 -2.8268973771849 +3L 17971611 17971781 -1.66324290730034 +3L 17971781 17971788 0 +3L 17971788 17973305 -3.74551537704315 +3L 17973305 17973538 -8.01023317745542 +3L 17973538 17973657 -7.03382603726108 +3L 17973657 17973661 -2.0877734257891 +3L 17973661 17973665 -2.0877734257891 +3L 17973665 17973813 -1.05910489381442 +3L 17973813 17973821 0 +3L 17973821 17973876 0 +3L 17973876 17973978 0 +3L 17973978 17974106 0 +3L 17974106 17974265 0 +3L 17974265 17975341 -2.15513305219635 +3L 17975341 17976227 -1.68798486539784 +3L 17976227 17977264 -0.324237692117883 +3L 17977264 17977269 -2.0877734257891 +3L 17977269 17978366 -1.78164771235138 +3L 17978366 17978477 -2.0877734257891 +3L 17978477 17978559 -2.0877734257891 +3L 17978559 17978566 -2.0877734257891 +3L 17978566 17979060 -2.23735919244861 +3L 17979060 17979249 -1.78164771235138 +3L 17979249 17980457 -1.95175728612821 +3L 17980457 17982440 -1.60732909864263 +3L 17982440 17982873 -0.825852039016662 +3L 17982873 17983195 -1.2591629460967 +3L 17983195 17983389 -3.53175247535673 +3L 17983389 17983578 -4.10890705137129 +3L 17983578 17984571 -2.86188032349347 +3L 17984571 17984876 -3.59133367705869 +3L 17984876 17984966 -2.55492586585599 +3L 17984966 17985966 -1.92778904339385 +3L 17985966 17986452 0.127982152511143 +3L 17986452 17988125 -0.37332161401315 +3L 17988125 17988321 -2.90724258636772 +3L 17988321 17988412 -2.0877734257891 +3L 17988412 17990803 -0.37355103300116 +3L 17990803 17991080 -3.5113324334873 +3L 17991080 17991166 -3.807732051728 +3L 17991166 17991678 -3.48014999994755 +3L 17991678 17991741 -2.0877734257891 +3L 17991741 17992029 -2.0877734257891 +3L 17992029 17992042 -2.0877734257891 +3L 17992042 17992122 1.19717008618447 +3L 17992122 17992423 0.34173270480825 +3L 17992423 17992433 -0.611227465071405 +3L 17992433 17992467 -0.882148566416277 +3L 17992467 17992737 -0.746639534218448 +3L 17992737 17993492 -0.397880698557361 +3L 17993492 17993796 -6.12762087661927 +3L 17993796 17993828 -5.9719432181294 +3L 17993828 17993834 -5.9719432181294 +3L 17993834 17993960 -4.43223310796149 +3L 17993960 17994035 0 +3L 17994035 17994085 0 +3L 17994085 17994211 -2.90724258636772 +3L 17994211 17995181 -0.722032901538044 +3L 17995181 17995510 -3.62972994085222 +3L 17995510 17996059 -3.62972994085222 +3L 17996059 17996144 -3.42663311503462 +3L 17996144 17996448 -1.20186343010074 +3L 17996448 17997007 -1.48854786879403 +3L 17997007 17997015 -4.10890705137129 +3L 17997015 17997022 -4.10890705137129 +3L 17997022 17997363 -4.1622711318164 +3L 17997363 17997406 -0.253525658789485 +3L 17997406 17997670 0.443079380820889 +3L 17997670 17997675 0.306575888386487 +3L 17997675 17997933 0.578218790261109 +3L 17997933 17998013 0.306575888386487 +3L 17998013 17998098 -2.0877734257891 +3L 17998098 17998113 0 +3L 17998113 17998264 -2.90724258636772 +3L 17998264 17998773 -4.07452348956241 +3L 17998773 17998862 -2.0877734257891 +3L 17998862 17998993 0 +3L 17998993 17999085 -4.10890705137129 +3L 17999085 17999429 0.203562944359285 +3L 17999429 17999479 2.15461189972513 +3L 17999479 17999536 0 +3L 17999536 17999644 -3.69153660144476 +3L 17999644 18000097 -5.993135824113 +3L 18000097 18000110 -5.72831892824831 +3L 18000110 18000240 -1.13540991625203 +3L 18000240 18000660 0.401664203993259 +3L 18000660 18000672 0.306575888386487 +3L 18000672 18000793 0.287406198768024 +3L 18000793 18000849 -1.57898014246747 +3L 18000849 18000878 -1.60743264278744 +3L 18000878 18001117 -1.66391406964304 +3L 18001117 18001122 -1.82119125748459 +3L 18001122 18001544 -1.80451591420941 +3L 18001544 18001556 0.491609840339092 +3L 18001556 18001569 0.491609840339092 +3L 18001569 18001777 -0.0392021537183818 +3L 18001777 18001906 -0.454542538543331 +3L 18001906 18002017 -3.96616443097239 +3L 18002017 18002650 -0.766942685560806 +3L 18002650 18004439 0.283862157594417 +3L 18004439 18004843 -1.31795045205158 +3L 18004843 18004848 -5.63702259370271 +3L 18004848 18004866 -5.63702259370271 +3L 18004866 18004928 -2.10787627925341 +3L 18004928 18005238 -1.57951549099154 +3L 18005238 18005522 -2.99210888828711 +3L 18005522 18005690 -4.69619873626645 +3L 18005690 18005717 0.138248241579792 +3L 18005717 18006280 0.87771536776208 +3L 18006280 18006454 -3.31324822703235 +3L 18006454 18006671 -2.4153674931314 +3L 18006671 18007140 -3.05563971869915 +3L 18007140 18007266 -2.55492586585599 +3L 18007266 18007284 -3.1901861920131 +3L 18007284 18008049 -3.4657502557624 +3L 18008049 18008194 -3.10189685114339 +3L 18008194 18008224 -3.42663311503462 +3L 18008224 18008257 -4.35793833134167 +3L 18008257 18008499 -4.91925594869484 +3L 18008499 18008596 -3.96616443097239 +3L 18008596 18008666 -1.39256885351249 +3L 18008666 18008843 -2.0877734257891 +3L 18008843 18008966 -1.05910489381442 +3L 18008966 18009043 2.71184381903969 +3L 18009043 18009131 2.97408106030375 +3L 18009131 18009888 0.557738985659875 +3L 18009888 18010077 -5.43502497407227 +3L 18010077 18010213 -5.20006743060371 +3L 18010213 18011410 -1.80157002655858 +3L 18011410 18011614 -2.68216688920208 +3L 18011614 18012196 -0.0219375265094411 +3L 18012196 18012451 -0.563840168048818 +3L 18012451 18012488 -1.24232725740425 +3L 18012488 18012494 -1.24232725740425 +3L 18012494 18013141 -2.22732684451591 +3L 18013141 18013467 -3.73914355418382 +3L 18013467 18013508 -2.55492586585599 +3L 18013508 18013608 -4.23878936063925 +3L 18013608 18013626 -5.20006743060371 +3L 18013626 18013846 -5.23161650549329 +3L 18013846 18014191 -3.00785017442921 +3L 18014191 18014264 -0.253525658789485 +3L 18014264 18014977 0.298348599527623 +3L 18014977 18015153 -2.01412904897499 +3L 18015153 18015359 -3.2536203089283 +3L 18015359 18015662 -2.32619471520932 +3L 18015662 18016765 -0.258112953962787 +3L 18016765 18016884 1.47301539422532 +3L 18016884 18016951 1.25412243436485 +3L 18016951 18017196 -0.309602338376222 +3L 18017196 18017297 1.43212454524063 +3L 18017297 18017386 2.01663924676308 +3L 18017386 18017492 -1.05910489381442 +3L 18017492 18017737 0 +3L 18017737 18018167 -1.3619660672629 +3L 18018167 18018993 1.650056502935 +3L 18018993 18019412 -2.31660683260668 +3L 18019412 18019434 -4.10890705137129 +3L 18019434 18019447 0 +3L 18019447 18019500 0 +3L 18019500 18019518 -2.0877734257891 +3L 18019518 18020949 -0.982053164881759 +3L 18020949 18021216 -3.1901861920131 +3L 18021216 18021389 -4.57024464575376 +3L 18021389 18021552 -4.43223310796149 +3L 18021552 18021652 -1.15807772649368 +3L 18021652 18021656 -0.666464442155458 +3L 18021656 18021872 -1.13453340848078 +3L 18021872 18022069 -2.99927831812144 +3L 18022069 18022553 -0.00467060013214633 +3L 18022553 18025392 0.369866769128838 +3L 18025392 18025523 0 +3L 18025523 18027877 0.0721974830852724 +3L 18027877 18028547 0.296680620824183 +3L 18028547 18028667 -2.4153674931314 +3L 18028667 18030169 -0.38913829674294 +3L 18030169 18030945 -2.72578845835579 +3L 18030945 18033431 -1.58876529241921 +3L 18033431 18033453 -2.90724258636772 +3L 18033453 18034212 -0.252181177778517 +3L 18034212 18036026 0.723969350973548 +3L 18036026 18036121 4.10441267255218 +3L 18036121 18037283 0.963921671753319 +3L 18037283 18037346 -3.807732051728 +3L 18037346 18037717 -0.906065716995537 +3L 18037717 18037887 -1.38430626412503 +3L 18037887 18037987 -2.90724258636772 +3L 18037987 18038134 -2.68216688920208 +3L 18038134 18038305 -3.56515643745284 +3L 18038305 18039112 -3.59780443920947 +3L 18039112 18040793 1.26829931482105 +3L 18040793 18041245 0.937058500033116 +3L 18041245 18041411 -2.74181239868136 +3L 18041411 18042269 -3.1901861920131 +3L 18042269 18042412 -2.68216688920208 +3L 18042412 18042603 -3.88912196962424 +3L 18042603 18042656 -4.75527859770636 +3L 18042656 18043216 -0.386591832371838 +3L 18043216 18043399 -1.63927861058929 +3L 18043399 18043781 -6.62245383238921 +3L 18043781 18044444 -0.137775903068753 +3L 18044444 18045188 -1.37273728153607 +3L 18045188 18045421 -2.74181239868136 +3L 18045421 18045459 -2.55492586585599 +3L 18045459 18045661 -3.56515643745284 +3L 18045661 18045668 -3.807732051728 +3L 18045668 18046077 -4.77987621570102 +3L 18046077 18046296 -6.39057025028289 +3L 18046296 18046398 -5.37977474696309 +3L 18046398 18047941 -0.798128779275947 +3L 18047941 18049103 0 +3L 18049103 18049116 0 +3L 18049116 18049683 0 +3L 18049683 18049728 0 +3L 18049728 18049887 0 +3L 18049887 18050778 0 +3L 18050778 18051009 0 +3L 18051009 18051531 0 +3L 18051531 18052023 0 +3L 18052023 18052104 0 +3L 18052104 18052408 0 +3L 18052408 18052430 0 +3L 18052430 18052490 0 +3L 18052490 18052613 0 +3L 18052613 18052737 0 +3L 18052737 18053448 0 +3L 18053448 18053565 0 +3L 18053565 18054284 -3.42663311503462 +3L 18054284 18054471 -2.55492586585599 +3L 18054471 18054640 0 +3L 18054640 18055202 0 +3L 18055202 18055284 0 +3L 18055284 18055782 0 +3L 18055782 18056512 0 +3L 18056512 18058159 -0.190733851261496 +3L 18058159 18058241 -3.1901861920131 +3L 18058241 18058498 -3.62972994085222 +3L 18058498 18058539 -3.1901861920131 +3L 18058539 18058635 -2.0877734257891 +3L 18058635 18058817 -0.858163070035469 +3L 18058817 18058851 0 +3L 18058851 18058873 0 +3L 18058873 18058965 0 +3L 18058965 18059182 -1.78164771235138 +3L 18059182 18059490 -3.59133367705869 +3L 18059490 18059506 -4.10890705137129 +3L 18059506 18059539 -4.10890705137129 +3L 18059539 18059718 -5.35133536722057 +3L 18059718 18060013 -6.96791737162587 +3L 18060013 18060128 -5.9719432181294 +3L 18060128 18060228 -4.27960886543453 +3L 18060228 18060522 -5.63702259370271 +3L 18060522 18063516 0.0430268835613151 +3L 18063516 18063857 -1.77407220597603 +3L 18063857 18064223 -0.880525970809334 +3L 18064223 18064588 -0.434548644168779 +3L 18064588 18065474 -0.938552525746051 +3L 18065474 18065704 -0.618655925434187 +3L 18065704 18065971 -0.693640388218044 +3L 18065971 18066779 -3.4975567130837 +3L 18066779 18067344 -0.548027260999803 +3L 18067344 18069091 -0.483684726021707 +3L 18069091 18069293 -4.27960886543453 +3L 18069293 18069444 -1.3058520039036 +3L 18069444 18069696 -2.4695409440924 +3L 18069696 18069899 -3.76379532813359 +3L 18069899 18070068 -3.56515643745284 +3L 18070068 18070081 -4.35793833134167 +3L 18070081 18070778 -3.93583934205124 +3L 18070778 18071197 -3.42663311503462 +3L 18071197 18071809 -4.63459580306002 +3L 18071809 18071826 -4.91925594869484 +3L 18071826 18072828 -1.93261292803933 +3L 18072828 18073004 -2.4153674931314 +3L 18073004 18073782 -4.57024464575376 +3L 18073782 18074410 -3.3180034436069 +3L 18074410 18074731 -1.81281910895569 +3L 18074731 18075001 -2.20506859775565 +3L 18075001 18075456 -5.35543262059093 +3L 18075456 18075585 -4.46799333170228 +3L 18075585 18075590 -3.42663311503462 +3L 18075590 18076472 -2.32225301248307 +3L 18076472 18077914 -1.47923225488917 +3L 18077914 18078699 -1.21431268730084 +3L 18078699 18079036 0.205970934384869 +3L 18079036 18079202 0.303204663390384 +3L 18079202 18079650 0.376429369331541 +3L 18079650 18079739 2.15461189972513 +3L 18079739 18079904 0.458891968871037 +3L 18079904 18080913 -0.270494190097244 +3L 18080913 18081017 -2.90724258636772 +3L 18081017 18081939 -3.04475103245643 +3L 18081939 18084040 -1.59350166000434 +3L 18084040 18084591 -6.17220279159045 +3L 18084591 18085069 -3.62972994085222 +3L 18085069 18085505 -4.1472245999285 +3L 18085505 18085519 -2.0877734257891 +3L 18085519 18086805 -1.02441121121987 +3L 18086805 18086834 -0.00462871813275598 +3L 18086834 18087339 -1.55310481000001 +3L 18087339 18087489 -4.91925594869484 +3L 18087489 18087936 -1.07010019845582 +3L 18087936 18087970 -0.477701727681987 +3L 18087970 18087976 -0.465831973201576 +3L 18087976 18089258 -1.89106755736324 +3L 18089258 18089295 -2.90724258636772 +3L 18089295 18089327 -2.90724258636772 +3L 18089327 18089363 -2.90724258636772 +3L 18089363 18089562 -3.62972994085222 +3L 18089562 18089587 -4.75527859770636 +3L 18089587 18089803 -4.7111970626043 +3L 18089803 18090086 -4.5302095028321 +3L 18090086 18090370 -4.10890705137129 +3L 18090370 18090601 -2.04681974110371 +3L 18090601 18090906 -1.49289054653492 +3L 18090906 18091539 -4.71616205407571 +3L 18091539 18091688 -3.2733826212464 +3L 18091688 18092433 -2.8492571817686 +3L 18092433 18092487 -3.1901861920131 +3L 18092487 18092560 -2.0877734257891 +3L 18092560 18093068 -2.23735919244861 +3L 18093068 18093510 -3.75907525084849 +3L 18093510 18093781 -3.02715636847465 +3L 18093781 18093982 -3.96616443097239 +3L 18093982 18094751 -4.4190064418029 +3L 18094751 18095041 -6.41825390098518 +3L 18095041 18095200 -4.91925594869484 +3L 18095200 18095415 -2.55492586585599 +3L 18095415 18095476 -3.62972994085222 +3L 18095476 18096437 -4.02274895893402 +3L 18096437 18096806 -3.91526259934797 +3L 18096806 18096858 -2.90724258636772 +3L 18096858 18097566 -4.40297269374249 +3L 18097566 18097876 -5.15689923582141 +3L 18097876 18098677 -4.45832811908613 +3L 18098677 18098831 -4.10890705137129 +3L 18098831 18098895 -2.0877734257891 +3L 18098895 18099809 0.465079398630601 +3L 18099809 18100374 -3.84447332790527 +3L 18100374 18100661 -3.5113324334873 +3L 18100661 18103984 -1.00819699474706 +3L 18103984 18105391 -3.13786351330112 +3L 18105391 18105423 -4.75527859770636 +3L 18105423 18105847 -2.32864184726809 +3L 18105847 18106276 -1.85565337727462 +3L 18106276 18107038 -0.00356355365022258 +3L 18107038 18107134 -4.10890705137129 +3L 18107134 18108071 -4.26717158112086 +3L 18108071 18108136 -2.90724258636772 +3L 18108136 18108218 0 +3L 18108218 18108316 0 +3L 18108316 18108337 0 +3L 18108337 18108392 -3.96616443097239 +3L 18108392 18108822 -5.60525984962007 +3L 18108822 18108854 -3.96616443097239 +3L 18108854 18109794 -2.58530095496211 +3L 18109794 18109992 -5.72831892824831 +3L 18109992 18110536 -4.81897392259593 +3L 18110536 18110880 -5.32232407845369 +3L 18110880 18110932 0.296680620824183 +3L 18110932 18110991 2.97408106030375 +3L 18110991 18111360 1.70553589342993 +3L 18111360 18111519 -1.66324290730034 +3L 18111519 18111660 0 +3L 18111660 18111788 0 +3L 18111788 18111958 0 +3L 18111958 18111973 0 +3L 18111973 18112166 0 +3L 18112166 18112224 -3.807732051728 +3L 18112224 18112606 -5.53955621377484 +3L 18112606 18113243 -3.52044370844044 +3L 18113243 18113412 -2.55492586585599 +3L 18113412 18113991 -3.67633139673871 +3L 18113991 18114297 -4.91925594869484 +3L 18114297 18114631 -1.21517046020173 +3L 18114631 18114873 -1.77155672248989 +3L 18114873 18114903 -3.1901861920131 +3L 18114903 18115050 0 +3L 18115050 18115188 0 +3L 18115188 18115227 0 +3L 18115227 18115384 -2.4153674931314 +3L 18115384 18115732 -4.75527859770636 +3L 18115732 18115875 -3.42663311503462 +3L 18115875 18117669 -0.746118358928867 +3L 18117669 18118073 -2.8150453150792 +3L 18118073 18119550 -0.982053164881759 +3L 18119550 18119830 -1.21234200048783 +3L 18119830 18119839 -1.97197155116823 +3L 18119839 18120240 -4.15656641309002 +3L 18120240 18120298 -5.8952170740068 +3L 18120298 18120337 -4.91925594869484 +3L 18120337 18120490 0 +3L 18120490 18120586 -1.39256885351249 +3L 18120586 18120729 -3.10189685114339 +3L 18120729 18121028 -3.87320717616847 +3L 18121028 18121118 -2.90724258636772 +3L 18121118 18121169 -2.90724258636772 +3L 18121169 18121766 -4.2249221066907 +3L 18121766 18121830 -5.20006743060371 +3L 18121830 18122118 -5.83075574626064 +3L 18122118 18122415 -4.35793833134167 +3L 18122415 18122731 -4.78959902445488 +3L 18122731 18122777 -5.13482089478576 +3L 18122777 18122807 -5.43502497407227 +3L 18122807 18123134 -5.49864801902651 +3L 18123134 18125416 -1.33563555857024 +3L 18125416 18125682 -2.0877734257891 +3L 18125682 18126068 -5.06648320422561 +3L 18126068 18126203 -3.807732051728 +3L 18126203 18126971 -2.65985131681495 +3L 18126971 18127447 -0.621593005588116 +3L 18127447 18128151 -5.80582032059082 +3L 18128151 18128367 -4.23878936063925 +3L 18128367 18128566 -3.88912196962424 +3L 18128566 18129562 -2.90724258636772 +3L 18129562 18129849 -4.64712862306423 +3L 18129849 18130219 -4.78393535212557 +3L 18130219 18130410 -4.10890705137129 +3L 18130410 18130431 -4.35793833134167 +3L 18130431 18130803 -4.83959573744425 +3L 18130803 18131537 -4.66572604956879 +3L 18131537 18131560 -3.42663311503462 +3L 18131560 18131639 2.71184381903969 +3L 18131639 18131669 2.97408106030375 +3L 18131669 18131854 3.28020677374147 +3L 18131854 18131873 5.06185448609285 +3L 18131873 18131987 4.10441267255218 +3L 18131987 18132144 0 +3L 18132144 18132315 -3.42663311503462 +3L 18132315 18132613 -5.36846617965385 +3L 18132613 18133074 -4.65246668120586 +3L 18133074 18133558 -3.42663311503462 +3L 18133558 18134814 -1.17902990514451 +3L 18134814 18135116 0 +3L 18135116 18135245 0 +3L 18135245 18136204 -2.3067112819453 +3L 18136204 18136587 -4.66572604956879 +3L 18136587 18136961 -2.4153674931314 +3L 18136961 18137813 -1.78164771235138 +3L 18137813 18139313 -3.02166662541891 +3L 18139313 18139356 -3.62972994085222 +3L 18139356 18139713 -3.42663311503462 +3L 18139713 18139828 -2.0877734257891 +3L 18139828 18139862 0 +3L 18139862 18139920 0 +3L 18139920 18141102 0 +3L 18141102 18142039 -1.17015064642504 +3L 18142039 18142143 -1.05910489381442 +3L 18142143 18142231 -4.35793833134167
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/peaks.gff Fri Apr 20 04:34:17 2018 -0400 @@ -0,0 +1,43 @@ +3L . . 17416196 17416656 2.06 . . FDR=0.0434782608695652 +3L . . 17428494 17428708 4.60 . . FDR=0.25 +3L . . 17434875 17435361 2.36 . . FDR=0.0434782608695652 +3L . . 17436191 17436240 1.55 . . FDR=0.0434782608695652 +3L . . 17454392 17455057 1.47 . . FDR=0.0434782608695652 +3L . . 17469263 17469357 1.76 . . FDR=0.25 +3L . . 17474212 17474436 2.67 . . FDR=0.1 +3L . . 17478792 17479470 2.71 . . FDR=0.25 +3L . . 17480437 17480672 4.80 . . FDR=0.1 +3L . . 17484190 17484263 1.55 . . FDR=0.0434782608695652 +3L . . 17493333 17493632 2.10 . . FDR=0.0434782608695652 +3L . . 17509604 17509706 2.15 . . FDR=0.0434782608695652 +3L . . 17548916 17549057 1.90 . . FDR=0.0434782608695652 +3L . . 17629032 17629240 1.43 . . FDR=0.0434782608695652 +3L . . 17637540 17637772 3.90 . . FDR=0.1 +3L . . 17639208 17639274 1.55 . . FDR=0.0434782608695652 +3L . . 17655832 17656401 4.47 . . FDR=0.0555555555555556 +3L . . 17663729 17663765 2.15 . . FDR=0.0434782608695652 +3L . . 17745743 17745838 1.47 . . FDR=0.0434782608695652 +3L . . 17748063 17748387 2.92 . . FDR=0.0434782608695652 +3L . . 17753481 17753917 4.56 . . FDR=0.1 +3L . . 17756873 17757959 1.63 . . FDR=0.0434782608695652 +3L . . 17801230 17802045 1.85 . . FDR=0.0434782608695652 +3L . . 17841131 17842320 2.31 . . FDR=1 +3L . . 17853549 17854148 1.90 . . FDR=0.0555555555555556 +3L . . 17858323 17858952 1.68 . . FDR=0.111111111111111 +3L . . 17859435 17859724 1.64 . . FDR=0.0555555555555556 +3L . . 17861050 17861152 2.02 . . FDR=0.0434782608695652 +3L . . 17869953 17870956 2.56 . . FDR=0.25 +3L . . 17872304 17873275 2.92 . . FDR=0.25 +3L . . 17874546 17875836 3.37 . . FDR=0.25 +3L . . 17910607 17910788 2.74 . . FDR=0.1 +3L . . 17918198 17918586 1.62 . . FDR=0.0434782608695652 +3L . . 17936813 17937334 1.85 . . FDR=0.0434782608695652 +3L . . 17999429 17999479 2.15 . . FDR=0.0434782608695652 +3L . . 18008966 18009131 2.85 . . FDR=0.1 +3L . . 18016765 18016884 1.47 . . FDR=0.0434782608695652 +3L . . 18017196 18017386 2.02 . . FDR=0.0555555555555556 +3L . . 18018167 18018993 1.65 . . FDR=0.0434782608695652 +3L . . 18036026 18036121 4.10 . . FDR=0.0434782608695652 +3L . . 18079650 18079739 2.15 . . FDR=0.0434782608695652 +3L . . 18110932 18111360 1.88 . . FDR=0.1 +3L . . 18131560 18131987 3.45 . . FDR=1