Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Forum
Planet
All »
Home
Log In
Question: Table With Gene Count Reads
0
gravatar for Luciano Cosme
6.7 years ago by
Luciano Cosme • 220
Luciano Cosme • 220 wrote:
Hi, I was wondering if there is any tool on Galaxy were I can obtain a table with how many reads have been mapped to a given sample and to a given gene (for example, use a Tophat output and use a GFF file to obtain the table). I am using HTSeq to get it (htseq-count). There is also GenomicRanges and easyRNASeq packages in bioconductor. Thank you. Luciano
gff • 1.7k views
ADD COMMENT • link •
modified 6.7 years ago by Jennifer Hillman Jackson ♦ 25k • written 6.7 years ago by Luciano Cosme • 220
0
gravatar for Jennifer Hillman Jackson
6.7 years ago by
Jennifer Hillman Jackson ♦ 25k
United States
Jennifer Hillman Jackson ♦ 25k wrote:
Hello Luciano, There is no single tool do to this operation (although there has been some discussion about including one in the Tool Shed), but the same information can be obtained by using a combination of existing tools. First, start by converting both starting datasets to interval format. mapped reads: - for TopHat output, "NGS: SAM Tools -> Convert SAM to interval" features: - for GFF file (convert to tabular if necessary), subtract "1" from the start position's value using tool "Text Manipulation -> Compute" - cut columns chrom, new start, stop, strand, name, and score from this result file using "Text Manipulation -> Cut" - set the data type to "interval" using the 'Edit attributes form (pencil icon) Next, use a tool in the group "Operate on Genomic Intervals" to compare these intervals for overlap. The tool "Cluster" with the option "Find" is mostly likely the one you will want to use. As a final step, summarize the data by feature using the tool "Join, Subtract and Group -> Group". Hopefully this helps, Best, Jen Galaxy team
ADD COMMENT • link written 6.7 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: 172 users visited in the last hour