Question: Tophat
gravatar for 杨继文
6.7 years ago by
杨继文210 wrote:
Dear all, This might be a silly question, but I couldn't figure it out by myself :-((. Could you please tell me how I can find out how many reads have been mapped to the genome after running Tophat for pair end RNA seq data? Thanks in advance. JIwen
rna-seq tophat • 770 views
ADD COMMENTlink modified 2.3 years ago by y.hoogstrate460 • written 6.7 years ago by 杨继文210
gravatar for Jennifer Hillman Jackson
6.7 years ago by
United States
Jennifer Hillman Jackson25k wrote:
Hello Jiwen, The tool "NGS: Picard (beta) -> SAM/BAM Alignment Summary Metrics" gives a nice set of statistics for paired end data. Hopefully this helps, Best, Jen Galaxy team -- Jennifer Jackson
ADD COMMENTlink written 6.7 years ago by Jennifer Hillman Jackson25k
gravatar for y.hoogstrate
2.3 years ago by
y.hoogstrate460 wrote:

I don't know how Tophat in Galaxy does this currently, but in CLI Tophat, unmapped reads were put in a separate BAM file making it hard/impossible to obtain statistics that involve unmapped reads.

Tophat in Galaxy should also produce a summary file that gives some basic statistics including the percentage mapped/unmapped reads. The file should appear in your galaxy history as e.g. "100 Tophat on data 16: align_summary"

ADD COMMENTlink written 2.3 years ago by y.hoogstrate460
Please log in to add an answer.


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