What Is The Fastest Alignment Tool To Map Reads To A Big Genome Over 4Gb?
5
3
Entering edit mode
13.1 years ago
Dejian ★ 1.3k

We assembled a genome around 5GB. Now we are trying to map the Illumina reads to the genome. There are many aligners that can do this job. But which is the fastest one? Have anyone benchmarked any of these tools?

The list of aligners is as follows: short-read aligners.

short aligner • 6.2k views
ADD COMMENT
6
Entering edit mode
13.1 years ago

Asking for the fastest tool is posing the question the wrong way.

What you need is finding a fast enough method for your data. As Jeremy points out there are many questions that need to be answered and you will need to run your own tests to identify which tool is fast enough and provides the necessary coverage as well. Your genome is likely to contain vast streches of low complexity regions - that could make some approaches perform differently.

Start with the known fast tools like bowtie and SOAP and see what these do. Also read the literature: Comparative analysis of algorithms for next-generation sequencing read alignment

ADD COMMENT
3
Entering edit mode
13.1 years ago
Ido Tamir 5.2k

There are some comparisons the latest to my knowledge is Ruffalo. But is being the fastest really an important criteria, if the difference to more accurate results is 2x running time?

ADD COMMENT
0
Entering edit mode

Currently we want to get a rough result of the coverage, so speed is the first consideration. Thanks for the paper recommended.

ADD REPLY
3
Entering edit mode
13.1 years ago

you need a lot more criteria to make this a fair fight
here are some just off the top of my head

  1. number of mismatches
  2. memory ceiling
  3. paired-end?
  4. quality-aware?
  5. sensitivity
ADD COMMENT
0
Entering edit mode

Currently we want to get a rough result of the coverage, so speed is the first consideration. Thanks for your suggestion.

ADD REPLY
0
Entering edit mode

your comment indicated that you do not understand that accuracy and sensitivity are not independent from coverage. Increasing these will necessarily slow down the method.

ADD REPLY
1
Entering edit mode
6.7 years ago
predeus ★ 2.0k

I think it would be interesting to bump it up in 2018.

What is the fastest aligners out there? I have heard (I think Heng Li wrote it somewhere) that with changes in computer architecture, hash-based aligners are coming to be more and more interesting, compared to FM-index based. I am about to try our SNAP, which promises up to 20x increase compared to bowtie2 and bwa.

Like Istvan rightly pointed out, aligner has to be fast for your application. But sometimes speed could be a limiting factor - e.g. when you are teaching a class and want to show something pretty quickly, and not wait for 30 minutes.

ADD COMMENT
0
Entering edit mode
13.1 years ago
Gjain 5.8k

I use Novoalign (http://www.novocraft.com/main/page.php?s=novoalign) which is good and fast.

ADD COMMENT
1
Entering edit mode

is novoalign really fast? I've found it to do a good job, but not especially fast.

ADD REPLY
0
Entering edit mode

it depends on your definition of what you call fast. If I run my mapping (HG19) pipeline on our cluster with around 800 nodes ... it normally takes 4-5 hours.

ADD REPLY
0
Entering edit mode

@Gjain: Do you use MPI version? or ...

ADD REPLY
0
Entering edit mode

just normal version

ADD REPLY
0
Entering edit mode

Hi,I‘m using a software named LSC,which need novoalign to align.But I can’t download novoalign.And I can't open tihs URL.Can you send it to me? My email address is liuyanhu005@gmail.com. Many thanks!

ADD REPLY
3
Entering edit mode

Hi, My apologies for reviving an old thead but I've receently become aware that the default options in LSC for Novoalign are far from ideal.

If you are using Novoalign with LSC you need to use -r Exhaustive not -r All. This will improve results dramatically and should be better than the other aligners.

For 100bp reads try.. novoalign_options = -c1 -r Ex 1000 -t 120 -F FA -g 0 -x 20 -o sam

If read length is different adjust -t accordingly. On Cerebellum reads used by Kin Fai results are 40% better with 30% reduced run time vs BWA using LSC 0.3.1

Best, Colin

ADD REPLY
0
Entering edit mode

This is very useful ... thanks Colin

ADD REPLY
0
Entering edit mode

hi, please use this page to review the licence agreement and you can download the one which suits your needs best. http://www.novocraft.com/main/page.php?id=968

ADD REPLY
0
Entering edit mode

Thank you a lot!

ADD REPLY

Login before adding your answer.

Traffic: 1543 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6