What is a tool to get the genome build of a VCF?
It shouldn't be too hard to create one, but if one exists already that's even better. I need it to be automatable / non-web based (assume no relevant info exists in the header).
• 1,318 views
•
link
updated
by
Ram
•
written
by
a5864557 •
•
written
by
a5864557 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Non-zero expression counts from deleted genes
written by loughrae
-
How do I put the color scales on heatmap.2?
written by Bioinf_Questions •
-
Error: mysql table already exists! Is it safe to delete the tables in mysql manually?
written by mthm •
-
Better place to download dbSNP than official ftp site?
written by MatthewP
-
How to calculate VAF from BAM or SAM(+pileups) files using pysam?
written by askif4 •
-
Combine multiple hard drives into one big volume in linux
written by ASid •
-
Is there an easier way to collect information from Cantatadb target interacting transcripts?
written by Arsenal •
-
visulaization of several strings
written by Learner
-
Integrating Genome Tracks Inside A Web Page
written by Abhi
-
Plotting Selection On Structures
written by Arjen Ten Have
I don't think you can get genome build unless it's written in VCF header. Maybe you can guess according to some contig IDs.
First, I would never guess on an analysis. If there is no code and documentation available for a file then I'd never use it. But if you are absolutely forced to, you might get positions for common variants, dbSNP for example with high AF, for hg18, hg19, hg38, and then intersect these files with your VCF. The correct build should have the best overlap with the VCF on these sites.
@MatthewP ATpoint It is definitely possible to deterministically determine the build, by cross-referencing dbSNP (e.g) as ATpoint mentioned. If you do this for a large number of both b37 and b38 matches and you get 1000 matches on b38 and none on b37, we can assign the build to b38. Unless I'm missing something?