Converting Thousand genome BCF file into VCF with BCFtools

Question

[always_learning]

Hi All.

I am trying to convert 1000 Genome BCF files into VCF files using bcftools and process is getting quit automatically after just converting 2-3 chromosomes. Just wondering if some one else is getting same problem as well ?

bcftools view ALL.wgs.phase3_shapeit2_mvncall_integrated_v5.20130502.genotypes.bcf | bgzip -c > ALL.wgs.phase3_shapeit2_mvncall_integrated_v5.20130502.vcf.gz

If problem is due to bgzip or bcftools ?

This one is version of BCFtools that I am using. bcftools 1.3 Using htslib 1.3

[dschika]

Is there an error message? How large is your bcf?

Answer 1

[Jorge Amigo]

instead of piping the output to bgzip, try creating the output file directly:

bcftools view -Oz \
-o ALL.wgs.phase3_shapeit2_mvncall_integrated_v5.20130502.vcf.gz \
ALL.wgs.phase3_shapeit2_mvncall_integrated_v5.20130502.genotypes.bcf

you do know you have the latest .vcf files already available, don't you?

[always_learning]

Yes I know but since I had .bcf file already to trying to convert that only. Thanks