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!

Emitting ref vs. alt allele counts on all sites

Hi,

I would like to use UnifiedGenotyper to produce a vcf file where all sites I specify via BED file will have a count of reads carrying the reference allele and a count of reads carrying an alternate allele, regardless of whether a variant has been called in the site.
I've tried to use EMIT_ALL_SITES, but this appears not to provide this count breakdown for sites that are not called as variable.
I would appreciate any pointers. I am using GATK 2.2.4.

Thank you very much.

Best Answer

Answers

Sign In or Register to comment.