cccnrc/plot-VCF

Using an other reference genome

Opened this issue · 0 comments

Dear,

I'm trying to use your tool to visualize Structural variants but I'm working with an other organism than human so I would like to change the reference genome to use.
I tried to do the following:

vcf_file <- './comparison_1kbp_typesafe_min3_30bp_strand-type1.vcf'
chrom=c("NC_058021.1","NC_058022.1","NC_058023.1","NC_058024.1","NC_058025.1",
        "NC_058026.1","NC_058027.1","NC_058028.1","NC_058029.1","NC_058030.1",
        "NC_058031.1","NC_058032.1","NC_058033.1","NC_058034.1","NC_058035.1",
        "NC_058036.1","NC_058037.1","NC_058041.1","NC_058038.1","NC_058039.1",
        "NC_058040.1")
createVCFplot(vcf_file, FASTA_FILE="./Solea_senegalensis/ncbi_dataset/data/GCF_019176455.1/GCF_019176455.1_IFAPA_SoseM_1_genomic.fna", ASSEMBLY="SoSeM",  CHR_NAMES=chrom)

But I end up with the following error message:

Erreur dans .Call2("new_input_filexp", filepath, PACKAGE = "XVector") : 
  cannot open file './Solea_senegalensis/ncbi_dataset/data/GCF_019176455.1/GCF_019176455.1_IFAPA_SoseM_1_genomic.fna'

I'm sure the path to the file is correct and the file is for sure in fasta format because it was downloaded directly from NCBI.

Can you help me with this issue?

Thank you in advance,
Enora