home consulting training resources documentation documentation links links admin & aux services admin & aux services

Equivalent Applications for GCG programs

GCG on the SACS system will be retired in Fall 2010. Most of the GCG programs are available in EMBOSS. This table lists the EMBOSS equivalent for each GCG program, or where applicable, an alternative program that is available on the SACS System.

Full list of all EMBOSS programs.

GCG program EMBOSS program Description/Comments
Assemble merger
Construct new sequences from pieces of existing sequences; merger only accepts 2 sequences while assemble and union accept several.
BackTranslate backtranseq
Backtranslate protein -> nucleotide sequence.
backtranambig backtranslates to ambiguous codons.
BestFit water
Bestfit uses the Smith-Waterman algorithm to find the best local alignment between 2 sequences. water uses Smith-Waterman, matcher uses Pearson's lalign algorithm.
dbiBlast NCBI homology search between query and database
Breakup splitter Splits a sequence into (overlapping) smaller sequences
Chopup - Helps to convert a non-GCG sequence format
Not needed in EMBOSS because it reads most sequence formats without conversion
CodonFrequency chips
CodonFrequency --tabulates codon usage.
compseq -- counts composition of dimer/trimer in sequence.
chips -- calculates codon usage stats
cusp -- creates a codon usage table.
CodonPreference syco
Recognize protein coding sequences
CoilScan pepcoil Predicts coiled-coil regions
Compare + DotPlot dottup + dotmatcher
2-sequence comparison.
dotpath does a non-overlapping wordmatch dotplot.
Composition compseq
Sequence composition
compresstext - Removes extra whitespace in text files. Can be done via Unix shell script.
comptable - Creates a scoring matrix
consensus prophecy Creates a consensus sequence
or matrices/profiles from multiple alignments
correspond codcmp Codon usage table comparison
corrupt msbar Randomly mutate sequence
dataset dbiflat
Creates searchable sequence database. GCG's Dataset requires sequences in GCG format, whereas dbiflat, dbiblast, dbigcg will take most formats between them.
detab - Replaces tabs with spaces in sequence files. Can be performed by Unix shell command.
distances - Calculates pairwise evolutionary distances between aligned sequences. The Phylip package can do this.
diverge - Estimates pairwise substitutions per site between 2 or more coding sequences. The Phylip package can do this.
dotplot dottup
2-sequence comparison
extractpeptide transeq ExtractPeptide takes the output of Map and can write one or more of the reading-frame translations. transeq translates one or more of the frames or specific regions directly from an input nucleotide sequence.
- Pearson's homology-search programs.
fetch seqret
Pull one or more sequences out of the databases. seqret/seqretsplit can save output in various sequence formats.
figure - Generates plots from other GCG programs. The equivalent EMBOSS programs usually generate plots (e.g. plotorf).
findpatterns fuzznuc
searches for patterns in a sequence or database
fingerprint - Finds the products of T1 ribonuclease digestion.
fitconsensus - Use after Consensus to find the best fits.
framealign - Finds best local alignment including frame shifts between a protein and nucleotide sequence.
frames plotorf
Show open reading frames. plotorf does this graphically
framesearch - Homology searches including frameshifts between protein and nucleotide sequences
- Converts from various formats to GCG sequence format. Unnecessary in EMBOSS because it can accept most sequence formats, but seqret can convert between formats if desired.
Gap needle
Needleman-Wunsch algorithm to compare 2 sequences. stretcher uses the Myers-Miller algorithm which is more memory-efficient. For sequences larger than 10kb, I would suggest you to use 'stretcher' program in EMBOSS which is also a global alignment program. If one of your sequence is genomic and you are trying to align an est sequence to it, you may want to consider the 'est2genome' program. On the other hand, water->matcher->supermatcher are local alignment programs for small, medium, and large sequences, respectively.
Gapshow plotcon Graphical representation of similarity of 2 sequences.
GCGtoBlast - Makes a Blast database. Use NCBI's 'formatdb' instead.
Parts of GCG's gel assembly suite.
GetSeq seqret Type in a new sequence
GrowTree - Creates phylogenetic tree. Can use Phylip or Clustal instead.
HelicalWheel pepwheel Plots peptide sequence as helical wheel to help recognize amphiphilic regions.
- Sean Eddy's HMMER package.
HTHScan helixturnhelix Finds HTH motifs in protein sequences.
IsoElectric iep Calculates isoelectric pt of protein.
Lineup - Edits multiple sequence alignments
ListFile - for printing. Can use Unix pcprint command instead.
Lookup - Versatile program for finding sequences in a database. whichdb in emboss can search for accession numbers, but GCG's lookup is much more sophisticated. Use NCBI Entrez instead.
finds restriction enzyme cleavage sites.
GCG & EMBOSS may display different isoschizomers of the same enzyme, but the results are equivalent. The EMBOSS remap program may not display a few of the available isoschizomers.
MeltTemp dan Computes melting temperature of oligos
MEME - Finds conserved motifs in a group of unaligned sequences.
MFold - Predicts nucleotide secondary structure. GCG's version is an old version of Zuker's MFOLD.
Moment pepnet,octanol
Makes a contour plot of the helical hydrophobic moment of a peptide sequence
hmoment prints the text output of the calculation.
Motifs patmatmotifs Finds common Prosite motifs in a sequence. Use '-full' tag to display abstract information when using EMBOSS patmatmotifs. Note that both these programs will only find Prosite 'Patterns' (e.g. CAMP Phosphorylation Site), and not Prosite 'Matrices' (e.g. Helix-turn-Helix). Use Interproscan to find all known domains and functional sites. (http://www.ebi.ac.uk/InterProScan/).
patmatmotifs can accept file containing multiple sequences or patterns.
Meme + Motifsearch prophecy + profit Search a sequence or database with a matrix or profile.
Names infoseq provides some info about sequence specifications.
- remote access to NCBI's Blast.
NoOverlap diffseq Finds differences between 2 sequences. NoOverlap can work with a group of sequences.
OldDistances - Makes a table of the pairwise similarities within a group of sequences.
onecase - converts sequence into lower or upper case. Can be performed by Unix shell command.
Overlap - Compares 2 sets of sequences using Wilbur-Lipman algorithm.
Paupdisplay + Paupsearch - PAUP Phylogenetic Analysis./td>
Pepdata getorf
Translates in all 6 reading frames. sixpack displays the DNA sequence with 6-frame translations and orfs.
Pepplot pepinfo Pepplot plots protein 2ndary structure and hydrophobicity. pepinfo plots hydrophobicity, and garnier does protein 2ndary structure prediction.
Peptidemap digest Enzyme/reagent cleavage map of a protein.
Peptidesort digest
GCG peptidesort sorts fragments from an enzyme/reagent cleavage of one or more proteins according to position, mol. wt., and HPLC retention. EMBOSS digest only processes one reagent cleavage at a time. EMBOSS pepstats can be used to determine the composition of the fragments afterwards.
The EMBOSS programs do not provide the elution times from HPLC. If you need this data, try the UCSF MS-Digest program which has an option for HPLC Indices.
Secondary structure prediction. Garnier does not include Jameson-Wolf antigenic indexing.
antigenic predicts potentially antigenic regions of a protein sequence, using the method of Kolaskar and Tongaonkar.
pepwindow displays Kyte-Doolittle protein hydropathy.
pepwindowall produces a set of superimposed Kyte & Doolittle hydropathy plots from an aligned set of protein sequences.
Pileup emma Multiple sequence alignment. emma is an interface to ClustalW.
PlasmidMap cirdna
Plot DNA constructs.
PlotFold - Plots MFold output.
PlotSimilarity plotcon Graphical representation of the similarity along a set of aligned sequences.
Calculates consensus sequence from a multiple sequence alignment, and displays them prettily.
Prime eprimer3 Selects oligonucleotide primers.
Creates matrices/profiles from multiple alignments. Gapped alignment for profiles and sequences.
PrimePair primersearch? Evaluates individual primers to determine their compatibility for use as PCR primer pairs.
Profilescan patmatdb Searches sequences or db for protein motifs. Profilescan uses Gribskov method.
Profilesearch profit Scans a sequence or database with a matrix or profile.
Profilesegments - Alignments for results of Profilesearch
Publish seqret
Makes publication-quality displays of sequences.
Reformat seqret GCG requires input sequences to be in GCG format, hence other formats need to be converted with 'reformat'. Emboss programs accept most sequence formats, so conversion is rarely required, but 'seqret' can be used to convert between formats if desired.
Repeat equicktandem

Finds tandem repeats in sequences. The equivalent group of Emboss programs will also look for inverted or palindromic repeats.
Replace biosed
Replaces characters in a text file. Degapseq is specific for replacing gap characters. Can be performed with Unix shell utilities like sed, awk or tr.
Reverse revseq Reverse/complement a sequence.
Sample extractseq Extract regions from a sequence.
Seg maskseq Masks off low-complexity regions from a sequence.
Seqed biosed, cutseq, degapseq, descseq, entret, extractfeat, extractseq, listor, maskfeat, maskseq, newseq, noreturn, notseq, nthseq, pasteseq, revseq, seqret, seqretsplit, skipseq, splitter, trimest, trimseq, union, vectorstrip, yank Sequence editor. EMBOSS has several tools for specific editing tasks. Or use a text editor (not word processor!).

Try the Jemboss alignment editor for editing multiple sequence alignments.
Other alternatives are BioEdit (Windows only) and Seaview (Mac, Windows, Unix).

SeqLab - X-windows interface to GCG.
Setkeys - Redefines keyboard keys, mainly used for GCG's gel assembly programs.
Shiftover - Moves text by column. Use the nedit editor instead.
Shuffle shuffleseq Shuffles a sequence.
Simplify - Reduce the number of symbols in a sequence.
Spew - Sends a sequence from a remote computer to your desktop. Use one of the File transfer mechanisms instead.
SPScan sigcleave Predicts signal peptides in protein sequences.
Ssearch - Part of Pearson's Fasta package.
StatPlot - Plotting program. Rarely used.
StemLoop palindrome
Finds inverted repeats.
Stringsearch textsearch Finds text phrases in sequence or database. Use NCBI's Entrez instead.
Terminator - searches for prokaryotic factor-independent RNA polymerase terminators according to the method of Brendel and Trifonov.
Testcode wobble Plots 3rd-position variability as an indicator of potential coding regions.
seqret Emboss accepts most sequence formats, therefore format conversion is rarely required. seqret can be used to convert between formats if desired.
Translate transeq Translates nucleotide -> Protein sequences
Transmem - predicts transmembrane helices.
Window + Statplot freak Residue/base frequency table or plot.
- Homology search using Wilbur/Lipman algorithm. Segments displays the result.
Xnu - Masks tandem repeats for future Blast search.
- abiview Reads ABI file and displays trace
- antigenic Finds antigenic sites in proteins
- banana Bending and curvature plot in B-DNA
- btwisted Calculates the twisting in a B-DNA sequence
- cai CAI codon adaptation index, to measure synonymous codon usage bias.
- chaos Create a chaos game representation plot for a sequence
- charge Protein charge plot.
- checktrans Reports STOP codons and ORF statistics of a protein
- coderet Extract CDS, mRNA and translations from feature tables
- cpgplot
Plots and reports CpG-rich regions.
seqed cutseq Removes a specified section from a sequence. seqed is interactive, cutseq is command-line.
seqed degapseq Alter name/description of sequence.
Findpatterns dreg Regular expression search of a sequence. Findpatterns is an approximate equivalent.
- emma interface to ClustalW program.
- emowse Protein identification by Mass spectrometry.
- epestfind Finds PEST motifs as potential proteolytic cleavage sites
- est2genome Align EST and genomic DNA sequences.
- extractfeat Extract features from a sequence.
- findkm Find Km and Vmax for an enzyme reaction by a Hanes/Woolf plot
- fuzztran Protein pattern search after translation
- geecee Calculates the fractional GC content of nucleic acid sequences
- isochore Plots isochores in large DNA sequences
- listor Writes a list file of the logical OR of two sets of sequences
- makenucseq
Create random nucleotide and protein sequences
- marscan Finds MAR/SAR sites in nucleic sequences
- maskfeat Mask off features of a sequence.
- mwcontam Shows molwts that match across a set of files
- mwfilter Filter noisy molwts from mass spec output
- noreturn remove carriage return from a ASCII files. Can be performed by Unix utilities like 'tr'.
Reformat nthseq Pulls one sequence out of a multiple set. Reformat will pull a sequence out of an MSF or RSF file.
- oddcomp Finds protein sequence regions with a biased composition
- polydot Displays all-against-all dotplots of a set of sequences
- printsextract Extract data from PRINTS
- pscan Scans proteins using PRINTS
- rebaseextract
Search and extract from REBASE.
- recoder Remove restriction sites but maintain the same translation
- seqmatchall all-against-all comparison of a set of sequences.
- showdb Shows info about currently available databases.
- showfeat Shows features of a sequence
- silent Silent mutation restriction enzyme scan
- sirna Finds siRNA duplexes in mRNA
- stssearch Searches a DNA database for matches with a set of STS primers
- supermatcher Finds a match of a large sequence against one or more sequences
- tfextract Extract data from TRANSFAC database.
gcghelp tfm shows documentation for a program.
- tfscan Scans DNA sequences for transcription factors
- tmap Displays membrane spanning regions
- tranalign Align nucleic coding regions given the aligned proteins
- trimest
Trim bits off ends of sequences. Can be done interactively with GCG's seqed.
- twofeat inds neighbouring pairs of features in sequences
- vectorstrip Strips out DNA between a pair of vector sequences
- wordcount Counts words of a specified size in a DNA sequence
- wordmatch Finds all exact matches of a given size between 2 sequences

Modified NIH page to meet SACS requirements

last modified 12/16/2010

Consulting | Training | Resources | Documentation | Links | Admin & Aux services