https://github.com/gantzgraf/vcfhacks
Tip revision: 3b569742cedc2eded1d60bdaba1c2ee91134c8e1 authored by David A. Parry on 10 July 2020, 10:56:40 UTC
Update readme.md
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";