Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Forum
Planet
All »
Home
Log In
Question: Mapping Paired End Reads With Tophat
0
gravatar for 杨继文
6.6 years ago by
杨继文 • 210
杨继文 • 210 wrote:
Hi all, After mapping, I used IGV to have a look at the mapping. There are a lot of mapped reads without pair reads. Should I keep these reads? or Is this a problem for cufflinks analysis? What I tried is: 1. BAM to SAM 2. Filter SAM: set Read mapped in a proper pair to Yes. The result is that only 1/5 reads were left. Can anybody tell me if this operation is proper?? How do you normally optimize the mapping rerults from Tophat? Which considerations should I take into account? Looking forward to your reply. Jiwen
rna-seq cufflinks • 2.3k views
ADD COMMENT • link •
modified 6.5 years ago by Jennifer Hillman Jackson ♦ 25k • written 6.6 years ago by 杨继文 • 210
0
gravatar for Jennifer Hillman Jackson
6.5 years ago by
Jennifer Hillman Jackson ♦ 25k
United States
Jennifer Hillman Jackson ♦ 25k wrote:
Hi Jiwen, The short answer is that if your data is not pairing, then there may be a quality problem. Or, there may be a problem with the TopHat mapping run. The best advice is to take a sample of your data and experiment with some TopHat alternate parameters (using the protocol suggestions in paper below or manual http://tophat.cbcb.umd.edu/) and see what works. If your overall goal is simply transcript/gene discovery/assembly, then filtering is probably OK. But if you are going to be doing any statistical expression analysis, then targeted filtering of the data (e.g. beyond general quality) should be done with caution, if at all, as you risk skewing the results. You may have seen this already, but the Cufflinks tool authors put out a new paper that covers best practice RNA-seq protocols: http://www.nature.com/nprot/journal/v7/n3/full/nprot.2012.016.html (also linked from http://cufflinks.cbcb.umd.edu/, 2nd item down) Apologies for the delayed reply. There were a few questions from you around this same time, but it wasn't clear if everything was addressed or not. And I don't think the paper link was sent out in reply, which will likely be the most helpful. Best, Jen Galaxy -- Jennifer Jackson http://galaxyproject.org
ADD COMMENT • link written 6.5 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: 179 users visited in the last hour