SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
MapView: a viewer for short reads alignment baohua100 Bioinformatics 99 08-23-2013 08:12 AM
Bowtie: Ultrafast and memory-efficient alignment of short reads to the human genome Ben Langmead Literature Watch 2 03-04-2013 02:06 AM
Slider - Maximum use of probability information for alignment of short sequence reads ECO Bioinformatics 17 09-21-2010 04:35 PM
SOAP: short oligonucleotide alignment program. ECO Bioinformatics 22 10-26-2009 01:55 PM
Haman genome alignment with short reads ptongyoo Bioinformatics 4 04-14-2009 05:27 PM

Reply
 
Thread Tools
Old 03-07-2011, 04:34 AM   #1
rboettcher
Member
 
Location: Berlin

Join Date: Oct 2010
Posts: 71
Default Short sequence alignment

Hi all,
I am currently looking for a software to do alignments of short reads against other short sequences (so not the whole genome). The only problem is: I want to do a lot of these alignments. So instead of aligning Xmio reads to one big sequence, my goal is to align them against 10k (or so) very short sequences. Does anybody know a suitable software?

Regards

€: I should add that I'm not aiming on multiple sequence alignments.

Last edited by rboettcher; 03-07-2011 at 04:58 AM.
rboettcher is offline   Reply With Quote
Old 03-07-2011, 05:01 AM   #2
henry.wood
Member
 
Location: Leeds, UK

Join Date: Apr 2010
Posts: 63
Default

We have done something like this using bwa. We made an artificial reference 'genome' from 3000 short (5-15Kb) fragments treating each fragment as a chromosome. We then aligned several lanes of Illumina data against this quite happily.
henry.wood is offline   Reply With Quote
Old 03-07-2011, 06:26 AM   #3
rboettcher
Member
 
Location: Berlin

Join Date: Oct 2010
Posts: 71
Default

I would prefer something not relying on prior indexing, so a fast implementation of Needleman-Wunsch or Hirschberg-algorithm would be sufficient I suppose, as long as it is accessible by console.
rboettcher is offline   Reply With Quote
Old 03-07-2011, 06:42 AM   #4
zee
NGS specialist
 
Location: Malaysia

Join Date: Apr 2008
Posts: 249
Default

Novoalign does full Needleman-Wunsch after seeding matches. You can also map short reads to small fragments.
zee is offline   Reply With Quote
Old 03-09-2011, 07:45 PM   #5
sparks
Senior Member
 
Location: Kuala Lumpur, Malaysia

Join Date: Mar 2008
Posts: 126
Default

Adding to Zees comment, Novoalign indexing step is very fast for small references so it shouldn't be an issue.

I was also wondering if you're interested in alignments that overlap rather than are contained in your short target sequences. If you are then you could add N's to both ends of the short sequences that you index.
sparks is offline   Reply With Quote
Old 03-11-2011, 06:00 AM   #6
rboettcher
Member
 
Location: Berlin

Join Date: Oct 2010
Posts: 71
Default

No overlaps, just aligning the sequences to a number of short reference sequences. Unfortunately, I couldn't get Novoalign to run on my system
rboettcher is offline   Reply With Quote
Old 03-11-2011, 08:58 PM   #7
sparks
Senior Member
 
Location: Kuala Lumpur, Malaysia

Join Date: Mar 2008
Posts: 126
Default

What problem did you have with Novoalign?
sparks is offline   Reply With Quote
Old 03-12-2011, 08:25 AM   #8
ericgdp
Junior Member
 
Location: Belo Horizonte, Brasil

Join Date: Mar 2010
Posts: 2
Default

May be you can use Abacas.
ericgdp is offline   Reply With Quote
Old 03-14-2011, 07:51 PM   #9
biznatch
Senior Member
 
Location: Canada

Join Date: Nov 2010
Posts: 124
Default

PatternHunter?

http://www.bioinformaticssolutions.com/products/ph/
biznatch is offline   Reply With Quote
Old 03-17-2011, 02:35 AM   #10
rboettcher
Member
 
Location: Berlin

Join Date: Oct 2010
Posts: 71
Default

Thank you all very much for your suggestions, but I guess none of the aligners mentioned is really meant to do what I intend it to do, as far as I can tell from studying the manuals/faqs.
So after some consideration, I decided to implement a simple NW-algorithm and after switching from Perl to Java the speed is acceptable.
I had hoped I could skip this step but it seems there is no real alternative.
Thanks anyway, I guess you will read about my next problems soon enough on the forums

Best Regards
rboettcher 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:50 AM.


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