Skip to content
Steve Bond edited this page Oct 8, 2015 · 4 revisions

--count_codons, -cc

Description

Generate codon frequency tables.

Argument

concatenate ( str )

Optional. Pass in the word 'concatenate' (or a truncated version of the word, like 'concat' or 'con') to combine all sequences in the file to create a single frequency table.

Examples

Input file: Mnemiopsis_cds.fa

>Mle-Panxα1 cDNA - ML078817.
ATGTACTGGATATTTGAGATTTGTCAAGAGATAAAGCGAGCTCAATCCTGCCGAAAGTTC
GCGATAGACGGACCATTCGACTGGACGAACCGGATTATCATGCCAACACTCATGGTAATC
TGCTGCTTTCTCCAAACCTTCACCTTCATGTTCGGCAGCAACATCAGCTGTATCGGCTTC
GAGAAGTTGGAAAGGAACTTTGTGGAGGAGTACTGCTGGACCCAGGGTATCTATACAAGC
AAGGCTGCGTATAACATGCCATTACATACTCCCTACCCGGGGATTGCCCCCTGTGTGCCC
GAGTATGATCCCGTGACTCAGAAGTATTGGTTACCCTGTGGGGTGGAGGAAGAAGACAAG
GCTTATCATTTGTGGTATCAGTGGGTTCCGTTTTACTTTCTCGCTGTGGCCGTGGGTTAT
TATTTGCCATTTCTTATCTTGAAGGGTTCAAAGCTGCATCAGGTGAAGCCGCTGATTACG
TATTTGATGAACCAGAGGAACCTGGAGACTGATCCTAACCATTTGGTAGGAAAGCTATCG
CATTGGATCTTCAGACAGCTTGTTTATTCAAGGTTTGCGGCCACCTCTACAATCAGAATG
TACTGGCACGACTGGGGGCTTGTCCTCCTTGTTTGCTCTGTAAAGATCCTCTACCTTACC
GTCTCTCTTATCCACCTCTTTGCCACTGCCAAGATGTTCCACATCGGCAACTGGTTTACG
TACGGGATCATGTTCGCGCGGCGCAGCAACAGTCACACTACCCACGTTAAGGATGTGTTC
TTCCCGAAGATGGTGGCCTGTAAGATCGAGACATGGAGTTTCACAGGGAAGAATCATCTT
CACGGGATGTGTGTTTTAGCTCTGAACGTGATGAACCAATATTTGTTTTTGATCGTGTGG
TACGTCAACGTAATCATCATCTTCCTCAACAGTATCAGCTGTATTTACACTATAGTCAAG
TTTTGTAGCCCTAACATCGTTCACCACCGGATAGTCAACTCCTCCTCCTTAGACGACCAC
CATGACTTTACCCGGATGTTTGGTTATGTGGGACCTTCCGGACGGATAATCCTGGCTAAA
ATGTCGGAACATATGCCGGGGTACATGCTGAAACAAGTAGCCAAGAAGGTGACAGAGAAG
ATAGATATAGAGAATGAAAAGAATAGAGGGAGAGCCCCGACTATTAAGTTTACCAAGGTT
AACGGTCAGCCCTCAGAGCTGGCCAGACAGCCTCTCATGCACCTGAACGCTCTGATGTTA
GGGATGGTTCCTCAGAATCTACCAGAACCTAAAATTCAGAATATCCAACGGTCGCAGAAA
AAAGTACGGTTTCTGGTTTAA
>Mle-Panxα2 cDNA - ML25998a.
ATGGTATTGGATCTCATTTCTGGAAGCTTGAATGGCTTTTTAAAGATCAAGTCAGTTAGC
ATCGACGATCAGTGGGACCAGATTAACAGAACCTATTTGGTCATGTTTTGTATTTTATCT
GGTACAATCATGACCTTTAAACAGAATTTAGGATCAATAATACACTGTATATCGGATGCA
AGAGGCGACGACAGTTCGTTTGCGGATGCTCATGCGACATTTGTGCAAGACTATTGTGCT
GCTCAAGGGCTGTACACTTTAAAAGAAGTGTATGACAAGTCTTGGCCAGATGAAATTCCT
TACCCAGGTATTCTCCAAATGAAAACAATCGGTTGTTTCCCGGGGAGACAGTTCAAAAAC
GGAACCCCCATCCAGTGCCCGGACGAGAAAGATCTGAAACCCTTCACAACGGTCTATCAT
GTCTGGTACATGTTCGTACCGTTCTACTTCTGCGCTGTTGGCATCGCTTTTTACTTCCCC
TACACGGTTTTCAGACACCTCAGCGGCATCTACGACATCAAGCCTATGTTGAACAGCCTT
GCCCTCGACATTGGGGCCTACACGGAGGAGGACATAAGTCGACGTATAGACAATGTCTCG
AGGTGGTTGTACATCAAGTTGGATCCCTACATGAACAACATGCTTCCTTATACTCAGATA
GTTCACAAACATTCCATCTTTTACACGGTGATGTTGGTGAAGGTGATGTACCTAGCTACC
AGTGTTTCTATTTTTTACGCCACTCACCGGATATTCGACCAAGGAAACTTTGCACTCTAC
GGATACGATGTTCTAATGAGCATACCACAGGAAACAAGCTATAAAGTGATGGACACAATC
TTCCCTAAAATGGTTGGCTGTGAGATCAACATGTGGGGCCGGACTGGCGAACAGAGCGAA
TCTCTTCTGTGTGTCCTCCCTCAAAACATCGGCAACCAATACTTCTTCCTTATATTCTGG
TTTCTCCTGATTCTCACCATACTTTCCAACTGTATCTCTGTAATAGTGACCATATTCAGA
TTTATATTCGTTAGTGGGAGCTACAAAAGGTTCCTGGCTACCAGCCTCTTGAATCACGAA
GAACGATACAAGCTGGTGTTTACACATGTCGGCACGACTGGAAGATACATTTTACTGCTC
TGTGCCGATCATAGCAACCCCAAAATATTCGAGGATCTTCTAGAGATCGTCTGTTCCCTT
CTCATAGCAAACTATCACAAAAGAAAGAGGAGTCGGGATAAGGGACACAGTCGAGCGGAG
GGGGTAGGGACTAAAGGGCGACACGGTCTGTCTTTTGTGGACTCAACCGTGTGA

