[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