Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #very simple vcf file filter
- #usage:
- #perl6 simple_filter.p6 filtered.vcf stock.vcf
- #VR
- my @info;
- #output:
- my $fh = open @*ARGS[0], :w;
- for @*ARGS[1].IO.lines -> $line {
- if $line ~~ /^\#/ {
- $fh.print("$line\n");
- } else {
- @info = split(':', $line);
- #all you need is proper number here instead of 7 to catch investigated values and comparison
- $fh.print("$line\n") if @info[7] > 50;
- }
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement