GATK allele frequencies
1
0
Entering edit mode
5.1 years ago
jyu429 ▴ 120

Hi,

I'm using the GATK Haplotype caller to call for variants, but I'd like to output the genotype fields. For this, I believe I have to use the -GF option, but is this an incorrect way?

java -jar $HOME/Tools/GenomeAnalysisTK-3.6/GenomeAnalysisTK.jar -GF -T HaplotypeCaller -R "imgt_VH_filtered.fasta" -I $PREFIX"realigned_bams/"$INDIVIDUALS_NAME"_realigned.bam" -ploidy 2 -stand_call_conf 30 -stand_emit_conf 10 -o $PREFIX"raw_vcfs/"$INDIVIDUALS_NAME"_raw.vcf"

GATK allele frequency • 1.4k views
ADD COMMENT
0
Entering edit mode
5.1 years ago
vakul.mohanty ▴ 240

I think you will have to carry out joint genotyping to obtain genotype calls.

ADD COMMENT

Login before adding your answer.

Traffic: 2224 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