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!
GATK reports wrong allele
I downloaded the following BAM file: ftp://ftp.sra.ebi.ac.uk/vol1/ERA596/ERA596386/bam/Burkhardtshohle.bam
I sorted and indexed the file with samtools.
Then I created a VCF with the following command:
java -Xmx4g -jar GenomeAnalysisTK.jar -l INFO -R resources/Homo_sapiens_assembly19.fasta -T UnifiedGenotyper -I Burkhardtshohle_sorted.bam -L Y -rf BadCigar -o Burkhardtshohle.vcf --output_mode EMIT_ALL_SITES -nct 4
At position Y:7190277 it reports reference allele A and no alt allele. However, the BAM file has G at this position. I also ran this with HaplotypeCaller and similar options but it still gave me A instead of G.
Am I doing something wrong, or...?