Question: Criteria for merging multiple VCF files
0
gravatar for seta
5 weeks ago by
seta1.2k
Sweden
seta1.2k wrote:

Dear all,

I merged many single VCF files to produce a multi-sample VCF file (human) using bcftools (bcftools merge -m non). Before merging, I split multiallelic sites with bctools, too. However, there are many sites with no genotype (call), just there is ./.. Could you please let me know if shall I apply any filtering step on each single VCF file before merging step?

Thanks

merge bcftool vcf • 111 views
ADD COMMENTlink modified 5 weeks ago by Brice Sarver3.3k • written 5 weeks ago by seta1.2k
0
gravatar for Brice Sarver
5 weeks ago by
Brice Sarver3.3k
United States
Brice Sarver3.3k wrote:

If you're combining non-filtered VCFs, it would make sense to expect a lot of false positive variant calls that are unique to a single sample, resulting in missing genotypes in other samples. You could avoid or reduce this in a number of ways:

  1. See the GATK's HaplotypeCaller documentation on how to generate and merge gVCFs; individuals that have reference calls will correctly have 0/0 when applicable.
  2. Filter the VCFs before so you only have confidently-called variation, then merge them.
  3. Additionally filter on fixed number or percentage of individuals with missing data in your multisample VCF.
  4. Replace missing genotypes with reference calls in your final VCF (has a potential to be problematic).
ADD COMMENTlink modified 5 weeks ago • written 5 weeks ago by Brice Sarver3.3k

Thanks for your comments. Actually, I have VCF files, not gvcf, that the minimum quality score (5th column of VCF file) is 30. I did a basic filtering on single VCF file to keep all variants (SNP and Indel) with the minimum DP of 10 and GQ of 20; then I'm going to merge all filtered single VCF files again followed by keeping just variants present in 80% of individuals within the merged (multi-sample) VCF file. Could you please share me your idea? is this filtering enough? Please kindly tell me if the additional filtering is needed.

Thanks

ADD REPLYlink written 5 weeks ago by seta1.2k
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: 866 users visited in the last hour