Handling single-sample VCF after haplotypecaller
0
0
Entering edit mode
13 months ago
Y.H.P • 0

Hi,

I'm working on QTL mapping of diploid plants and succeeded in first GATK haplotypecaller run.

Results are about 300 VCF files, each containing a single sample of the plant population.

I found out that hard-filtering and BQSR are recommended after first haplotypecaller run, and thus I'm planning to run VariantFiltration command line of GATK.

My question is that, if each raw VCF file is per-sample file, is it recommended to run further procedures by merging all files into a single file, or do I have to go on with separate single-sample files?

Additionally, now every file has a different header and will it be helpful to go on for merged, multi-sample single VCF file for further studies?

Waiting for help,

Thank you.

VCF • 681 views
ADD COMMENT
1
Entering edit mode

300 VCF files should be called using the GVCF mode followed by CombineGVCFs and GenotypeGVCFs

ADD REPLY
0
Entering edit mode

Does it mean that before going through combining and final genotyping, all procedures regarding handling VCF files (ex: hard-filtering, BQSR) should be conducted with separate single-sample file, using -ERC GVCF option in every procedure?

ADD REPLY

Login before adding your answer.

Traffic: 2657 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6