Handle selenoprotein transcripts
Created by: mutalyzerbot
Original ticket: https://humgenprojects.lumc.nl/trac/mutalyzer/ticket/106 Original date: 2012/07/20 Original reporter: martijn
Currently Mutalyzer aborts translation of selenoprotein transcripts due to the appearent in-frame stop-codon.
In these transcripts, however, this codon does not stop translation and instead a selenocysteine is inserted. The GenBank format has a standard way of annotating this with /transl_except=(pos:X..Y,aa:Sec)
. There is also a way of annotating this in LRG references.
Attached is an example GenBank file.