High scoring segment selection for pairwise whole genome sequence alignment with the maximum scoring subsequence and GPUs Online publication date: Thu, 13-Feb-2020
by Abdulrhman Aljouie; Ling Zhong; Usman Roshan
International Journal of Computational Biology and Drug Design (IJCBDD), Vol. 13, No. 1, 2020
Abstract: Whole genome alignment programs use string matching with hash tables to identify high scoring fragments between a query and target sequence around which a full alignment is then built. A recent study comparing alignment programs showed that while evolutionary similar genomes were easy to align, divergent genomes still posed a challenge to existing methods. To fill this gap we explore the use of the maximum scoring subsequence to identify high scoring fragments. We split the query genome into several fragments and align them to the target with a previously published parallel algorithm for short read alignment. We then pass such high scoring fragments on to the LASTZ program to obtain a more complete alignment. On simulated data we obtain an average of at least 20% higher accuracy than the alignment given by LASTZ at the expense of few hours of additional runtime. Our source code is freely available at http://web.njit.edu/usman/MSGA
Existing subscribers:
Go to Inderscience Online Journals to access the Full Text of this article.
If you are not a subscriber and you just want to read the full contents of this article, buy online access here.Complimentary Subscribers, Editors or Members of the Editorial Board of the International Journal of Computational Biology and Drug Design (IJCBDD):
Login with your Inderscience username and password:
Want to subscribe?
A subscription gives you complete access to all articles in the current issue, as well as to all articles in the previous three years (where applicable). See our Orders page to subscribe.
If you still need assistance, please email subs@inderscience.com
Our Blog 
Follow us on Twitter 
Visit us on Facebook