Mercurial > repos > elixir-it > covacs_freebayes
annotate filter.fb.pl @ 1:cbe203c9bc3a draft default tip
Uploaded
| author | elixir-it |
|---|---|
| date | Thu, 15 Nov 2018 15:55:26 -0500 |
| parents | 03593410f057 |
| children |
| rev | line source |
|---|---|
| 0 | 1 #!/usr/bin/perl -w |
| 2 # | |
| 3 $f=shift; | |
| 4 $out=shift; | |
| 5 open(IN,$f); | |
| 6 open(OUT,">$out"); | |
| 7 while(<IN>) | |
| 8 { | |
| 9 if ($_=~/^#/) | |
| 10 { | |
| 11 print OUT; | |
| 12 next; | |
| 13 }else{ | |
| 14 $vl=(split())[5]; | |
| 15 $gt=(split())[-1]; | |
| 16 $gt=(split(/\:/,$gt))[0]; | |
| 17 next if $gt eq "0/0"; | |
| 18 print OUT if $vl>20; | |
| 19 } | |
| 20 } | |
| 21 close(OUT); |
