SOAP3-dp: Fast, Accurate and Sensitive GPU-based Short Read Aligner

Authors:
  • RLRuibang Luo
    Ruibang Luo
  • TWThomas Wong
    Thomas Wong
  • JZJianqiao Zhu
    Jianqiao Zhu
  • CLChi-Man Liu
    Chi-Man Liu
  • EWEdward Wu
    Edward Wu
  • LLLap-Kei Lee
    Lap-Kei Lee
  • HLHaoxiang Lin
    Haoxiang Lin
  • WZWenjuan Zhu
    Wenjuan Zhu
  • DCDavid W. Cheung
    David W. Cheung
  • HTHing-Fung Ting
    Hing-Fung Ting
  • SYSiu-Ming Yiu
    Siu-Ming Yiu
  • CYChang Yu
    Chang Yu
Preprint
Publisher website
Google Scholar
Abstract
To tackle the exponentially increasing throughput of Next-Generation Sequencing (NGS), most of the existing short-read aligners can be configured to favor speed in trade of accuracy and sensitivity. SOAP3-dp, through leveraging the computational power of both CPU and GPU with optimized algorithms, delivers high speed and sensitivity simultaneously. Compared with widely adopted aligners including BWA, Bowtie2, SeqAlto, GEM and GPU-based aligners including BarraCUDA and CUSHAW, SOAP3-dp is two to tens of times faster, while maintaining the highest sensitivity and lowest false discovery rate (FDR) on Illumina reads with different lengths. Transcending its predecessor SOAP3, which does not allow gapped alignment, SOAP3-dp by default tolerates alignment similarity as low as 60 percent. Real data evaluation using human genome demonstrates SOAP3-dp's power to enable more authentic variants and longer Indels to be discovered. Fosmid sequencing shows a 9.1 percent FDR on newly discovered deletions. SOAP3-dp natively supports BAM file format and provides a scoring scheme same as BWA, which enables it to be integrated into existing analysis pipelines. SOAP3-dp has been deployed on Amazon-EC2, NIH-Biowulf and Tianhe-1A.