SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
help calculating unique mapping rate for paired end emolinari Bioinformatics 2 03-22-2015 02:51 AM
Multiple Pair-end Reads phynex92 Illumina/Solexa 10 07-08-2014 12:22 PM
Low mapping rate for RNAseq 2x150 trimmed data Markovia Bioinformatics 0 08-30-2013 07:55 PM
TopHat pair-end RNAseq statistics Xiaomin Yu Bioinformatics 0 02-20-2012 09:02 AM
Low pairing rate in SOLiD 4 pair-end transcriptome sequencing amurocw Bioinformatics 7 04-13-2011 08:32 AM

Reply
 
Thread Tools
Old 06-01-2015, 08:42 AM   #1
TTYE
Junior Member
 
Location: Maryland

Join Date: Jun 2015
Posts: 7
Default High multiple mapping rate Drosophila pair-end 50 bp RNAseq

Hi, I used Tophat2 for alignment. The mapping rate is ~92% while the multiple mapping rate is ~25%. Based on the FastQC result, it has some problem in the per tile quality, I see many red strips. So I trimmed the low quality head and tail, but the mapping rates and multiple mapping rates are not changed much. What would be the reason? How can I fix it? Thanks!!!
TTYE is offline   Reply With Quote
Old 06-01-2015, 08:49 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,989
Default

Post the FastQC plots (quality, adapter) to give us an idea of the "problem" you are referring to. Did you do trimming based on quality (what was the cutoff used) or just blanket trimming for a certain number of bases at beginning/end of read?
GenoMax is offline   Reply With Quote
Old 06-01-2015, 01:50 PM   #3
TTYE
Junior Member
 
Location: Maryland

Join Date: Jun 2015
Posts: 7
Default per tile picture attached

Hi GenoMax, Thanks for your quick reply.

I should be more accurate, I trimmed the sequences based on the quality score. Do you think filter the whole sequence with low quality segment will be better?

Just looking at the per tile figure, what could be the reason of such situation. I checked the FastQC webpage, the author suggest that maybe the bubbles in the flow cell. What do you think?
Attached Images
File Type: png per_tile_quality.png (11.4 KB, 13 views)
TTYE is offline   Reply With Quote
Old 06-01-2015, 02:11 PM   #4
TTYE
Junior Member
 
Location: Maryland

Join Date: Jun 2015
Posts: 7
Default trimming parameters

I used trimmomatic with parameters:
LEADING:20 \
TRAILING:20 \
SLIDINGWINDOW:4:20

Thanks!

Quote:
Originally Posted by GenoMax View Post
Post the FastQC plots (quality, adapter) to give us an idea of the "problem" you are referring to. Did you do trimming based on quality (what was the cutoff used) or just blanket trimming for a certain number of bases at beginning/end of read?
TTYE is offline   Reply With Quote
Old 06-01-2015, 03:55 PM   #5
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,989
Default

Have you asked the sequencing facility if there was any problem with the run this sample was on? What kind of sequencer (MiSeq/HiSeq) is this data from? Particular tiles have been affected across entire run so unless there were multiple bubbles stuck at certain positions in the flowcell (seems unlikely) ...
GenoMax is offline   Reply With Quote
Old 06-02-2015, 08:17 AM   #6
TTYE
Junior Member
 
Location: Maryland

Join Date: Jun 2015
Posts: 7
Default

Just confirmed by the sequencing core, their machine does go wrong around that time. So what should I do with it? Thanks!
TTYE is offline   Reply With Quote
Old 06-02-2015, 08:41 AM   #7
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,989
Default

They should not have released the data in the first place, if there was a known problem with the run. I suppose you can ask them to re-run your sample (if this was a machine/reagent problem then Illumina will generally provide free replacements, if your facility has a maintenance agreement).

Going back to the original question of multi-mapping .. that may be a characteristic of your sample. I would imagine that run related problems should not change the composition of your sequences (this can be validated if you do get a new run done from your facility).
GenoMax is offline   Reply With Quote
Old 06-03-2015, 06:22 AM   #8
TTYE
Junior Member
 
Location: Maryland

Join Date: Jun 2015
Posts: 7
Default

Got you. Thanks a lot!
TTYE is offline   Reply With Quote
Reply

Tags
multiple alignment, pair end, per tile quality

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:34 AM.


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