Question: Compare Vcf Files
0
gravatar for win
5.2 years ago by
win730
India
win730 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.2k views
ADD COMMENTlink modified 5.2 years ago by Pierre Lindenbaum106k • written 5.2 years ago by win730
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.2 years ago • written 5.2 years ago by SES8.1k
3
gravatar for Ashutosh Pandey
5.2 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.2 years ago by Ashutosh Pandey11k
3
gravatar for Pierre Lindenbaum
5.2 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum106k 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.2 years ago by Pierre Lindenbaum106k
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: 1031 users visited in the last hour