call haplotypes in GATK HaplotypeCaller
0
0
Entering edit mode
2.5 years ago
Peter Chung ▴ 190

I am new in bioinformatics. I would like to call haplotypes in GATK HaplotypeCaller, however the output result is 0kb. Can anyone give me some advice how to do the haplotypes ? I am doing the WGS.

REF="refs/ucsc.hg19.fasta"
name="samples/sample-1"
gatk --java-options "-Xmx16g" HaplotypeCaller \
-R $REF \
-I ${name}.addRG.mkdup.recal.bam \
-ERC GVCF \
-O ${name}.g.vcf.gz \
-bamout ${name}.haplotypes.bam \
--bam-writer-type CALLED_HAPLOTYPES \
--do-not-run-physical-phasing false

I am wondering the haplotypes output file. should it be a vcf file or a bam file ? Thanks.

sequencing GATK Haplotype HaplotypeCaller • 851 views
ADD COMMENT
0
Entering edit mode

Hello,

which version of GATK are you running? Are there any warning/error messaging?

fin swimmer

ADD REPLY
0
Entering edit mode

Hi finswimmer,

I am using GATK4, there are not any messages but the output is 0kb. Thanks.

ADD REPLY
0
Entering edit mode

GATK (like most stuff from Broad), is probably the most talky (verbosity) tool I ever came across. There must be something unless you redirected stderr somehow.

ADD REPLY

Login before adding your answer.

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