SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Illumina reads failing to map to bacterial genome with Tophat C9r1y Bioinformatics 2 12-07-2015 06:51 AM
does tophat rescue PE reads that do not map by aligning a single end of the PE? zorph Bioinformatics 0 09-09-2011 02:14 PM
Tophat - Left Kept Reads jkozubek Bioinformatics 2 07-27-2011 05:40 AM
Reads map in Tophat, but Cufflinks doesn't find transcripts janec Bioinformatics 0 07-07-2011 01:42 AM
Best tool to map 454 reads onto sanger reads? dan Bioinformatics 3 07-27-2009 08:51 AM

Reply
 
Thread Tools
Old 08-14-2013, 07:13 AM   #1
airad22
Junior Member
 
Location: Europe

Join Date: Aug 2013
Posts: 2
Unhappy Tophat RNASeq mapping - right reads map and left reads map 50% less

Hi!
I am trying to map RNA Seq (chick) with Tophat2 and I am getting low overall mapping rates (~30%). In the align_summary.txt file generated by Tophat2 - it appears that my right reads are mapping at over 90% whereas the left reads are mapping at 30%. This is happening to all but one sample.
Has anyone ever encountered this or know how to resolve this?

This is the command I used:

tophat -p 8 -G gene_model.gtf -o Sample1 genome4 Sample1_1.fastq.gz Sample1_2.fastq.gz > Sample1.log 2>&1 &

Thanks in advance

Adi
airad22 is offline   Reply With Quote
Old 08-14-2013, 07:17 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

Did you run fastQC (or any other QC software) on these reads?
dpryan is offline   Reply With Quote
Old 08-14-2013, 07:26 AM   #3
airad22
Junior Member
 
Location: Europe

Join Date: Aug 2013
Posts: 2
Default

Quote:
Originally Posted by dpryan View Post
Did you run fastQC (or any other QC software) on these reads?
Yes I ran FastQC. The Per base sequence quality is in the green zone for all of them. In some the right reads are slightly lower than the left, but not much.

The only other thing I notice with FastQC is that 'Sequence Duplication Levels' are very high (>=55.7%) . I don't know whether this could influence mapping though.

Thanks.
airad22 is offline   Reply With Quote
Reply

Tags
rnaseq, tophat

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:13 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO