SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samstat versus tophat overall alignment rate % alyamahmoud RNA Sequencing 1 12-26-2013 07:11 AM
Tophat2 normal alignment rate JJP RNA Sequencing 7 08-12-2013 09:16 AM
Determining alignment rate for Tophat all_your_base Bioinformatics 0 12-12-2012 07:04 AM
Tophat2 Alignment Rate SHeaph Bioinformatics 5 11-21-2012 01:05 AM
Alignment Cover for Simulated Illumina Reads BertieWooster Bioinformatics 1 09-13-2010 07:15 PM

Reply
 
Thread Tools
Old 01-15-2014, 08:23 PM   #1
alyamahmoud
Member
 
Location: Egypt, Saudi Arabia

Join Date: Nov 2013
Posts: 29
Default simulated 50 bp SE gives higher alignment rate compared to 100 bp PE

Hi All

I have 100 bp PE ss RNASeq reads. I tried to simulate 50 bp SE out of those by trimming using trimmomatic and then mapping both 100 PE and simulated 50 bp SE using tophat. I used -g 1 to make sure I get reads counted only once.

Comparing the alignment rate; the 50 bp SE have higher number of mapped reads in the accepted_hits.bam output of tophat. I also checked the correlation between log2 fold change of DEG and they seem to be very well correlated between the two experiments.

Can anyone explain to me why would 50 bp SE reads map better than 100 bp PE ?

Thanks a lot
Alyaa
alyamahmoud is offline   Reply With Quote
Old 01-16-2014, 03:00 PM   #2
Dario1984
Senior Member
 
Location: Sydney, Australia

Join Date: Jun 2011
Posts: 166
Default

It's because for single-end reads you only have 50 bases to map, whereas for paired-end you have to map two sets of 50 bases correctly. Assuming sequencing errors happen at the same rate, there would be twice as many errors to overcome for mapping the paired-end data.
Dario1984 is offline   Reply With Quote
Reply

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 09:24 PM.


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