Dear all,

I have been given to analyse a dataset like this, where I have the counts for each transcript. Since I want to use DESeq I need to find a way to convert the transcript counts to gene counts. The out put comes from HOMER, and I am not sure if tximport in R supports it(?)

# A tidytable: 54,013 × 9
   transcript             chr   start   end strand Length Copies `Annotation/Divergence` counts
   <chr>                  <chr> <int> <int> <chr>   <dbl>  <int>                   <dbl>  <dbl>
 1 transcript:AT1G01010.1 1      3631  5899 +        2268      1                       0      4
 2 transcript:AT1G01020.1 1      6788  9130 -        2342      1                       0     51
 3 transcript:AT1G01020.2 1      6788  8737 -        1949      1                       0     50
 4 transcript:AT1G01020.3 1      6788  9130 -        2342      1                       0     51
 5 transcript:AT1G01020.4 1      6788  9130 -        2342      1                       0     51
 6 transcript:AT1G01020.5 1      6788  9130 -        2342      1                       0     51
 7 transcript:AT1G01020.6 1      6788  8737 -        1949      1                       0     50
 8 transcript:AT1G01030.1 1     11649 13714 -        2065      1                       0     39

Any idea or suggestion is appreciated!
Would it be just enough to sum the transcripts of each in a new column?

Homer is not really standard software for RNA-seq afaik. The easiest would be to use something like salmon to quantify your reads and then use tximport as in the manual. Then you can be sure data are processed correctly in terms of how tximport expects them.