I have a bacterial chromosome. I am struggling to understand how I can get protein sequences from genes annotated in the following way:
where sstart = 9 and send = 515
This is actually fine, as I can index the forward sequence using biopython:
faa = fasta[sstart:send].seq.translate(table=11)
However, when a gene is annotated in the reverse way, where sstart > send:
I am unsure how to get the corresponding protein sequence.
It would be much appreciated if someone could explain this.
Thanks for the clarification. That makes sense.
I also realised that the reason I was confused, is because I need to also figure out which frame the amino acid sequence is transcribed in, as the genes may not start with the right frame. Can you suggest what may be the best way to do this?