Open cgjosephlee opened 5 years ago
@cgjosephlee, thanks to reporting this issue, can I have the full sequence to work on this issue. It will be greatly appreciated.
Hi @cgjosephlee,
This branch issue_16 should fix this issue. At this point I only have a singularity container here containing the new code.
I have some other work to do before merging the code on the main branch, but if you want to test with the singularity version it will be appreciate.
Thanks for your comments and your help.
I'm using the latest mfannot via docker image.
I'm annotating a fungal mt genome with genetic code 4, and some of the ORFs were problematic. I have cross-validated with ncbi ORFfinder, mostly they are ORFs with non-ATG start codon.
e.q.
The start codon should be ATT (42151 bp) and stop codon should be TAA (44497 bp). The protein sequence in
.sqn
file isncbieaa "-ITFNKNIAFI...
.AAT is not a valid start codon and it should be ATA (56567 bp). The protein sequence in
.sqn
isncbieaa "-KIINTNMRI...
.I have got an error message:
It seems to be dealing with frame-shift annotations, I don't known if this is related to the problem.