Uploaded image for project: 'Jalview'
  1. Jalview
  2. JAL-3771

CDS aligned in Aligned CDS/Protein view shown via 'show crossrefs->Uniprot' does not match protein alignment

    XMLWordPrintable

    Details

    • Type: Bug
    • Status: Open
    • Priority: Major
    • Resolution: Unresolved
    • Affects Version/s: 2.11.0, 2.11.1
    • Fix Version/s: 2.11.2
    • Component/s: na
    • Labels:
      None

      Description

      There's a subtle problem in the 'show cross-refs->Protein' workflow for a nucleotide alignment means that the CDS view doesn't include gaps inserted into the protein alignment due to misaligned codons in the original nucleotide alignment.

      e.g.
      >nucA
      AATGAG
      >nucB
      -ATG

      produces a protein alignment:
      >pA
      N-A
      >pB
      -M-

      When the protein alignment is shown (via show crossrefs->Uniprot for the case of Ensembl: ENSG00000113924), its corresponding CDS alignment is just the spliced exons from the original Ensembl locus view. Inserts reflecting the final aligned peptides need to be propagated back to the CDS prior to display.

       

      But the CDS reconstruction for this alignment should be:
      >nucA
      AAT---GAG
      >nucB
      ---ATG---

        Attachments

          Activity

            People

            Assignee:
            jprocter James Procter
            Reporter:
            jprocter James Procter
            Votes:
            0 Vote for this issue
            Watchers:
            1 Start watching this issue

              Dates

              Created:
              Updated: