[EMBOSS] Transeq question, frame phases
David Mathog
mathog at caltech.edu
Wed Feb 16 20:54:05 UTC 2011
Test case fasta file
>8Achars
AAAAAAAA
all 6 frames for transeq, standard mode emits:
>_1
KKX
>_2
KKX
>_3
KK
>_4
FF
>_5
FFX
>_6
FFX
But...
AAAAAAAA Forward
TTTTTTTT Reverse
abc cba <--- codons in diagram
^a^^b^^c^ phase 1 1 KKX 4 XFF
x^a^^b^^c^ phase 2 2 KKX 5 XFF
xx^a^^b^ phase 3 3 KK 6 FF
That is, frames 4->6 are supposed to use, respectively, the same
set of codons as 1->3, but translate on the opposite strand, shouldn't
the number of residues returned be like in the table above, also with
the X at the beginning rather than the end? Or to put it another way,
shouldn't the little "x" bases be ignored on the - strand if they were
also ignored on the +?
Thanks,
David Mathog
mathog at caltech.edu
Manager, Sequence Analysis Facility, Biology Division, Caltech
More information about the EMBOSS
mailing list