Mercurial > repos > portiahollyoak > temp
comparison scripts/filterFalsePositive.ex.pl @ 0:28d1a6f8143f draft
planemo upload for repository https://github.com/portiahollyoak/Tools commit 132bb96bba8e7aed66a102ed93b7744f36d10d37-dirty
| author | portiahollyoak |
|---|---|
| date | Mon, 25 Apr 2016 13:08:56 -0400 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:28d1a6f8143f |
|---|---|
| 1 #! /usr/bin/perl | |
| 2 | |
| 3 use strict; | |
| 4 | |
| 5 open (input, "<$ARGV[0]") or die "Can't open $ARGV[0] since $!\n"; | |
| 6 my %leng=(); | |
| 7 my %trans=(); | |
| 8 my %coordinate=(); | |
| 9 while (my $line=<input>) { | |
| 10 chomp($line); | |
| 11 my @a=split(/\t/, $line); | |
| 12 if (defined $leng{$a[9]}) { | |
| 13 $trans{$a[9]} += $a[2]-$a[1]; | |
| 14 } | |
| 15 else { | |
| 16 $leng{$a[9]}=$a[8]-$a[7]-10; | |
| 17 $trans{$a[9]}=$a[2]-$a[1]; | |
| 18 $coordinate{$a[9]}="$a[6]\:$a[7]\:$a[8]"; | |
| 19 } | |
| 20 } | |
| 21 close input; | |
| 22 | |
| 23 open (output, ">>$ARGV[2]") or die "Can't open $ARGV[2] since $!\n"; | |
| 24 while ((my $key, my $value) = each (%coordinate)) { | |
| 25 if ((($leng{$key}-$trans{$key}) <= $ARGV[1])&&(($leng{$key}-$trans{$key}) >= 0)) { | |
| 26 # if (($leng{$key}-$trans{$key}) <= 500) { | |
| 27 my @b=split(/\:/, $value); | |
| 28 print output "$b[0]\t$b[1]\t$b[2]\t$key\n"; | |
| 29 } | |
| 30 } | |
| 31 close output; |
