Opened 11 years ago

Closed 5 years ago

Last modified 4 years ago

#145 closed defect (fixed)

bug when realign nucleic acid according to aligned protein

Reported by: guest Owned by: westram
Priority: normal Milestone: arb6.1
Component: no idea Version: arb-4+older
Keywords: Cc: syan@…

Description

Hello,

For a protein database, when there is "AA-" or "TGN" in the nucleic acid, it is translated to an "X" in aa sequence.

When I need to realign the nucleotide sequence according to the protein alignment (even I've just translated from the nucleic acid sequence to the aa sequence), it results either "Not a codon (not all IUPAC-combinations of 'AAT' translate and Not all IUPAC-combiations of 'GAT' translate)", or "Error while reading 'transl_table' (Illegal (or unsuppported) value (0) in 'transl_table' (species=…))".

I know that's much easier in Bioedit, since it can keep nucleotide codons and aa residues synchronized. I can easily align or change alignment of a nucleotide-based sequence in the protein mode. However, due to the realignment problem, it makes it not possible in arb.

My thought to solve this problem is: since in ARB, nucleotide and aa are in two different alignment fields, it might cause algorithmic problem concerning the gaps, dots or 'N's in nucleotide sequence. Then you can offer an option: to keep codon reading frames not changed. Since for most of the sequences, the user has corrected these Ns and gaps themselves manually, if only the protein alignment rather than the basic sequences are changed, the program can do the realignment without having to adjust de novo, where of the nucleotide correspond to a residue.

for example: "AAA AAT AAN GGT GG- —- GGA GGG GG." (the user is sure the reading frame should be arranged like that) will be translated as "L L X G X - G G X " When the user changed the protein alignment (but no sequence change made), he should be still able to realign the nucleotide according to the protein sequence.

Best greets and thanks!

Yan Shi (syan@…)

Change History (4)

comment:1 Changed 6 years ago by westram

  • Version changed from release_20071207 to oldVersion

comment:2 Changed 5 years ago by westram

  • Owner changed from devel to westram
  • Status changed from new to _started

comment:3 Changed 5 years ago by westram

  • Resolution set to fixed
  • Status changed from _started to closed

by [12667]

comment:4 Changed 4 years ago by westram

  • Milestone set to arb6.1

mark changes that got fixed after arb 6.0.x

Note: See TracTickets for help on using tickets.