Hello,

We are working on annotating phage Pochaccho (Cluster GK) and we're uncertain as to whether we should keep a gene call made by Glimmer but not GeneMark. Glimmer predicts a gene starting at 16,587 to 16,730. This would give an overlap of 8 nucleotides with the previous gene and a total size of 111 bp. Given that it is less than 120 bp we were about to delete it but there's coding potential and DeepTMHMM shows that there's a transmembrane domain present is there a chance that this is a potential real gene? The gene is an orpham so there's nothing else we can compare it to.

Thank you for your advice.