hello,
I wonder if anybody else has come across this issue or has some insight into the same. I have RNA-seq data and I usually align to the genome using GSNAP or Tophat. Thus far I have found Cufflinks to be the only major tool to estimate Isoform FPKMs using this approach (I know there are a couple of other tools out there). For comparison, I wanted to estimates FPKMs using RSEM, eXPress or BitSeq with the same bams, but am stuck as they require the alignments to be made to the reference transcriptome not the genome?? Any suggestions to circumvent this issue other than complete realignment??. Some thoughts along this line are :
Any thoughts, insights in this context will be much appreciated
cheers
Ashok
I wonder if anybody else has come across this issue or has some insight into the same. I have RNA-seq data and I usually align to the genome using GSNAP or Tophat. Thus far I have found Cufflinks to be the only major tool to estimate Isoform FPKMs using this approach (I know there are a couple of other tools out there). For comparison, I wanted to estimates FPKMs using RSEM, eXPress or BitSeq with the same bams, but am stuck as they require the alignments to be made to the reference transcriptome not the genome?? Any suggestions to circumvent this issue other than complete realignment??. Some thoughts along this line are :
- Is there a tool to convert genome space to transcriptome space and vice versa
- would such a tool be useful
- what would be needed in this context?
- sam flags when reads map to the same exons in two isoforms?
- will realignment be needed?
Any thoughts, insights in this context will be much appreciated
cheers
Ashok
Comment