Are there any online services or biopython scripts that will calculate the codon usage average for all sequences in a file not individually for each sequence. I only know python, unix and a little bit of perl.
7 answers
The EMBOSS program cusp takes one or more nucleotide sequences as input and outputs codon usage data, looking like this (first few lines):
#CdsCount: 1
#Coding GC 67.79%
#1st letter GC 67.88%
#2nd letter GC 46.89%
#3rd letter GC 88.60%
#Codon AA Fraction Frequency Number
GCA A 0.077 7.772 3
GCC A 0.462 46.632 18
GCG A 0.462 46.632 18
GCT A 0.000 0.000 0
....
There are a number of EMBOSS servers if you want to run the analysis online.
For the web services, you can find some services in the BioCatalogue.
You can then run those services with Taverna.
Hi Mr.Pierre
Can u guide me how to do codon usage biase - COA using codonW. I have downloaded source file of CodonW 1.4.2 and I used in window platform by offline.. I tried to to do.. I got some out put and I am not able to interpret. I need SCUO vs CAI SCUO vs GC1,2,3 ENC, but I am not able see...understand these thing in result.. can u plz guide me..
what is code I should type in window cmd.. to get these things
Thanks..
You could try http://www.bioinformatics.org/sms2/codon_usage.html and their codon usage tool. There is a mechanism to use the tool off-line. See link on the above page. I have no experience with this tool, but you may need to concatenate your individual sequences into one.
Leaving this here for the future reference... Following up on the concatenation approach, you could make use of the Bioconductor packages in R to do this: concatenate the sequences using Biostrings, and then analyse codon usage with coRdon.
https://bioconductor.org/packages/release/bioc/html/Biostrings.html
https://bioconductor.org/packages/release/bioc/html/coRdon.html
Log in to answer this question.
Not sure what is meant by "codon usage average". Perhaps you mean "frequency" ?