Could someone suggest a fast and efficient implementation of a glocal (semi-global) sequence aligner?
3
0
Entering edit mode
5.5 years ago
Tim ▴ 130

I am looking for a program that would be able to perform semi-global (global in query, local in reference sequence) protein alignments fast and efficiently. My query sequences are short (20-30 aa), while reference sequences 10-20 times longer. So far I have been able to find GLSEARCH (from FASTA package) and align, but they seem to be somewhat slow for my task. I would really appreciate it if someone could suggest another alternative. enter image description here

alignment glocal semi-global protein • 3.0k views
ADD COMMENT
1
Entering edit mode
5.5 years ago
cmdcolin ★ 2.5k

At least one software that strongly uses the "glocal" terminology is LAGAN, specifically Shuffle-LAGAN

e.g. this website http://genome.lbl.gov/cgi-bin/VistaInput?align_pgm=slagan&num_seqs=2 uses the phrase

Shuffle-LAGAN Glocal pair-wise alignment of finished sequences (detects rearrangements)

ADD COMMENT
0
Entering edit mode

Thanks, I have checked Shuffle-LAGAN, seems to be a good program and I will test it later. My only concern is that was designed to do glocal alignments of long DNA sequences, while I have short-medium protein sequences (20-500 aa long).

ADD REPLY
1
Entering edit mode
5.5 years ago
l.souza ▴ 80

Parasail (https://github.com/jeffdaily/parasail)

ADD COMMENT
0
Entering edit mode

Thanks, I will check it.

ADD REPLY
0
Entering edit mode
5.5 years ago
anin.gregory ▴ 110

Nucmer, part of the Mummer package (http://mummer.sourceforge.net/).

ADD COMMENT
0
Entering edit mode

Please correct me if I am wrong, but I had an impression that nucmer performs global alignments of long sequences, while I need semiglobal (glocal) alignments of short to moderately long sequences. I updated my original post to reflect the difference between different methods of alignment.

ADD REPLY

Login before adding your answer.

Traffic: 2171 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