MEGABLAST

# MEGABLAST Search

Mega BLAST uses the greedy algorithm for nucleotide sequence alignment search. This program is optimized for aligning sequences that differ slightly as a result of sequencing or other similar "errors". When larger word size is used (see explanation below), it is up to 10 times faster than more common sequence similarity programs. Mega BLAST is also able to efficiently handle much longer DNA sequences than the blastn program of traditional BLAST algorithm.

## Parameters

### Word size

Word size is roughly the minimal length of an identical match an alignment must contain if it is to be found by the algorithm. Mega BLAST is most efficient with word sizes 16 and larger, although word size as low as 8 can be used.
If the value W of the word size is divisible by 4, it guarantees that all perfect matches of length W + 3 will be found and extended by Mega BLAST search, however perfect matches of length as low as W might also be found, although the latter is not guaranteed. Any value of W not divisible by 4 is equivalent to the nearest value divisible by 4 (with 4i+2 equivalent to 4i).

### Percent identity

If this parameter P is set, only the alignments with identity percentage higher than P will be retained. Also the default match reward and mismatch penalty scores are chosen in this case close to the log-odds (i.e. the most statistically effective) scores for the PAM distance corresponding to a sequence conservation level somewhat higher than P. The following table shows the relation between the percent identity cut-off values, the target conservation levels and the corresponding log-odds match and mismatch scores used by Mega BLAST:

Percent Identity | Target | Match score | Mismatch score
None | 95 | 1 | -2
>= 95 | 99 | 1 | -3
85, 90 | 95 | 1 | -2
80 | 88 | 2 | -3
75 | 83 | 4 | -5

### Gapping parameters

By default, non-affine gapping parameters are assumed. This means that the gap opening penalty is 0, and gap extension penalty E can be computed from match reward r and mismatch penalty q by the formula: E = r/2 - q. The non-affine version of Mega BLAST requires significantly less memory and is also significantly faster, however affine gapping parameters can also be used, preferably with larger word sizes. Non-affine gapping parameters tend to yield alignments with more gaps, but the gap lengths are shorter.

### X-dropoff value

As in BLAST, this value provides a cutoff threshold for the extension algorithm tree exploration. When the score of a given branch drops below the current best score minus the X-dropoff, the exploration of this branch stops. However the actual values of the X-dropoff for Mega BLAST and for traditional nucleotide BLAST algorithms are not necessarily compatible, i.e. with the same word size, match, mismatch and gapping penalties and with the same X-dropoff, the two algorithms might produce different results, which can be remedied by changing the X-dropoff value for one of the algorithms.