I am trying to incorporate variants in my ref.fa file using bcftools consensus. But the output consensus.fa has only one chromosome (the first one from the ref.fa file). Why are the other chromosomes not showing up in consensus .fa?
The command used was :
bcftools consensus -f hg38.fa compressed.vcf.gz > consensus.fa
compressed.vcf.gz is my variants file.