![]() |
contactalign |
Standard (Mandatory) qualifiers: [-querysequence] sequence First (query) polypeptide sequence to be aligned [-templatesequence] sequence Second (template) polypeptide sequence to be aligned [-contactscoringfile] matrixf Matrix of scores for frequency of contacts between specific amino acid residue types [-cmaptemplate] infile contact map of template structure -gapopen float The gap opening penalty is subtracted from the total score every time a gap is inserted in the alignment. The default value is chosen for the EBLOSUM62 scoring matrix when you are aligning protein sequences, and the EDNAFULL matrix when you are aligning nucleotide sequences. If you sequence contains many small errors you should set a low gap opening penalty (near zero) to avoid over-penalizing those errors. -gapextend float The gap extension penalty is the score substracted when more spaces are added to an existing alignment gap. Often a single insertion or deletion affects multiple residues/nucleotides; separate opening and extension penalties reduce the excessive penalty such events would otherwise carry; generally the gap extension penalty should be significantly lower than the gap opening penalty. [-verticaltrace] outfile vertical trace of initial alignment [-alignedoutfile] seqout Output sequence USA Additional (Optional) qualifiers: -substitutionscoringfile matrixf Matrix of scores for specific amino acid residue type substitution within one-dimensional sequence strings Advanced (Unprompted) qualifiers: -[no]zeroend boolean If this option is true then no penalty is added for gaps at either ends of a trace. Associated qualifiers: "-querysequence" associated qualifiers -sbegin1 integer Start of the sequence to be used -send1 integer End of the sequence to be used -sreverse1 boolean Reverse (if DNA) -sask1 boolean Ask for begin/end/reverse -snucleotide1 boolean Sequence is nucleotide -sprotein1 boolean Sequence is protein -slower1 boolean Make lower case -supper1 boolean Make upper case -sformat1 string Input sequence format -sdbname1 string Database name -sid1 string Entryname -ufo1 string UFO features -fformat1 string Features format -fopenfile1 string Features file name "-templatesequence" associated qualifiers -sbegin2 integer Start of the sequence to be used -send2 integer End of the sequence to be used -sreverse2 boolean Reverse (if DNA) -sask2 boolean Ask for begin/end/reverse -snucleotide2 boolean Sequence is nucleotide -sprotein2 boolean Sequence is protein -slower2 boolean Make lower case -supper2 boolean Make upper case -sformat2 string Input sequence format -sdbname2 string Database name -sid2 string Entryname -ufo2 string UFO features -fformat2 string Features format -fopenfile2 string Features file name "-verticaltrace" associated qualifiers -odirectory5 string Output directory "-alignedoutfile" associated qualifiers -osformat6 string Output seq format -osextension6 string File name extension -osname6 string Base file name -osdirectory6 string Output directory -osdbname6 string Database name to add -ossingle6 boolean Separate file for each entry -oufo6 string UFO features -offormat6 string Features format -ofname6 string Features file name -ofdirectory6 string Output directory General qualifiers: -auto boolean Turn off prompts -stdout boolean Write standard output -filter boolean Read standard input, write standard output -options boolean Prompt for standard and additional values -debug boolean Write debug output to program.dbg -verbose boolean Report some/full command line options -help boolean Report command line options. More information on associated and general qualifiers can be found with -help -verbose -warning boolean Report warnings -error boolean Report errors -fatal boolean Report fatal errors -die boolean Report deaths |
Standard (Mandatory) qualifiers | Allowed values | Default | |
---|---|---|---|
[-querysequence] (Parameter 1) |
First (query) polypeptide sequence to be aligned | Readable sequence | 1hmp_a.fasta |
[-templatesequence] (Parameter 2) |
Second (template) polypeptide sequence to be aligned | Readable sequence | 1hgx_a.fasta |
[-contactscoringfile] (Parameter 3) |
Matrix of scores for frequency of contacts between specific amino acid residue types | Comparison matrix file in EMBOSS data path | EBLOSUM62 for protein EDNAFULL for DNA |
[-cmaptemplate] (Parameter 4) |
contact map of template structure | Input file | 1hgx_a.con |
-gapopen | The gap opening penalty is subtracted from the total score every time a gap is inserted in the alignment. The default value is chosen for the EBLOSUM62 scoring matrix when you are aligning protein sequences, and the EDNAFULL matrix when you are aligning nucleotide sequences. If you sequence contains many small errors you should set a low gap opening penalty (near zero) to avoid over-penalizing those errors. | Floating point number from 1.0 to 100.0 | 10.0 for any sequence |
-gapextend | The gap extension penalty is the score substracted when more spaces are added to an existing alignment gap. Often a single insertion or deletion affects multiple residues/nucleotides; separate opening and extension penalties reduce the excessive penalty such events would otherwise carry; generally the gap extension penalty should be significantly lower than the gap opening penalty. | Floating point number from 0.0 to 10.0 | 0.5 for any sequence |
[-verticaltrace] (Parameter 5) |
vertical trace of initial alignment | Output file | vertical.trace |
[-alignedoutfile] (Parameter 6) |
Output sequence USA | Writeable sequence | <sequence>.format |
Additional (Optional) qualifiers | Allowed values | Default | |
-substitutionscoringfile | Matrix of scores for specific amino acid residue type substitution within one-dimensional sequence strings | Comparison matrix file in EMBOSS data path | ../../../emboss/data/EBLOSUM62 |
Advanced (Unprompted) qualifiers | Allowed values | Default | |
-[no]zeroend | If this option is true then no penalty is added for gaps at either ends of a trace. | Boolean value Yes/No | Yes |
Under development, no example input file(s) yet
Under development, no example output file(s) yet
Program name | Description |
---|---|
alignrunner | Align sequence pairs in a directory |
comparator | Compare contact scores of two sequence alignments |
esim4 | Align an mRNA to a genomic DNA sequence |
est2genome | Align EST and genomic DNA sequences |
nawalign | Damian Counsell's NW implementation |
nawalignrunner | Nawalign all sequence pairs in a directory |
needle | Needleman-Wunsch global alignment |
needlerunner | Needle all sequence pairs in a directory |
scorer | Score alignments using structural alignments |
scorerrunner | SCORER for ordered pairs of substituted seqs |
stretcher | Finds the best global alignment between two sequences |
substitute | Substitute matches into a template |
substituterunner | Run SUBSTITUTE on a directory of traces |