We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

HaplotypeCaller and vcf file

Hi, I made the haplotype caller but most of my file looks like that:

CHROM POS ID REF ALT QUAL FILTER INFO FORMAT L001

chr1 1 . N . . END=10055 GT:DP:GQ:MIN_DP:PL 0/0:0:0:0:0,0,0
chr1 10056 . A . . END=10074 GT:DP:GQ:MIN_DP:PL 0/0:4:6:3:0,6,90
chr1 10075 . A . . END=10077 GT:DP:GQ:MIN_DP:PL 0/0:3:9:3:0,9,81
chr1 10078 . C . . END=10078 GT:DP:GQ:MIN_DP:PL 0/0:3:9:3:0,9,85
chr1 10079 . T . . END=10082 GT:DP:GQ:MIN_DP:PL 0/0:3:9:3:0,9,82

with no filter and no QUAL number. I have a few lines like that:
chr1 14604 . A G, 20 . DP=1;MLEAC=1,0;MLEAF=0.500,0.00;MQ=21.00 GT:AD:DP:GQ:PGT:PID:PL:SB 1/1:0,1,0:1:3:0|1:14590_G_A:45,3,0,45,3,45:0,0,1,0

Is that normal?

I run the haplotypeCaller as follows:
java -jar /usr/local/Modules/modulefiles/tools/gatk/3.4.0/GenomeAnalysisTK.jar -T HaplotypeCaller -R reference.fa -I BaseRecalibrated.bam --emitRefConfidence GVCF --genotyping_mode DISCOVERY --dbsnp dbsnp.vcf -o out.g.vcf

Best Answer

Answers

Sign In or Register to comment.