[Bioperl-l] Translation routine
Robson Francisco de Souza
rfsouza@citri.iq.usp.br
Wed, 24 Apr 2002 22:57:50 -0300 (BRT)
Hi,
I've found that translate function from PrimarySeqI apparently
is not tranlating bacterial start codons to 'M' when the start codon
is different from ATG. Please could someone take a look at my code
and data set (below) and see if I'm doing something wrong?
Thanks,
Robson
CODE:
use Bio::SeqIO;
$in = new
Bio::SeqIO(-file=>"/home/rfsouza/20020424.fasta",
-format=>"fasta");
$out = new Bio::SeqIO(-fh=>\*STDOUT,
-format=>fasta);
while ($s = $in->next_seq) {
print $o->write_seq($s->translate("*","X",0,11));
}
INPUT:
>NC_002488 Xylella fastidiosa 9a5c, gene:XF0677
TTGATGAGTACAACAAGTACGCGCCAGGGTATTCTGTCGCTCACCCTGAAAGACAGGGCG
GCGCTCTATAGCGCCTATATGATGTATATGCAAAATGGAGCGATTTTCGTCCCTACGCCC
AAACGTTATTTCCTTGGTAACGAAGTATTTTTGTTGCTGACCCTCCCTGATTCCAGCGAG
CGCTTACCGGTTGCTGGAAAAGTGGTCTGGTGTACTCCCGTTGGAGCGCAAGGCAATCGT
GTTGCGGGTATCGGTGTGCAACTGCCCGATGGTCCAGGAGGGGAAGTGGTGCGCAATCGC
ATTGAGACGTTGCTTGCCGGCATGTTCAACTCAGATAAGCCAACGCACACAATG
OUTPUT:
>NC_002488 Xylella fastidiosa 9a5c, gene:XF0677
LMSTTSTRQGILSLTLKDRAALYSAYMMYMQNGAIFVPTPKRYFLGNEVFLLLTLPDSSE
RLPVAGKVVWCTPVGAQGNRVAGIGVQLPDGPGGEVVRNRIETLLAGMFNSDKPTHTM
EXPECTED OUTPUT:
>NC_002488 Xylella fastidiosa 9a5c, gene:XF0677
MMSTTSTRQGILSLTLKDRAALYSAYMMYMQNGAIFVPTPKRYFLGNEVFLLLTLPDSSE
RLPVAGKVVWCTPVGAQGNRVAGIGVQLPDGPGGEVVRNRIETLLAGMFNSDKPTHTM