Using an other reference genome
Opened this issue · 0 comments
Enorya commented
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