我从gnomAD数据库下载了人群频率的数据,是一个bgz压缩的文件,“gnomad.exomes.r2.1.1.sites.10.liftover_grch38.vcf.bgz”,然后想解压提取亚洲人群基因变异频率数据,但是在提取的时候出现了报错
我用的代码如下
zcat gnomad.exomes.r2.1.1.sites.10.liftover_grch38.vcf.bgz | grep -v "##"| grep -c "chr10" | head -n 5
运行结果及报错内容
报错如下:
gzip: gnomad.exomes.r2.1.1.sites.10.liftover_grch38.vcf.bgz: invalid compressed data--crc error
gzip: gnomad.exomes.r2.1.1.sites.10.liftover_grch38.vcf.bgz: invalid compressed data--length error
134809
请问一下是什么原因?该怎么解决才能提取到这些变异位点的频率?