RNA-seq analysis, wrong genome build

RNA-seq analysis, wrong genome build


Hi, I have been doing a feature count and subsequent differential gene expression analysis on some RNA-seq samples which I now suspect is giving me poor results because I used a GRCm39 feature file from Ensembl but bam files which I suspect were aligned to mm10 (based on examination of bw files I produced from the bam files I was given from the sequencing facility).

Is there an Ensembl repository of feature files based on previous genome builds that I can access? I can’t seem to find any on the Ensembl website! Alternatively, I believe I will have to re-align my fastq files to the new genome build, because as far as I can see, the Ensembl assembly converter tool doesn’t support .bam files. Or is it ‘safe’ to just convert my gtf feature file from mm39 to mm10 using this Ensembl assembly converter, and then re-run the feature count? What would you do? Any input on how to best approach this will be much appreciated!






Read more here: Source link