AL2CO - calculation of position conservation for a multiple sequence alignment
Options: sequence weighting scheme: Sequence weighting corrects for biased representation of similar sequences in an alignment. Options: 0=unweighted, 1=weighted by the modified method of Henikoff & Henikoff, 2=independent-count based Default = 2 conservation calculation method: Options: 0=entropy-based C(i)=sum_{a=1}^{20}f_a(i)*ln[f_a(i)], where f_a(i) is the frequency of amino acid a at position i, 1=variance-based C(i)=sqrt[sum_{a=1}^{20}(f_a(i)-f_a)^2], where f_a is the overall frequency of amino acid a, 2=sum-of-pairs measure C(i)=sum_{a=1}^{20}sum_{b=1}^{20}f_a(i)*f_b(i)*S_{ab}, where S_{ab} is the element of a scoring matrix for amino acids a and b Default = 0 scoring matrix This is only for the conservation calculation method of sum-of-pairs measure. Options: 0=BLOSUM62 matrix 1=identity matrix Default = 0 scoring matrix transformation This is only for the conservation calculation method of sum-of-pairs measure. Options: 0=no transformation, 1=normalization S'(a,b)=S(a,b)/sqrt[S(a,a)*S(b,b)], 2=adjustment S"(a,b)=2*S(a,b)-(S(a,a)+S(b,b))/2 Default = 0 window size used for averaging conservation (for smoothing purpose): Default = 1 Recommended value for motif analysis: 3 normalize conservation values: Subtract the mean from each conservation index and divide by the standard deviation. Default = True exclude the first sequence from calculation: If set to true, the first sequence in the alignment will not be included in the conservation calculation. Default = False gap fraction above which conservation calculation is not performed: The value should be more than 0 and no more than 1. Conservation indices are calculated only for positions with gap fraction less than the specified value. Otherwise, conservation indices will be set to M-S, where M is the mean conservation value and S is the standard deviation. Default = 0.5 Block size of the output alignment file with conservation indices: The number of residues in a line of each block in output alignment with conservation indices. Default = 70 pdb file for which b-factor field is replaced with conservation (optional): The sequence in the pdb file should match exactly the first sequence of the alignment. Results include a link to a pdb file with b-factor field replaced by conservation values so that conservation information can be mapped onto a structure and displayed.