SEQanswers

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



Similar Threads
Thread Thread Starter Forum Replies Last Post
bowtie index problem (bowtie-build and then bowtie-inspect) tgenahmet Bioinformatics 4 09-10-2013 11:51 AM
454 GS FLX+ long reads question eladSeq 454 Pyrosequencing 19 01-17-2012 10:16 AM
bowtie 80 mer-long mapping kenosaki Bioinformatics 5 07-03-2011 02:24 AM
problem aligning SOLiD reads with Bowtie pgalante SOLiD 2 07-09-2010 07:29 AM
Velvet and long Illumina reads Peter Bjarke Olsen De novo discovery 15 05-13-2010 10:13 PM

Reply
 
Thread Tools
Old 06-22-2010, 07:10 AM   #1
cp229
Junior Member
 
Location: Leicester

Join Date: Jun 2010
Posts: 5
Default bowtie with long reads, is it a problem?...help!

Hi everyone, I am a very new beginner with bowtie. And I have some trouble about reads...When I run bowtie with long reads(the read length is nearly 200), the misaligned percentage is nearly 98% (I tried human, yeast and e_coli...)
The result is very weird and I cannot find the reason for it...Is there anyone who knows how to solve this problem? Please help me


Thanks!
cp229 is offline   Reply With Quote
Old 06-22-2010, 11:31 PM   #2
simonandrews
Simon Andrews
 
Location: Babraham Inst, Cambridge, UK

Join Date: May 2009
Posts: 870
Default

Bowtie is, by design, a short read aligner. You can make it work with longer reads, but the default setup will not work well there. The major limitation is that it doesn't allow the insertion of gaps in alignments, so that any indel errors will probably make your alignment fail.

Bowtie also has a cutoff for the sum of the mismatched quality scores. The way it's set up you only need around 3 mismatches for this to reject an alignment. If you are using long reads then you will probably want to increase this value to allow for more flexibility in your alignments.
simonandrews is offline   Reply With Quote
Old 06-23-2010, 03:42 AM   #3
cp229
Junior Member
 
Location: Leicester

Join Date: Jun 2010
Posts: 5
Default

Thank you very much...as you say, I tried another e_coli sample which is 35 reads, it`s working nicely. I should use another algorithm to solve the long read problem.
cp229 is offline   Reply With Quote
Old 06-23-2010, 05:13 PM   #4
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

I think for reads longer than 100bp, gapped alignment become important. for long reads (>150bp), try ssaha2/mosaik/bwasw.
lh3 is offline   Reply With Quote
Old 06-24-2010, 02:13 AM   #5
cp229
Junior Member
 
Location: Leicester

Join Date: Jun 2010
Posts: 5
Default

Thx!
When we use mosaik, we found it failed to align long genome, such as yeast and human, it seems that mosaik only works well with very short genome...Is that true??
cp229 is offline   Reply With Quote
Old 07-11-2012, 04:11 AM   #6
Mark
Member
 
Location: Raleigh, NC

Join Date: Nov 2008
Posts: 48
Default

Is bowtie2 more effective given that it permits gaps?
Mark 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 12:38 PM.


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