Alignment method

Blastp

SmithWaterman


Substitution matrix

BLOSUM62

BLOSUM62

DirAtPf100

Statistics
   
Pvalue (KarlinAltschul)

0.008

NA

NA

Zvalue (PearsonLipman)

10

11

12

Tvalue (TULIP theorem)

0.01

8.10^{3}

7.10^{3}

Pvalue (this work)

1.5.10^{6}

3.7.10^{7}

1.10^{7}

 TFIIA gamma sequences from Plasmodium (UniProtKB Q8I4S7_PLAF7) and Arabidopsis (UniProtKB T2AG_ARATH) were aligned with Blastp and SmithWaterman methods. Statistics were computed following the KarlinAltschul model (as implemented in the Blastp algorithm) or the LipmanPearson Zvalue model. The upper bound for the Pvalue based on the TULIP theorem is given following the formula: Tvalue = 1/Zvalue^{2}. The Pvalue deduced from the Zvalue Gumbel distribution was computed following the model presented here. Substitution matrices were either BLOSUM62, or the asymmetric DirAtPf100 matrix specified for Plasmodium vs. Arabidopsis comparisons. NA: not applicable.