Question: Compare Vcf Files
0
gravatar for win
5.5 years ago by
win780
India
win780 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.3k views
ADD COMMENTlink modified 5.5 years ago by Pierre Lindenbaum110k • written 5.5 years ago by win780
1

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 5.5 years ago • written 5.5 years ago by SES8.1k
3
gravatar for Ashutosh Pandey
5.5 years ago by
Philadelphia
Ashutosh Pandey11k wrote:

Try vcftools (http://vcftools.sourceforge.net/). 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 5.5 years ago by Ashutosh Pandey11k
3
gravatar for Pierre Lindenbaum
5.5 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum110k 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 SAMPLE.name,VCFCALL.value,count(*) from VCFCALL,SAMPLE where SAMPLE.id=VCFCALL.sample_id and prop='GT' group by SAMPLE.id,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 5.5 years ago by Pierre Lindenbaum110k
Please log in to add an answer.

Help
Access

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