Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

show not allow setting translation start in the wrong direction of translation end #1838

Closed
khowe opened this issue Jan 4, 2018 · 1 comment
Assignees
Milestone

Comments

@khowe
Copy link

khowe commented Jan 4, 2018

Hi there,

Apollo lets you create transcripts with "backwards" CDSs. That is, for a forward strand transcript, you can set the translation end to a position upstream of the translation start. I've created one of these on the demo server:

ctgA:13251..28300 (15.05 Kb)
(annotation Apple3a-00001)

If you use the UI to fetch the cds/peptide sequence for this transcript, you get nothing back.

These annotations have a number of knock-on effects:

  • CDS features in the GFF3 having start > end (illegal GFF3)
  • REST IO/write method completely fails (for all sequences/transcripts) for cds and peptide FASTA dumps.

The problem similar to the second occurs also if there is no CDS for a protein-coding transcript (peptide dump failure for all transcripts). Although we're still trying to work out how our curators created CDS-less mRNAs!

A couple of suggested improvements:

  • Apollo could warn (or better, prevent) on attempt to create a CDS in the opposite orientation to the parent transcript
  • An "ignorebad" (or similar) option for the IO/write end-point that will return data for all valid transcripts, ignoring invalid ones
@nathandunn nathandunn added this to the 2.0.9 milestone Jan 4, 2018
@nathandunn nathandunn changed the title Backwards / missing CDSs cause problems for dumping show not allow setting translation start in the wrong direction of translation end Jan 25, 2018
@nathandunn nathandunn assigned nathandunn and unassigned deepakunni3 Jan 25, 2018
@nathandunn
Copy link
Contributor

@khowe I think it should not just not allow you to create a negative translation region, ever.

The strand is correct, but it has an impossible start / stop position:

ctgA	.	exon	18000	18800	.	+	.	Parent=8bcfd1c2-c4a6-4ca7-a41f-25e9aa5a1d9b;ID=2e542650-310e-42bf-9c9b-5182478826c8;Name=2e542650-310e-42bf-9c9b-5182478826c8
ctgA	.	CDS	18023	18013	.	+	0	Parent=8bcfd1c2-c4a6-4ca7-a41f-25e9aa5a1d9b;ID=eea9deba-b57b-4231-b727-a8c47427fe69;Name=eea9deba-b57b-4231-b727-a8c47427fe69

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants