blastn Parameters sets

Parameter sets for blastn, with notes. Please refer to Posting 4 ("blastn Parameters for noncoding queries")

A -r 10 -q -11 -W 8 -G 100 -E 50

B -r 10 -q -11 -W 7 -G 5 -E 20

C -r 10 -q -11 -W 7 -G 15 -E 4

D -r 7 -q -14 -W 7 -G 7 -E 3

E -r 4 -q -5 -W 8 -G 10 -E 2

F -r 5 -q -4 -W 7 -G 10 -E 4

G -r 8 -q -5 -W 7 -G 20 -E 5

H -r 7 -q -8 -W 7 -G 15 -E 2

I -r 8 -q -7 -W 7 -G 16 -E 2

N -r 1 -q -3 -W 11 -G 5 -E 2

Notes (parameters are listed above, comments are provided below):

A Use these if you are sure that there will be no gaps in your alignment

B Use these to permit only very short gaps.

C Standard snRNA parameters (similar to G and H)

D Very low gap penalties, but high mismatch penalty

E

F my recommendation (start here) . When the Arabidopsis genome is searched with human U12 snRNA using these parameters, the correct top hit is found, and the ratio between the E value of the true hit and of the top false hit is approximately 300).

G Chau's initial set. These are promiscuous (extreme r/q). For example, when the Arabidopsis genome is searched with human U12 snRNA using these parameters, the correct top hit is found, but many random hits score nearly as well (The ratio between the E value of the true hit and of the top false hit is only 10).

H Chau's "best" set. These are also promiscuous (low G and especially low E). When the Arabidopsis genome is searched with human U12 snRNA using these parameters, the correct top hit is again found, but many random hits score nearly as well (The ratio between the E value of the true hit and of the top false hit is only 10).

I Very similar to H. I have recommended these defaults to a lot of people over the years.

N NCBI defaults. You will never find noncoding RNA matches using these parameters.

Good luck! If you have experience that bears on this, or can cite relevant literature, please let me know and I'll update this posting.