Hi,
Thank u for your reply
The chromosome names are the same in both reference and gtf file: "chromosome".
The out put of running IDxstats on my bam file was:
1 2 3 4
chromosome 6537648 1021248 0
- 0 0 2584368
I am pretty sure that I have mapped the reads since I had previously used the result of mapping in cuffdiff and I got the results.
Since I have to remove the batch effects between my samples I have to run deseq2 and I needs the raw read counts
Have you had a look at your .bam file in IGV? This will allow you to visualise the alignments from Bowtie and indicate if the problem lies at your mapping stage.