Question: Comparing samtools,gatk and varscan against gold standard 1000G snps and need to decide the best method
0
gravatar for bioinforesearchquestions
2.3 years ago by
United States
bioinforesearchquestions160 wrote:

Dear All,

For three samples, I generated vcf files using samtools, gatk, and varscan.

I tried accessing the files 1000G_phase1_high_confidence_snps.vcf, it has only following columns 1) Chrom, 2) Pos, 3) ID, 4) Ref, 5) Alt, 6) Qual, 7) Filter and 8) Info.

I have vcf files with following columns samtools: 1) Chrom, 2) Pos, 3) ID, 4) Ref, 5) Alt, 6) Qual, 7) Filter, 8) Info and 9) Sample1. varscan: 1) Chrom, 2) Pos, 3) ID, 4) Ref, 5) Alt, 6) Qual, 7) Filter, 8) Info and 9) Sample1. gatk:1) Chrom, 2) Pos, 3) ID, 4) Ref, 5) Alt, 6) Qual, 7) Filter, 8) Info and 9) Sample1.

Now I want to compare those vcf files generated by samtools, gatk and varscan against with 1000G_phase1_high_confidence_snps.vcf to find out which variant caller gives the same list of high confident mutations from 100G_phase1?

vcfs snps dnaseq • 1.3k views
ADD COMMENTlink modified 2.2 years ago by Biostar ♦♦ 20 • written 2.3 years ago by bioinforesearchquestions160
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: 687 users visited in the last hour