Hi, all!
I've tried to apply standard Vcf.pm example to SNP data file from 1000 genomes project:
use Vcf;
my $vcf = Vcf->new(file=>'ALL.2of4intersection.20100804.genotypes.vcf.gz',region=>'1:1000-2000');
$vcf->parse_header();
The file is bgzipped and tabix'ed.
And there is an error:
"Broken VCF header, no column names?"
With another small vcf file all is ok.
What is wrong?
Thank you!
SN
I've tried to apply standard Vcf.pm example to SNP data file from 1000 genomes project:
use Vcf;
my $vcf = Vcf->new(file=>'ALL.2of4intersection.20100804.genotypes.vcf.gz',region=>'1:1000-2000');
$vcf->parse_header();
The file is bgzipped and tabix'ed.
And there is an error:
"Broken VCF header, no column names?"
With another small vcf file all is ok.
What is wrong?
Thank you!
SN
Comment