Visualize Vcf In Igv
when i was trying to use R to visualize my VCF file I first loaded the file and compressed it:
vcf <- readVcf(str.vcf, "Genome")
compressVcf <- bgzip(vcf, tempfile())
but when i run indexTabix
idx <- indexTabix(compressVcf, "vcf")
Error in value[3L] : internal: samtools invoked 'exit(1)'; see warnings() and restart R file: /tmp/RtmpTzJzHz/filea305fdaaade In addition: Warning message: In doTryCatch(return(expr), name, parentenv, handler) : [ti_index_core] the file out of order at line 19
It seems my file is not in order, but what kind of order or how should i sort my vcf or genome ?
Thanks!
• 3,673 views
•
link
1 answer
never mind, the answer seems lie in this post Tabix -p vcf ERROR
• 67 views
•
link
Log in to answer this question.
Do the chromosome names match EXACTLY between the vcf and the genome?
It seems my file is not in order, but what kind of order or how should i sort my vcf or genome ?