Fast Analysis of DNA and Protein Sequence on Apple IIe: Restriction Sites Search, Alignment of Short Sequence and Dot Matrix Analysis
Overview
Authors
Affiliations
A fast restriction sites search algorithm using a quadruplet look-ahead feature has been written in 6502 assembly language code. The search time, tested on the sequence of pBR322, is 4.1 s/kilobase using a restriction site library including 112 specificities corresponding to a total site length of over 700 bases. The search for a short sequence (less than 36 bases) within a longer one (up to 9999 bases) with a given number of mismatches or gaps allowed has also been written in assembly language. Typical run time for the search of a 12 base sequence with 1, 2 or 3 gaps allowed are 6.2, 9.4 or 13.6 s/kilobase, respectively. The dot matrix analysis needs 7.5 minutes per square kilobase when using a stringency of 15 matched bases out of 25. A 7/21 matrix of two 500 amino acid proteins is obtained in 3 minutes. These three routines are included in DPSA, a general package of programs allowing manipulation and analysis of DNA and protein sequences.
Hagege J, Brasch M, Cohen S J Bacteriol. 1999; 181(19):5976-83.
PMID: 10498709 PMC: 103624. DOI: 10.1128/JB.181.19.5976-5983.1999.
Hagege J, Pernodet J, Sezonov G, Gerbaud C, Friedmann A, Guerineau M J Bacteriol. 1993; 175(17):5529-38.
PMID: 8366038 PMC: 206609. DOI: 10.1128/jb.175.17.5529-5538.1993.
Warren T, Pasternak J Nucleic Acids Res. 1988; 16(22):10833-47.
PMID: 3205722 PMC: 338942. DOI: 10.1093/nar/16.22.10833.
Michaelis U, Schlapp T, Rodel G Mol Gen Genet. 1988; 214(2):263-70.
PMID: 3070350 DOI: 10.1007/BF00337720.
Chang Y, Cronan Jr J J Bacteriol. 1988; 170(9):3937-45.
PMID: 3045082 PMC: 211393. DOI: 10.1128/jb.170.9.3937-3945.1988.