Bug Bulletin: The recent 3.2 release fixes many issues. If you run into a problem, please try the latest version before posting a bug report, as your problem may already have been solved.

Missing data in HaplotypeCaller

digitonindigitonin Los AngelesPosts: 1Member

I have an issue with HaplotypeCaller and its use to call SNPs in RAD data.

I recently used HaplotypeCaller to call SNPs in 600 samples + 50 subspecies samples. Worked fine. After adding more data to 1200 samples + 120 subspecies samples, ~100 of these subspecies results in 0 calls and just missing data "./." at these loci. Recall that some of these samples and sites were analyzed and called in the first analysis. Any ideas why?

My setup: ~1200 Read Reduced bam files Merged into 1 bam file with proper RG headers. Should I not have merged? Should I not have Reduced Reads?

Thanks!

Tagged:

Answers

Sign In or Register to comment.