Visualizing VCF file on IGV

Hi all,
This is my first time trying to understand and visualize a VCF file on IGV, I make my vcf as follows:

samtools mpileup -I -u -f genome.fa alignment.bam > pileup.bcf
bcftools call -cv pileup.bcf > variant.bcf

Everything worked, the final result looks like the manual said it should be, but IGV says:

Error loading /variant.bcf: Unable to parse header with error: Input stream does not contain a BCF encoded file; BCF magic header info not found, at record 0 with position 0:, for input source: ./mysource

Any suggestion?


BCFTOOLS


VCF


IGV

• 84 views



Source link