SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
why low mapping rates for RNAseq? NGSfan RNA Sequencing 49 09-09-2014 12:34 PM
reason for low mapping rate?? miaom RNA Sequencing 3 05-10-2014 07:25 AM
The low mapping rate vivienne_lovely Bioinformatics 7 06-05-2013 05:45 PM
TopHat extremely low paired mapping rate. PLS HELP! Felix.Lee RNA Sequencing 0 05-22-2013 01:00 PM
General mapping rate of human resequencing data against reference in GAiix/Hiseq cybog337 Illumina/Solexa 2 01-12-2011 08:43 AM

Reply
 
Thread Tools
Old 08-30-2013, 06:55 PM   #1
Markovia
Junior Member
 
Location: United States

Join Date: Aug 2013
Posts: 1
Post Low mapping rate for RNAseq 2x150 trimmed data

I just got a data set from my collaborator. 2 X 150 bps, Illumina RNAseq data for human samples. We did QC on the data and trimmed it to remove adaptors and low quality bases. Then I use tophat2 to map them to hg19 in three ways:

1. paired end mapping with default parameters;
2. single end mapping for each end separately;
3. paired end mapping with inner distance between the two ends to have a mean of 0 and sd=50;

Mapping rate for 1 and 3 are both around 45%. But mapping rate for 2 is around 75%. So most of the reads are mappable but cannot be paired even when the inner distance between the two ends can be 0 or negative (if the distance can go negative). I took a look at the BAM files for the accepted hits. It looks like read pairs with alignment

------------>>>----------------
------------<<<<----------------

can be paired. But read pairs with alignment

-------------->>>>--------------
-----------<<<<---------------

can not be paired. Does any one know why and how I can gain those unpaired alignment back? Thanks in advance!
Markovia is offline   Reply With Quote
Reply

Tags
rnaseq low mapping rate

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 12:28 AM.


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