Question: Compare Vcf Files
gravatar for win
6.6 years ago by
win810 wrote:

HI all, I was wondering if someone could help. What would be the best way to compare two VCF files? I was thinking of using BED tools but not sure if there is anything more advanced out there which can output statistics as well.

thanks in advance.

vcf • 4.8k views
ADD COMMENTlink modified 6.6 years ago by Pierre Lindenbaum122k • written 6.6 years ago by win810

There is a program called vcf-compare in vcftools for comparing files. That's probably a good start, though you'll have to be more specific if that is not exactly what you want (vcf-stats may also be of interest).

ADD REPLYlink modified 6.6 years ago • written 6.6 years ago by SES8.2k
gravatar for Ashutosh Pandey
6.6 years ago by
Ashutosh Pandey11k wrote:

Try vcftools ( You will have to use tabix and bgzip that comes with tabix. These tools are basically for compressing and indexing the vcf file. It has to be done beofre using vcf tools.

ADD COMMENTlink written 6.6 years ago by Ashutosh Pandey11k
gravatar for Pierre Lindenbaum
6.6 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum122k wrote:

I wrote a tool named vcf2sqlite: it put one ore more VCFs into a local sqlite3 database. The schema is simple and you can query the database using some simple SQL statements.

$ sqlite3 -column -header  db.sqlite \
   "select,VCFCALL.value,count(*) from VCFCALL,SAMPLE where and prop='GT' group by,VCFCALL.value"

name         value       count(*)  
-----------  ----------  ----------
rmdup_1.bam  0/1         545       
rmdup_1.bam  1/1         429       
rmdup_2.bam  0/1         625       
rmdup_2.bam  1/1         349       
rmdup_3.bam  0/1         595       
rmdup_3.bam  1/1         379       
rmdup_4.bam  0/1         548       
rmdup_4.bam  1/1         426       
rmdup_5.bam  0/1         564       
rmdup_5.bam  1/1         410       
rmdup_6.bam  0/1         724       
rmdup_6.bam  1/1         250
ADD COMMENTlink written 6.6 years ago by Pierre Lindenbaum122k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1161 users visited in the last hour