SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA amount for RNA seq, Illumina basager Illumina/Solexa 12 04-09-2014 08:48 AM
Directional RNA-seq: Illumina Tru-seq versus dUTP based method jazz Sample Prep / Library Generation 35 06-06-2013 11:50 AM
RNA-Seq: High-Throughput Illumina Strand-Specific RNA Sequencing Library Preparation. Newsbot! Literature Watch 1 01-22-2013 12:48 AM
RNA-Seq: Single Read and Paired End mRNA-Seq Illumina Libraries from 10 Nanograms Tot Newsbot! Literature Watch 0 11-09-2011 03:10 AM
RNA-Seq: De Novo Transcriptome Sequencing in Anopheles funestus Using Illumina RNA-Se Newsbot! Literature Watch 0 12-15-2010 03:00 AM

Reply
 
Thread Tools
Old 02-09-2010, 07:46 AM   #1
ik76
Junior Member
 
Location: London, UK

Join Date: Jan 2010
Posts: 9
Default problem with Illumina PE RNA seq alingment

Hi All,
We recently started to do PE RNA seq and we encountered a problem the cause of which have eluded us so far. In summary, following cycle 30 we see rapid increase in alignment error. The error rates per cycle go up to 20% and the % reads with 1 or more mismatches goes up to 35-50%.
What we have found is:
1) The error is alingner-nonspecific ( we tested eland, bowtie and even commercial software).


2) The error is not due to contamination with amplified adapters or low quality scores for the bases in the second half of the amplicons.

3) The problem is RNA-specific. Both PhiX and genomic DNA samples had normal error and alignment plots.

4) The problem seems to depend on both the length of the insert and the length of the run. During the same run, shorter insert libraries have more misaligned reads than libraries consisting of longer fragments. A given insert size library will have more misaligned reads during a longer run.

5) There is nothing visibly and strongly wrong with the IVC plots of the affected samples. The only striking observation is that both the % basecall as well as (even more strongly) the % all intensities from the affected samples form some type of ‘fork’ structures in the later cycles of each read.

Has anybody heard about such a problem or encountered it? Any suggestions are welcomed.

We know that the problem is related to the sample prep and taking longer reads ( >200bp insert size) will help alleviate it but we wonder what is the cause of it?
ik76 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 05:33 PM.


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