Usage example 1

$: sb Mnemiopsis_cds.fa -cc

Output

#### Mle-Panxα1 ####
Codon	AA	Num	Percent
AAA	K	5	1.119
AAC	N	17	3.803
AAG	K	23	5.145
AAT	N	5	1.119
ACA	T	6	1.342
ACC	T	8	1.79
ACG	T	3	0.671
ACT	T	7	1.566
AGA	R	5	1.119
AGC	S	6	1.342
AGG	R	3	0.671
AGT	S	3	0.671
ATA	I	8	1.79
ATC	I	21	4.698
ATG	M	19	4.251
ATT	I	7	1.566
CAA	Q	6	1.342
CAC	H	10	2.237
CAG	Q	11	2.461
CAT	H	8	1.79
CCA	P	5	1.119
CCC	P	6	1.342
CCG	P	6	1.342
CCT	P	6	1.342
CGA	R	2	0.447
CGC	R	1	0.224
CGG	R	7	1.566
CTA	L	2	0.447
CTC	L	8	1.79
CTG	L	10	2.237
CTT	L	7	1.566
GAA	E	6	1.342
GAC	D	7	1.566
GAG	E	12	2.685
GAT	D	4	0.895
GCC	A	9	2.013
GCG	A	4	0.895
GCT	A	7	1.566
GGA	G	4	0.895
GGC	G	3	0.671
GGG	G	9	2.013
GGT	G	5	1.119
GTA	V	6	1.342
GTC	V	5	1.119
GTG	V	13	2.908
GTT	V	9	2.013
TAA	*	1	0.224
TAC	Y	10	2.237
TAT	Y	12	2.685
TCA	S	3	0.671
TCC	S	5	1.119
TCG	S	3	0.671
TCT	S	3	0.671
TGC	C	5	1.119
TGG	W	12	2.685
TGT	C	8	1.79
TTA	L	5	1.119
TTC	F	13	2.908
TTG	L	8	1.79
TTT	F	15	3.356

