Hello,

I am currently trying to call a consensus sequence for a genome that is 1.2 Gbp in length.

I have a VCF file that has only ~7 million records, including variant and invariant sites.

When I run bcftools consensus Test.vcf.gz -f Reference.fasta -a N;

I get an output consensus with ONLY 5% N content.

This is confusing to me because I wouldn't expect 7 million records to contain enough information to get a consensus for a genome of this size.

I have confirmed on a test file that consensus is working as expected and not including reference bases that are not covered in the VCF file.

Any idea what could be going on?