Tagged with #bamout
0 documentation articles | 0 announcements | 1 forum discussion


No posts found with the requested search criteria.
No posts found with the requested search criteria.
Comments (0)

I learn that HaplotypeCaller looks for candidate regions and builds a denovo assembly. So how is this region defined? How big is the region?

Recently we found a six base pair insertion in an interesting gene using HaplotypeCaller which UnifiedGenotyper missed. We sanger sequenced it and found it to be true. Happy!!! But when I tried to write a bam file with bamout option using a 1MB flanking window and the alignment looked totally new and the variant disappeared. In fact there were no reads covering that position in the new assembly.