SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Optimizing settings for mapping 100bp RNA-seq reads with tophat seqhorn Bioinformatics 3 10-21-2011 11:40 PM
RNA-Seq: Duplex-specific nuclease efficiently removes rRNA for prokaryotic RNA-seq. Newsbot! Literature Watch 0 09-02-2011 03:00 AM
Combined mapping of RNA-Seq reads originating from multiple species schelhorn RNA Sequencing 7 11-05-2010 09:55 AM
RNA-Seq: MapSplice: Accurate mapping of RNA-seq reads for splice junction discovery. Newsbot! Literature Watch 2 10-14-2010 09:35 AM
RNA-seq reads mapping to coding regions m!x RNA Sequencing 0 02-17-2010 01:04 PM

Reply
 
Thread Tools
Old 04-30-2011, 07:00 PM   #1
nivea
Member
 
Location: NC

Join Date: Apr 2011
Posts: 17
Question Can BLAST do the RNA-seq reads mapping efficiently?

Hi all,

I'm dealing with low-quality reads of a bacteria which are 100 base long. Thus, most of the software cannot map back most of the reads to the reference genome with no more than three mismatches. So can I use BLAST to do the mapping work? Will it to slow?

Thanks
nivea is offline   Reply With Quote
Old 05-02-2011, 01:19 AM   #2
Thorondor
Member
 
Location: Heidelberg

Join Date: Feb 2011
Posts: 69
Default

amount of reads? size of your reference genome? what software you tried already?
Thorondor is offline   Reply With Quote
Old 05-02-2011, 01:56 AM   #3
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

if you think you have more than three mismatches, why do you limit your aligner to it? (Your questions sounds a bit as if you think that most short-read aligners had a hard-coded limit to the maximum number of mismatches allowed. They don't.)
Simon Anders is offline   Reply With Quote
Old 05-02-2011, 02:08 AM   #4
Thorondor
Member
 
Location: Heidelberg

Join Date: Feb 2011
Posts: 69
Default

if you have illulmina reads you might want to read this review:
http://seqanswers.com/forums/showthread.php?t=11045
Thorondor is offline   Reply With Quote
Old 05-02-2011, 04:06 AM   #5
doc.ramses
Member
 
Location: Planet Earth

Join Date: Jan 2011
Posts: 26
Default

Quote:
Originally Posted by Thorondor View Post
if you have illulmina reads you might want to read this review:
http://seqanswers.com/forums/showthread.php?t=11045
Can someone send me the paper via PM ?
Thank you !
doc.ramses is offline   Reply With Quote
Old 05-02-2011, 07:55 AM   #6
nivea
Member
 
Location: NC

Join Date: Apr 2011
Posts: 17
Default

Quote:
Originally Posted by Thorondor View Post
amount of reads? size of your reference genome? what software you tried already?
Hi Thorondor,

I used bowtie to map32129789 reads (76bp) to the bacteria E.coli K12 with a tolerance of mismatches <=3. However, only 0.02% of the reads mapped uniquely to the genome.
After using BLAST check the result I found for most of the reads, only the previous 35b of reads got hit in the genome. Thus, it's probably because of the sequencing error that I cannot use bowtie to map them back.

So I want to try BLAST. But do you think it's too slow? Or does this program work accurately enough for reads mapping?

Thanks for reply~
nivea is offline   Reply With Quote
Old 05-02-2011, 08:07 AM   #7
nivea
Member
 
Location: NC

Join Date: Apr 2011
Posts: 17
Default

Quote:
Originally Posted by Simon Anders View Post
if you think you have more than three mismatches, why do you limit your aligner to it? (Your questions sounds a bit as if you think that most short-read aligners had a hard-coded limit to the maximum number of mismatches allowed. They don't.)
Hi Simon,


The problem is the quality of the reads, which you can see from the upper post. The sequencing error made only smaller than 35bp of most of the reads can map back to the reference genome. If I want to take advantage of the current reads, do you think it's good enough for me to use BLAST?

Thanks~
nivea is offline   Reply With Quote
Old 05-02-2011, 09:09 AM   #8
Simon Anders
Senior Member
 
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 994
Default

Actually, Bowtie should disregard mismatches on low-quality bases, but this does not work that well.

However, why don't you simply trim off the bad-quality parts and give the trimmed reads to Bowtie? In a bacterial genome, 35 bp should be more than enough to get good matches.

regarding BLAST: It is way too slow for millions of reads, and its model of scoring mismatches is based on assuming mutations rather than read errors to be the cause, which is inappropriate.
Simon Anders is offline   Reply With Quote
Old 05-02-2011, 09:15 AM   #9
Thorondor
Member
 
Location: Heidelberg

Join Date: Feb 2011
Posts: 69
Default

so why not trim the reads first and then try to map without quality values? You don't want to assemble bad reads anyway and they mostly won't overlap if you don't trim. Or you could just cut all your fastq reads to a length of 40 bp and then try to map with BWA.

Maybe you also want to have a look at Cufflinks for your assembly? http://cufflinks.cbcb.umd.edu/

edit: simon was faster :-/
Thorondor is offline   Reply With Quote
Old 05-02-2011, 09:16 AM   #10
nivea
Member
 
Location: NC

Join Date: Apr 2011
Posts: 17
Default

Quote:
Originally Posted by Simon Anders View Post
Actually, Bowtie should disregard mismatches on low-quality bases, but this does not work that well.

However, why don't you simply trim off the bad-quality parts and give the trimmed reads to Bowtie? In a bacterial genome, 35 bp should be more than enough to get good matches.

regarding BLAST: It is way too slow for millions of reads, and its model of scoring mismatches is based on assuming mutations rather than read errors to be the cause, which is inappropriate.
Thanks Simon! That's very helpful.
nivea is offline   Reply With Quote
Old 05-02-2011, 09:33 AM   #11
nivea
Member
 
Location: NC

Join Date: Apr 2011
Posts: 17
Default

Quote:
Originally Posted by Thorondor View Post
so why not trim the reads first and then try to map without quality values? You don't want to assemble bad reads anyway and they mostly won't overlap if you don't trim. Or you could just cut all your fastq reads to a length of 40 bp and then try to map with BWA.

Maybe you also want to have a look at Cufflinks for your assembly? http://cufflinks.cbcb.umd.edu/

edit: simon was faster :-/

Thanks Thorondor!
nivea 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 07:17 AM.


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