We are trying to sort out a pipeline for doing allele specific expression. Our plan is to call SNPs from RNA-seq data and combine with known SNP annotations. A well known problem in ASE is reference bias, where reads are more like to map if they have the reference sequence than the alt sequence.
We have come accross two possible ways to deal with this. In many papers (e.g. Zhuo et al, Sigurdsson et al) a new reference is generated, substitute N bases where there are SNPs using FastaAlternativeReferenceMaker from GATK and reads are remapped to that.
HISAT2 offers the possibility to provide SNP calls to the index builder, and build an index that accounts for these SNPs.
Does anyone have any advice on which of these might be better?