Hello everyone,
How to compare two VCF files called with the same software (FreeBayes), but with completely different reference genomes based on the same specie?
The question has been asked here already: C: Comparing VCFs constructed from mapping to different reference genomes
But it seems that LiftoverVcf is made to adjust VCF sequence called with similar genomes rather than with completely different reference sequences (with different headers).
Would you have any trick to suggest to compare these VCF files?
Thanks for your answer igor. I was planning to compare a VCF file mapped to a genome to target-capture VCFs mapped to the baits sequence.
Eventually and as you quote, it is easier to map the target-capture dataset to the same reference genome than the first VCF file. Thanks for the help!
It's generally a good idea to map to the full genome, not just the target regions. Your target capture protocol does not capture just the baits.