#### Mle-Panxα2 ####
Codon	AA	Num	Percent
AAA	K	13	2.968
AAC	N	12	2.74
AAG	K	9	2.055
AAT	N	4	0.913
ACA	T	7	1.598
ACC	T	8	1.826
ACG	T	5	1.142
ACT	T	6	1.37
AGA	R	7	1.598
AGC	S	10	2.283
AGG	R	3	0.685
AGT	S	6	1.37
ATA	I	15	3.425
ATC	I	15	3.425
ATG	M	14	3.196
ATT	I	9	2.055
CAA	Q	6	1.37
CAC	H	8	1.826
CAG	Q	8	1.826
CAT	H	5	1.142
CCA	P	3	0.685
CCC	P	5	1.142
CCG	P	3	0.685
CCT	P	5	1.142
CGA	R	4	0.913
CGG	R	3	0.685
CGT	R	1	0.228
CTA	L	3	0.685
CTC	L	11	2.511
CTG	L	8	1.826
CTT	L	7	1.598
GAA	E	7	1.598
GAC	D	14	3.196
GAG	E	7	1.598
GAT	D	11	2.511
GCA	A	3	0.685
GCC	A	4	0.913
GCG	A	3	0.685
GCT	A	7	1.598
GGA	G	7	1.598
GGC	G	9	2.055
GGG	G	7	1.598
GGT	G	4	0.913
GTA	V	4	0.913
GTC	V	7	1.598
GTG	V	10	2.283
GTT	V	8	1.826
TAC	Y	19	4.338
TAT	Y	7	1.598
TCA	S	3	0.685
TCC	S	3	0.685
TCG	S	3	0.685
TCT	S	7	1.598
TGA	*	1	0.228
TGC	C	2	0.457
TGG	W	6	1.37
TGT	C	9	2.055
TTA	L	5	1.142
TTC	F	17	3.881
TTG	L	8	1.826
TTT	F	13	2.968

Usage example 2

$: sb Mnemiopsis_cds.fa -cc 'conc'

Output

#### concatination ####
Codon	AA	Num	Percent
AAA	K	18	2.034
AAC	N	29	3.277
AAG	K	32	3.616
AAT	N	9	1.017
ACA	T	13	1.469
ACC	T	16	1.808
ACG	T	8	0.904
ACT	T	13	1.469
AGA	R	12	1.356
AGC	S	16	1.808
AGG	R	6	0.678
AGT	S	9	1.017
ATA	I	23	2.599
ATC	I	36	4.068
ATG	M	33	3.729
ATT	I	16	1.808
CAA	Q	12	1.356
CAC	H	18	2.034
CAG	Q	19	2.147
CAT	H	13	1.469
CCA	P	8	0.904
CCC	P	11	1.243
CCG	P	9	1.017
CCT	P	11	1.243
CGA	R	6	0.678
CGC	R	1	0.113
CGG	R	10	1.13
CGT	R	1	0.113
CTA	L	5	0.565
CTC	L	19	2.147
CTG	L	18	2.034
CTT	L	14	1.582
GAA	E	13	1.469
GAC	D	21	2.373
GAG	E	19	2.147
GAT	D	15	1.695
GCA	A	3	0.339
GCC	A	13	1.469
GCG	A	7	0.791
GCT	A	14	1.582
GGA	G	11	1.243
GGC	G	12	1.356
GGG	G	16	1.808
GGT	G	9	1.017
GTA	V	10	1.13
GTC	V	12	1.356
GTG	V	23	2.599
GTT	V	17	1.921
TAA	*	1	0.113
TAC	Y	29	3.277
TAT	Y	19	2.147
TCA	S	6	0.678
TCC	S	8	0.904
TCG	S	6	0.678
TCT	S	10	1.13
TGA	*	1	0.113
TGC	C	7	0.791
TGG	W	18	2.034
TGT	C	17	1.921
TTA	L	10	1.13
TTC	F	30	3.39
TTG	L	16	1.808
TTT	F	28	3.164

Main Toolkit Pages





Further Reading

Clone this wiki locally