https://github.com/gantzgraf/vcfhacks
Raw File
Tip revision: 3b569742cedc2eded1d60bdaba1c2ee91134c8e1 authored by David A. Parry on 10 July 2020, 10:56:40 UTC
Update readme.md
Tip revision: 3b56974
countVariants.pl
#!/usr/bin/env perl
use strict;
use warnings;
use FindBin qw($RealBin);
use lib "$RealBin/lib/dapPerlGenomicLib";
use VcfReader 0.3;

die "Usage: countVariants.pl input1.vcf input2.vcf ...\n" if @ARGV <  1;
my $total = 0;
while (my $vcf = shift){ 
    my $variants = VcfReader::countVariants($vcf) ;
    print "$vcf:\t$variants variants\n";
    $total += $variants;
}
print "Total:\t$total variants\n";
back to top