I am hoping to run an allele specific expression analysis on a set of RNAseq samples I have. I need to obtain the genotypes for all samples to determine heterozygosity of each variant which is needed for my final analysis.
How can I obtain the variant genotypes per sampple?
I have tried GATK HaplotypeCaller on my filtered .bams however this takes an immense amount of time on the HPC I have access to. Is there another method besides variant calling that can produce genotypes, or any way reduce the compute time for the variant calling.
Any help would be immensly appreciated.
Note: I am new to this kind of analysis, so please excuse my lack of knowledge of some terminology.