Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Forum
Planet
All »
Home
Log In
Question: Alternative To Fpkm From Cufflinks
0
gravatar for James Chitwood
7.4 years ago by
James Chitwood • 10
James Chitwood • 10 wrote:
Hi all, Is there any way to find out the number of reads aligning to a transcript rather than the FPKM calculated by Cufflinks? I'm also interested in obtaining summary statistics for mapping analyses or RNA-Seq data, such as % of reads aligned, % uniquely aligned, mapped to exons, introns, etc. Is there a tool that would provide a summary table with this information? Thanks for your help and providing this great resource, James
rna-seq cufflinks • 1.1k views
ADD COMMENT • link •
modified 7.4 years ago by Jennifer Hillman Jackson ♦ 25k • written 7.4 years ago by James Chitwood • 10
0
gravatar for Jennifer Hillman Jackson
7.4 years ago by
Jennifer Hillman Jackson ♦ 25k
United States
Jennifer Hillman Jackson ♦ 25k wrote:
Hello James, For coverage, a rough coverage attribute is available in the Cufflinks' GTF file. Absolute coverage is not calculated because RNA-seq reads can map to multiple isoforms and it's not possible to say which one it truly maps to (hence the need for Cufflinks' statistical models). To address some of the other statistics, the tool "NGS: SAM Tools -> flagstat" provides some basic counts for BAM files. For the rest, you will need to use a combination of the tools available for file conversions and data manipulation found in tool groups such as "Text Manipulation", "Filter and Sort", "Join, Subtract and Group", and the "NGS*" groups. As an example, to identify the reads mapping to particular elements: convert Tophat's mapped reads to intervals (BAM-to-SAM, then SAM-to-interval) and intersect those coordinates with subsets of a GFF file (created using GFF filtering by keyword, e.g. "exon"). Then use the Group tool to generate counts. Most sorts of manipulations one would do with simple scripting to gather counts are possible. Once you work out how to generate a particular statistic, be sure to save the steps into a workflow for future use (and consider publishing it for others to use). Hopefully this helps! Best, Jen Galaxy team -- Jennifer Jackson http://usegalaxy.org/ http://galaxyproject.org/
ADD COMMENT • link written 7.4 years ago by Jennifer Hillman Jackson ♦ 25k
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 173 users visited in the last hour