Using Tophat "accepted_hits.bam" file for RNAseq variant calling with GATK
Hi, I have some doubts about variant calling in RNAseq data.
So far, I have aligned the reads to a reference (human hg38) genome using Tophat. However, I have seen in your workflow for variant discovery in RNAseq that you recommend to use STAR. Is this really necessary? Could I use the "accepted_hits.bam" file from Tophat instead of the STAR output? would the rest of the workflow remain the same?
Thanks for your help