From bam files to isoform read counts
2
0
Entering edit mode
7.5 years ago
ddzhangzz ▴ 90

I downloaded RNASeq bam files from TCGA data portal. I am interested in the expression analysis at isofrom level and wondering how to translate these bam files into counts at each isoform. I have done some analysis at gene level using star or htseq-count but fresh to isoform analysis. Could somebody have me some suggestions. I knew there are softwares such as rMATs for alternative splicing (AS) analysis but what I wanted are the counts for each isoform rather than the AS analysis.

RNA-Seq • 2.8k views
ADD COMMENT
0
Entering edit mode

You would want to carefully check how multi-mappers were handled if you are starting with BAM files and considering isoforms.

ADD REPLY
0
Entering edit mode
7.5 years ago
Benn 8.3k

There is an example (case study) in the edgeR user's guide. Chapter 4.5. They use featureCounts first instead of htseq.

https://bioconductor.org/packages/release/bioc/vignettes/edgeR/inst/doc/edgeRUsersGuide.pdf

ADD COMMENT
0
Entering edit mode
7.1 years ago
firatuyulur ▴ 320

Try sailfish They are not the exact counts but would work for many purposes.

ADD COMMENT

Login before adding your answer.

Traffic: 1402 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6