I am a master student in Bioinformatics and this is my first post in Biostars, so please be gentle. :)
I started using vg as part of a project I have for the university, and already got into troubles. I must say here, that i have never used vg before and I am just trying to make it work and learn more about it.
So, I downloaded the version 38 for the human reference genome along with common_all.vcf.gz in order to construct my graph. I created a index file for my VCF file and then i executed the following command:
vg construct -r Homo_sapiens.GRCh38.dna.alt.fa -v 00-common_all.vcf.gz > reference_grap.vg
It runs for a while, trying to create a index file for the reference genome and then it is crashing, giving the following error:
vg: src/constructor.cpp:2095: void vg::Constructor::construct_graph(const std::vector<FastaReference*>&, const std::vector<vcflib::VariantCallFile*>&, const std::vector<FastaReference*>&, std::function<void(vg:$id(vg::Graph&)>): Assertion `reference_for.count(fasta_contig)' failed
I tried with different reference genomes and different VCF files (00-All.vcf.gz) but again im getting the same error.
I would appreciate any help regarding this matter,
Thank you, Nikiforos