CDS 64 - 612
    /note=Phagesdb Function: hits with cluster CT (20 counts) at 40% frequency
    /note=Phagesdb BLAST: top hits with Faja and Richie, e-value significant, call terminase small subunit
    /note=HHPRED: to hit terminase small subunit, e-value not significant
    /note=NCBI BLAST: top hits with Faja and Richie, e-value significant, call terminase small subunit
    /note=Pham Maps: compared with BillyTP, call terminase small subunit
    /note=
    /note=Start codon: original

CDS 599 - 1678
    /note=Phagesdb Function: hits with cluster S (17 counts) 24% frequency 
    /note=Phagesdb BLAST: top hits with Seahorse and Auxilium, e-value significant, call terminase large subunit ATPase domain
    /note=HHPRED: top hit terminase large subunit e-value significant
    /note=NCBI BLAST: top hits with Seahorse and Auxilium, e-value significant, call terminase large subunit ATPase domain
    /note=Pham Maps: compared to BillyTP matched terminase large subunit ATPase domain
    /note=
    /note=Start codon: original

CDS 1708 - 2037
    /note=Phagesdb Function: no data available 
    /note=Phagesdb BLAST: top hits with Attis and Bjanes7, e-value significant, function unknown 
    /note=HHPRED: top hit RNA polymerase I subunit A43, e-value not significant
    /note=NCBI BLAST: top hits with Actinomycetota bacterium and Arthrobacter sp. PsM3, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 2128 - 2829
    /note=Phagesdb Function: hits with cluster A1 (23 counts) 26% frequency 
    /note=Phagesdb BLAST: top hits with Vine and IDyl, e-value significant, call HNH endonuclease
    /note=HHPRED: top hit 
    /note=NCBI BLAST: top hits with Dermacoccus nishinomiyaensis and Rhodococcoides fascians, e-value significant, call HNH endonuclease
    /note=
    /note=Start codon: edited to position 2128 based on Z score and 1 to 1 hits in NCBI

CDS 2989 - 3453
    /note=Phagesdb Function: hits with cluster CR2, (11 counts) 16% frequency 
    /note=Phagesdb BLAST: top hits with Auxilium and BenchScraper, e-value significant, call terminase large subunit
    /note=HHPRED: top hit Terminase large subunit, endonuclease domain, e-value significant
    /note=NCBI BLAST: top hits with Raphaella and BenchScraper, e-value significant, call terminase large subunit nuclease domain
    /note=Pham Maps: compared to BillyTP matched terminase large subunit nuclease domain
    /note=
    /note=Start codon: original

CDS 3470 - 4807
    /note=Phagesdb Function: hit with cluster A1 (59 counts) 78% frequency
    /note=Phagesdb BLAST: top hits with Hestia and Auxilium, e-value significant, call portal protein
    /note=HHPRED: top hit portal protein, e-value significant
    /note=NCBI BLAST: top hits with Hestia and Auxilium, e-value significant, call portal protein
    /note=Pham Maps: compared with BillyTP matched portal protein
    /note=
    /note=Start codon: original

CDS 4810 - 6024
    /note=Phagesdb Function: hits with cluster F1 (30 counts) 50% frequency 
    /note=Phagesdb BLAST: top hits with MidnightRain and Faja, e-value significant, call capsid maturation protease
    /note=HHPRED: top hit Phage minor capsid protein 2, e-value significant, 
    /note=NCBI BLAST: top hits with MidnightRain and Faja, e-value significant, call capsid maturation protease
    /note=
    /note=Start Codon: Original

CDS 6147 - 6653
    /note=Phagesdb Function: hit with cluster L1 (19 counts) 28% frequency 
    /note=Phagesdb BLAST: top hits with Gorpy and Sakai, e-value significant, call scaffolding protein
    /note=HHPRED: top hit DUF/scaffolding protein, e-value significant
    /note=NCBI BLAST: top hits with Gorpy and Sakai, e-value significant, call scaffolding protein
    /note=Pham Maps: compared to BillyTP matched scaffolding protein
    /note=
    /note=Start codon: original

CDS 6685 - 7080
    /note=Phagesdb Function: hit with cluster CR2 15(counts) 22% frequency 
    /note=Phagesdb BLAST: top hits with Hestia and EvePickles, e-value significant, call capsid decoration protein
    /note=HHPRED: top hit lambda head decoration protein D, e-value significant
    /note=NCBI BLAST: top hits with Hestia and RadFad, e-value significant, call capsid decoration protein
    /note=Pham Maps: compared to BillyTP matched capsid decoration protein, LamD-like
    /note=
    /note=Start codon: original

CDS 7097 - 8137
    /note=Phagesdb Function: hit with cluster AY (17 counts) 26% frequency
    /note=Phagesdb BLAST: top hits with Faja and Gorpy, e-value significant, call major capsid protein
    /note=HHPRED: top hit major capsid protein, e-valye. significant
    /note=NCBI BLAST: top hits with Faja and Gorpy, e-value significant, call major capsid protein
    /note=Pham Maps: compared to BillyTP matched major capsid protein
    /note=
    /note=Start codon: original

CDS 8150 - 8488
    /note=Phagesdb Function: hits with clusters CR2 and AY (1 count) 50% frequency 
    /note=Phagesdb BLAST: top hits with BenchScraper and Gorpy, e-value significant, unknown function 
    /note=HHPRED: e-value not significant
    /note=NCBI BLAST: top hits with BillyTP and Auxilium, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 8505 - 8861
    /note=Phagesdb Function: hits with cluster A4 (64 counts) 79% frequency 
    /note=Phagesdb BLAST: top hits with Anekin and Persistence, e-value significant, call head-to-tail adaptor 
    /note=HHPRED: top hit with head-to-tail adaptor, e-value significant 
    /note=NCBI BLAST: top hits with Anekin and Persistence, e-value significant, call head-to-tail adaptor
    /note=Pham Maps: compared to BillyTP matches head-to-tail adaptor
    /note=
    /note=Start codon: original

CDS 8858 - 9256
    /note=Phagesdb Function: hits with cluster AY (16 counts) 57% frequency
    /note=Phagesdb BLAST: top hits with Anekin and Hestia, e-value significant, call head-to-tail stopper
    /note=HHPRED: top hit DUF, e-value significant
    /note=NCBI BLAST: top hits with Anekin and Hestia, e-value significant, call head-to-tail stopper
    /note=Pham Maps: compared to BillyTP matched head-to-tail stopper
    /note=
    /note=Start codon: original

CDS 9253 - 9636
    /note=Phagesdb Function: hits with cluster F1 (2 counts) 40% frequency 
    /note=Phagesdb BLAST: top hits with Hestia and Seahorse, e-value significant, call function unknown
    /note=HHPRED: top hit minor protein capsid, e-value significant
    /note=NCBI BLAST: top hits with Hestia and Seahorse, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 9633 - 10031
    /note=Phagesdb Function: hit with cluster AY (16 counts) 57% frequency
    /note=Phagesdb BLAST: top hits with Hestia and Seahorse, e-value significant, call tail terminator 
    /note=HHPRED: top hit with tail terminator protein, e-value significant 
    /note=NCBI BLAST: top hits with Hestia and Seahorse, e-value significant, call tail terminator 
    /note=Pham Maps: compared to BillyTP matched tail terminator
    /note=
    /note=Start codon: original

CDS 10079 - 10264
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hits with Hestia and Persistence, e-value significant, call unknown function 
    /note=HHPRED: to hit TraT complement resistance protein, e-value not significant 
    /note=NCBI BLAST: top hits with Hestia and Persistence, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 10257 - 10766
    /note=Phagesdb Function: hits with cluster EA1 (18 counts) 25% frequency 
    /note=Phagesdb BLAST: top hits with Hestia and ChuckDuck, e-value significant, call major tail protein 
    /note=HHPRED: hits with major tail protein, e-value not significant
    /note=NCBI BLAST: top hits with Hestia and ChuckDuck, e-value significant, call major tail protein
    /note=Pham Maps: compared to BillyTP match major tail protein
    /note=
    /note=Start codon: original

CDS 10868 - 11404
    /note=Phagesdb Function: hits with cluster AY (15 counts) 27% frequency 
    /note=Phagesdb BLAST: top hits with Seahorse and Hestia, e-value significant, call tail assembly chaperone
    /note=HHPRED: hit with tail assembly chaperone, e-value not significant
    /note=NCBI BLAST: top hits with Seahorse and Hestia, e-value significant, call tail assembly chaperone
    /note=Pham Maps: compared with BillyTP matched tail assembly chaperone
    /note=
    /note=Start codon: original

CDS 11491 - 11754
    /note=Phagesdb Function: hits with cluster AY (14 counts) 48% frequency 
    /note=Phagesdb BLAST: top hits with Hestia and Persistence, e-value significant, call tail assembly chaperone
    /note=HHPRED: top hit zinc finger protein, e-value not significant 
    /note=NCBI BLAST: top hits with Hestia and Persistence, e-value significant, call tail assembly chaperone
    /note=
    /note=Slippery sequence: CCCTTTC

CDS complement (11759 - 11974)
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: tops hits with Hestia and Seahorse, e-value significant, call unknown function 
    /note=HHPRED: top hit DUF, e-value not significant
    /note=NCBI BLAST: top hit Paenarthrobacter sp, e-value significant, call hypothetical protein

CDS 12043 - 16074
    /note=Phagesdb Function: hits with cluster EF (35 counts) 51% frequency 
    /note=Phagesdb BLAST: top hits with Seahorse and Auxilium, e-value significant, call tape measure protein 
    /note=HHPRED: top hit tape measure protein, e-value significant
    /note=NCBI BLAST: top hits with Seahorse and ChuckDuck, e-value significant, call tape measure protein
    /note=Pham Maps: compared to BillyTP matched tape measure protein

CDS 16074 - 16901
    /note=Phagesdb Function: hits with cluster AY (15 counts) 25% frequency 
    /note=Phagesdb BLAST: top hits with Hestia and Gorpy, e-value significant, call minor tail protein 
    /note=HHPRED: top hit baseplate centerpiece, e-value not significant 
    /note=NCBI BLAST: top hits with Hestia and Gorpy, e-value significant, call minor tail protein
    /note=Pham Maps: compared to BillyTP matched minor tail protein
    /note=
    /note=Start codon: original

CDS 16913 - 18214
    /note=Phagesdb Function: hits with cluster AY (15 counts) 25% frequency 
    /note=Phagesdb BLAST: top hits with Seahorse and Hestia, e-value significant, call minor tail protein 
    /note=HHPRED: top. hit baseplate centerpiece, e-value not significant
    /note=NCBI BLAST: top hits with Hestia and Gorpy, e-value significant, call minor tail protein
    /note=Pham Maps: compared to BillyTP matched minor tail protein
    /note=
    /note=Start codon: original

CDS 18223 - 19182
    /note=Phagesdb Function: hits with cluster AY (15 counts) 28% frequency 
    /note=Phagesdb BLAST: top hits with Seahorse and Auxilium, e-value significant, call minor tail protein 
    /note=HHPRED: top hit DUF, e-value not significant 
    /note=NCBI BLAST: top hits with Seahorse and Isolde, e-value significant, call minor tail protein
    /note=Pham Maps: compared to BillyTP matched minor tail protein
    /note=
    /note=Start codon: original

CDS 19182 - 19838
    /note=Phagesdb Function: hits with cluster AK (40 counts) 56% frequency 
    /note=Phagesdb BLAST: top hits with Seahorse and Eileen, e-value significant, call minor tail protein 
    /note=HHPRED: top hit DUF, e-value not significant 
    /note=NCBI BLAST: top hits with Seahorse and Eileen, e-value significant, call minor tail protein 
    /note=Pham Maps: compared to BillyTP matched minor tail protein
    /note=
    /note=Start codon: original

CDS 19849 - 20175
    /note=Phagesdb Function: hits with cluster O (7 counts) 78% frequency
    /note=Phagesdb BLAST: top hits with Anekin and Hestia, e-value significant, call function unknown 
    /note=HHPRED: top hit DUF, e-value not significant 
    /note=NCBI BLAST: top hits with Hestia and Seahorse, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original
    /note=
    /note=One membrane hit on DeepTMHMM; confirmed with SOSUI

CDS 20150 - 20383
    /note=Phagesdb Function: hits to cluster A9 (1 count). 100% frequency 
    /note=Phagesdb BLAST: top hits with Hestia and Seahorse, e-value significant, call function unknown 
    /note=HHPRED: top hit fullerene organizing protein, e-value not significant
    /note=NCBI BLAST: top hits with Hestia and Seahorse, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 20380 - 20691
    /note=Phagesdb Function: hits with cluster AO2 (6 counts) 86% frequency 
    /note=Phagesdb BLAST: top hits with Anekin and Seahorse, e-value significant, call function unknown
    /note=HHPRED: top hit DUF, e-value not significant
    /note=NCBI BLAST: top hits with Seahorse and Anekin, e-value significant, call hypothetical protein
    /note=*called membrane protein in NCBI for Anekin*
    /note=
    /note=Start codon: original
    /note=
    /note=3 red membrane peaks in TMHMM changed to membrane protein

CDS 20688 - 21674
    /note=Phagesdb Function: hits with cluster. AZ1 (20 counts) 31% frequency
    /note=Phagesdb BLAST: top hits with Seahorse and Anekin, e-value significant, Seahorse call lysin A, Anekin call endolysin
    /note=HHPRED: top hit membrane protein, e-value significant
    /note=NCBI BLAST: top hits with Seahorse and Anekin, e-value significant Seahorse call lysin (CDS notes peptidase), Anekin call endolysin
    /note=Pham Maps: compared to BillyTP matched endolysin
    /note=
    /note=Start codon: original

CDS 21686 - 22126
    /note=Phagesdb Function: hits with cluster EF (16 counts) 52% frequency
    /note=Phagesdb BLAST: top hits with Anekin and Seahorse, e-value significant, call function unknown
    /note=HHPRED: top hit holin, e-value not significant
    /note=NCBI BLAST: top hits with Seahorse and Anekin, e-value significant, call hypothetical protein
    /note=*CDS notes membrane protein*
    /note=3 red membrane peaks in TMHMM changed to membrane protein
    /note=
    /note=3 red membrane peaks in TMHMM changed to membrane protein
    /note=
    /note=Start codon: original

CDS 22188 - 22478
    /note=Phagesdb Function: hits with cluster A11 (5 counts) 33% frequency 
    /note=Phagesdb BLAST: top hits with RadFad and BenchScraper, e-value significant, call function unknown
    /note=HHPRED: e-value not significant
    /note=NCBI BLAST: top hits with RadFad and MidnightRain, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 22488 - 22922
    /note=Phagesdb Function: hits with cluster AU6 (4 counts) 21% frequency
    /note=Phagesdb BLAST: top hits with Persistence and Judy, e-value significant, call function unknown 
    /note=HHPRED: top hit with minor tail protein, e-value significant
    /note=NCBI BLAST: top hits with Judy and Persistence, e-value significant, call hypothetical protein 
    /note=*CDS notes minor tail protein/tail protein*
    /note=changed to minor tail protein with evidence from HHPRED and NCBI
    /note=Start codon: original

CDS 22926 - 23225
    /note=Phagesdb Function: hits with cluster EC (8 counts) 29% frequency 
    /note=Phagesdb BLAST: top hits with Seahorse and Zany, e-value significant, call function unknown/ minor tail protein 
    /note=HHPRED: top hit Rhs-family protein, e-value significant 
    /note=NCBI BLAST: top hits with NicSoilC5 and Seahorse, e-value significant, call hypothetical protein

CDS 23225 - 24667
    /note=Phagesdb Function: hits with cluster DV (5 counts) 21% frequency 
    /note=Phagesdb BLAST: top hits with Sarge and MacGully, e-value significant, call glycoside hydrolase/ minor tail protein
    /note=HHPRED: top hit glycoside hydrolase, e-value significant 
    /note=NCBI BLAST: top hits with Sphingomonas sp. and Microbacterium sp., e-value significant, call hypothetical protein

CDS complement (24723 - 26159)
    /note=Phagesdb Function: hits with cluster CA (33 counts) 54% frequency 
    /note=Phagesdb BLAST: top hits with Damascus and DizzyRudy, e-value not significant, call tape measure protein
    /note=HHPRED: top hit O-antigen polymererase, e-value significant
    /note=NCBI BLAST: top hit hypothetical protein, e-value significant, call hypothetical protein

CDS 26729 - 27340
    /note=Phagesdb Function: hits with cluster DE1 (16 counts) 22% frequency
    /note=Phagesdb BLAST: top hits with AxeJC and HFrancette, e-value significant, call HNH endonuclease
    /note=HHPRED: top hit endonuclease, e-value significant
    /note=NCBI BLAST: top hits with hypothetical protein and HNH endonuclease, e-value significant, call HNH endonuclease

CDS complement (26149 - 26730)
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hits with Beagle and Odyssey395, e-value significant, call function unknown 
    /note=HHPRED: hit protein, e-vlaue not significant, 
    /note=NCBI BLAST: no data
    /note=
    /note=Start codon: original (called by glimmer nothing to compare it to)

CDS 27344 - 27700
    /note=Phagesdb Function: hits with cluster CA (1 count) 100% frequency
    /note=Phagesdb BLAST: top hits with Frank15 and RadFad, e-value significant, call function unknown 
    /note=HHPRED: DUF, e-value not significant
    /note=NCBI BLAST: top hits hypothetical protein, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS complement (27689 - 27925)
    /note=Phagesdb Function: hits with cluster AY (8 counts) 13% frequency
    /note=Phagesdb BLAST: tops hit with Anekin, e-value significant, call helix-turn-helix DNA binding domain
    /note=HHPRED: DUF, e-value significant
    /note=NCBI BLAST: top hit helix-turn-helix DNA binding domain, e-value significant, call helix-turn-helix DNA binding domain
    /note=Pham Maps: compared to BillyTP matched helix-turn-helix DNA binding domain
    /note=*CDS notes HTH DNA binding protein*
    /note=
    /note=Start codon: start edited to 27925 based on 1 to 1 hits

CDS complement (27957 - 29135)
    /note=Phagesdb Function: hits with cluster K1 (21 counts) 28% frequency
    /note=Phagesdb BLAST: top hits with Seahorse and Anekin, e-value significant, call (tyrosine) integrase
    /note=HHPRED: hits with intergrase, e-value significant
    /note=NCBI BLAST: top hits with Seahorse and Anekin, e-value significant, call (tyrosine) integrase
    /note=Pham Maps:,comapred to BillyTP matche tyrosine integrase
    /note=*CDS notes Y-int*
    /note=
    /note=Start codon: original

CDS 29326 - 29664
    /note=Phagesdb Function: hits with cluster EA2 (2 counts) 67% frequency
    /note=Phagesdb BLAST: top hits with CherryTomatoes and Magpie, e-value not significant 
    /note=HHPRED: e-value not significant 
    /note=NCBI BLAST: no data
    /note=
    /note=Start codon: original (called by glimmer nothing to compare with)

CDS complement (29128 - 29685)
    /note=Phagesdb Function: hits with cluster BK1, (12 counts) 14% frequency
    /note=Phagesdb BLAST: top hits with Seahorse and Reed, e-value significant, call (metallo)phosphoesterase
    /note=HHPRED: top hits with metallophosphoesterase, e-value significant
    /note=NCBI BLAST: top hit with Seahorse, e-value significant,
    /note=
    /note=start edited to 29685 based on 1 to 1 hits, gene edit to metallo due to hits in HHPRED

CDS complement (29682 - 29888)
    /note=Phagesdb Function: hits with cluster C1 (1 count) 50% frequency
    /note=Phagesdb BLAST: top hits with Seahorse and EvePickles, e-value significant, call function unknown 
    /note=HHPRED: e-value not significant 
    /note=NCBI BLAST: top hits with Seahorse and EvePickles, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS complement (29885 - 30037)
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hit with Phrank15 and Sucker, e-value significant, call function unknown 
    /note=HHPRED: top hit biosynthesis protein, e-value not significant
    /note=NCBI BLAST: top hits with Phrank15 and Sucker, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS complement (30034 - 30585)
    /note=Phagesdb Function: hits with cluster CZ4 (9 counts) 53% frequency 
    /note=Phagesdb BLAST: top hits with Phrank15 and Seahorse, e-value significant, call function unknown 
    /note=HHPRED: e-value not significant
    /note=NCBI BLAST: top hits with Phrank15 and Seahorse, e-value significant, call hypothetical protein

CDS 31017 - 31184
    /note=Phagesdb Function: hits with cluster A1 (1 count) 100% frequency 
    /note=Phagesdb BLAST: top hits with BillyTP and Isolde, e-value significant, call function unknown 
    /note=HHPRED: top hit terminal domain, e-value not significant 
    /note=NCBI BLAST: top hits with BillyTP and Isolde, e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 31181 - 31543
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hits with BenchScraper and Richie, e-value significant, call function unknown 
    /note=HHPRED: e-value not significant 
    /note=NCBI BLAST: top hits with BenchScraper and Richie, e-value significant, call hypothetical protein 
    /note=Pham Maps: compared to BillyTP matched

CDS 31540 - 31731
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hit with Zucker, e-value significant, call function unknown 
    /note=HHPRED: top hits protein, e-value not significant 
    /note=NCBI BLAST: top hit with Zucker, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 31728 - 32027
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hits with TripleJ and Faja, e-value significant, call function unknown 
    /note=HHPRED: top hits protein, e-value not significant
    /note=NCBI BLAST:  top hits with TripleJ and Faja, e-value significant, call hypothetical protein

CDS 32024 - 32257
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hits with Anekin and Persistence, e-value significant, call function unknown 
    /note=HHPRED: top hit DNA glycosylase, e-value not significant
    /note=NCBI BLAST: top hits with Anekin and Faja, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 32254 - 32454
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hit MellowYellow, e-value significant, call function unknown 
    /note=HHPRED: top hit DUF, e-value not significant
    /note=NCBI BLAST: top hit Pseudarthrobacter sp, e-value significant, call hypothetical protein
    /note=
    /note=start edited to 32269 based on 1 to 1 target and guery hit in NCBI

CDS complement (32506 - 32880)
    /note=Phagesdb Function: hits cluster K3 (5 counts) 15% frequency
    /note=Phagesdb BLAST: top hits with Faja and Isolde, e-value significant, call helix-turn-helix DNA-binding domain protein
    /note=HHPRED: top hits hydrolase, e-value significant
    /note=NCBI BLAST: top hits with Isolde and Seahorse, e-valeu significant, call helix-turn-helix DNA-binding domain protein
    /note=*CDS notes metallo-protease*
    /note=
    /note=changed position to match with BillyTP

CDS complement (32897 - 33355)
    /note=Phagesdb Function: hits with cluster F1 (35 counts) 52% frequency
    /note=Phagesdb BLAST: top hits with MidnightRain and Anekin, e-value significant, call helix-turn-helix DNA binding domain protein
    /note=HHPRED: top hit DNA-binding domain, e-value not significant
    /note=NCBI BLAST: top hits with MidnightRain and Anekin, e-value significant, call helix-turn-helix DNA binding domain protein
    /note=
    /note=Start codon: original

CDS 33695 - 33883
    /note=Phagesdb Function: hits with cluster AK (22 counts) 31% frequency
    /note=Phagesdb BLAST: top hits with Phrank15 and BenchScraper, e-value significant, call helix-turn-helix DNA binding domain protein
    /note=HHPRED: top hit DNA-binding domain, e-value significant
    /note=NCBI BLAST: top hits with Phrank15 and RadFad, e-value significant, call helix-turn-helix DNA binding domain protein
    /note=Pham Maps: compared to BillyTP matched

CDS 33888 - 34112
    /note=Phagesdb Function: hits with cluster EE (28 counts) 43% frequency 
    /note=Phagesdb BLAST: top hits with Phrank15 and Anekin, e-value significant, call helix-turn-helix DNA binding domain protein
    /note=HHPRED: top hit DNA-binding domain, e-value significant 
    /note=NCBI BLAST: top hits with Phrank15 and EvePickles, e-value significant, call helix-turn-helix DNA binding domain protein
    /note=
    /note=Kept as helix-turn-helix DNA binding domain; tyrosine integrase present, but have not identified the immunity repressor
    /note=
    /note=tyrosine integrase and immunity repressor present; HHPRED hits with excisionase

CDS 34194 - 34430
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hits with Phrank15 and AbbyDaisy, e-value significant, call function unknown
    /note=HHPRED: DUF, e-value not significant
    /note=NCBI BLAST: top hit with Phrank15, e-value significant, call hypothetical protein

CDS complement (34166 - 34285)
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: no data, only draft for Bhageatrice 
    /note=HHPRED: hit hydrolase, e-value not significant
    /note=NCBI BLAST: no data
    /note=
    /note=Start codon: original (called by glimmer, nothing to compare with)

CDS 34427 - 34927
    /note=Phagesdb Function: hits with cluster, C1 (28 counts) 42% frequency 
    /note=Phagesdb BLAST: tops hits with Phrank15 and BillyTP, e-value significant, call HNH endonuclease
    /note=HHPRED: top hit homing endonuclease, e-value significant
    /note=NCBI BLAST: tops hit with Phrank15, e-value significant, call HNH endonuclease
    /note=Pham Maps: compared to BillyTP matched

CDS 35193 - 35471
    /note=Phagesdb Function: hits with cluster AY (7 counts) 39% frequency
    /note=Phagesdb BLAST: top hit with Phrank15, e-value significant, call helix-turn-helix DNA binding domain protein
    /note=HHPRED: top hit hypothetical protein, e-value not significant
    /note=NCBI BLAST: top hit with Phrank15, e-value significant, call helix-turn-helix DNA binding domain protein
    /note=Pham Maps: compared to BillyTP matched
    /note=*CDS notes as DNA binding domain protein also*
    /note=
    /note=Start codon: original

CDS 35468 - 35626
    /note=Phagesdb Function: hits with cluster AY (2 counts) 100% frequency
    /note=Phagesdb BLAST: top hits with Raphaella and BIllyTP, e-value significant, call function unknown 
    /note=HHPRED: DUF, e-value not significant
    /note=NCBI BLAST: top hits with Raphaella and Isolde, e-value significant, call membrane protein/oxidoreductase
    /note=Pham Maps: compared to BillyTP matched as hypothetical protein
    /note=
    /note=Start codon: original
    /note=
    /note=one large membrane hit with DeepTMHMM; confirmed by SOSUI

CDS 35623 - 36084
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hits with Raphaella and Richie, e-value significant, call function unknown 
    /note=HHPRED: top hit DUF, e-value significant
    /note=NCBI BLAST: top hit with Raphaella, e-value significant, call hypothetical protein
    /note=
    /note=start edited to 35623 more 1 to 1 hits with others off by 1-2

CDS 36081 - 36239
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hits with Raphaella and BillyTP, e-value significant, call function unknown
    /note=HHPRED: e-value not significant
    /note=NCBI BLAST: top hit with EvePickles, e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BillyTP matched

CDS 36239 - 36382
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hits with Auxilium and BillyTP, e-value significant, call function unknown
    /note=HHPRED: e-value not significant
    /note=NCBI BLAST: top hits with Auxilium and Hestia and other Arthrobacter phages , e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BIllyTP matched
    /note=
    /note=Start codon: original

CDS 36379 - 36504
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hit with RadFad, e-value significant, call function unknown
    /note=HHPRED: hit minor tail protein, e-value not significant
    /note=NCBI BLAST: top hit with RadFad, e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 36501 - 36650
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hit Seahorse, e-value significant, call function unknown 
    /note=HHPRED: hit uncharacterized protein, e-value not significant 
    /note=NCBI BLAST: top hit Seahorse, e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 36647 - 36949
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: no data only Bhageatrice draft
    /note=HHPRED: hit DUF, e-value not significant
    /note=NCBI BLAST: top hit Paenarthrobacter sp., e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 36946 - 37134
    /note=Phagesdb Function: hit cluster CF (1 count) 100% frequency
    /note=Phagesdb BLAST: top hits RadFad, e-value significant, call function unknown 
    /note=HHPRED: top hit Conserved nitrate reductase-associated protein, e-value not significant
    /note=NCBI BLAST: top hit RadFad, e-value significant, call hypothetical protein 
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 37131 - 37268
    /note=Phagesdb Function: no data
    /note=Phagesdb BLAST: top hit Seahorse, e-value significant, call function unknown 
    /note=HHPRED: top hit DUF, e-value significant
    /note=NCBI BLAST: top hit Seahorse, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 37265 - 37486
    /note=Phagesdb Function: hits cluster AY (3 counts) 19% frequency 
    /note=Phagesdb BLAST: top hits Seahorse and Auxilium, e-value significant, call function unknown (Seahorse)/ DNA binding protein (Auxilium)
    /note=HHPRED: e-value not significant
    /note=NCBI BLAST: top hits with Seahorse, e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 37483 - 37770
    /note=Phagesdb Function: hits with cluster F1 (1 count) 100% frequency 
    /note=Phagesdb BLAST: top hits Seahorse, e-value significant, call function unknown 
    /note=HHPRED: hits domain-containing protein, e-value not significant
    /note=NCBI BLAST: top hit Seahorse, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 37736 - 38218
    /note=Phagesdb Function: hits cluster C1 (2 counts) 25% frequency
    /note=Phagesdb BLAST: top hit Seahorse, e-value significant, call function unknown 
    /note=HHPRED: tits protein, e-value not significant
    /note=NCBI BLAST: top hit Seahorse, e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 38390 - 38560
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hits BenchScraper and BillyTP, e-value significant, call function unknown 
    /note=HHPRED: e-value not significant
    /note=NCBI BLAST: top hits Richie and BillyTP (other Arthrobacter phages), e-value significant, call hypothetical protein 
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 38560 - 38709
    /note=Phagesdb Function: non data 
    /note=Phagesdb BLAST: top hit BrayBeast, e-value significant, call function unknown 
    /note=HHPRED: hits DUF, e-value not significant
    /note=NCBI BLAST: top hit BrayBeast, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 38750 - 39430
    /note=Phagesdb Function: hits cluster N (16 counts) 22% frequency 
    /note=Phagesdb BLAST: top hits Isolde and Gorpy, e-value significant, call exonuclease
    /note=HHPRED: top hit exonuclease, e-value significant 
    /note=NCBI BLAST: top hits Isolde and Gorpy, e-value significant, call exonuclease
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 39430 - 39975
    /note=Phagesdb Function: hits cluster BP (8 counts) 89% frequency 
    /note=Phagesdb BLAST: top hits Gorpy and Sakai, e-value significant, call function unknown 
    /note=HHPRED: top hit binding protein, e-value not significant 
    /note=NCBI BLAST: op hits Gorpy and Sakai, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 39972 - 41192
    /note=Phagesdb Function: hits with cluster CR2 (20 counts)  25% frequency 
    /note=Phagesdb BLAST: top hits with Anekin and Zucker, e-value significant, call unknown function/AAA-ATPase
    /note=HHPRED: top hit DNA replication and repair protein, e-value significant 
    /note=NCBI BLAST: top hit with Arthrobacter sp, e-value significant, call AAA family ATPase
    /note=*DOUBLE CHECK*
    /note=
    /note=Start codon: original

CDS 41192 - 41716
    /note=Phagesdb Function: hits with cluster AY (6 counts) 9% frequency 
    /note=Phagesdb BLAST: top hits with Persistence and Anekin, e-value significant, call ssDNA binding protein(Persistence)/ SSB protein(Anekin) 
    /note=HHPRED: top hit single-stranded DNA-binding protein, e-value significant
    /note=NCBI BLAST: top hit with Persistence, e-value significant, call ssDNA binding protein
    /note=Pham Maps: compared to BillyTP matched as SSB protein
    /note=
    /note=Start codon: original

CDS 41825 - 42142
    /note=Phagesdb Function: hits with cluster N (29 counts) 44% frequency
    /note=Phagesdb BLAST: top hits with Phrank15 and EvePickles, e-value significant, call NrdH-like glutaredoxin
    /note=HHPRED: hits protein, e-value significant
    /note=NCBI BLAST: top hits with Phrank15 and EvePickles, e-value significant, call NrdH-like glutaredoxin
    /note=Pham Maps: compared to BillyTP and matched
    /note=
    /note=Start codon: original

CDS 42139 - 42372
    /note=Phagesdb Function: hits with cluster BB2 ( 2 counts) 100% frequency
    /note=Phagesdb BLAST: top hit with Persistence, e-value significant, call function unknown 
    /note=HHPRED: top hit hydrolase inhibitor, e-value not significant  
    /note=NCBI BLAST: top hit with Persistence, e-value significant, call hypothetical protein 
    /note=Pham Maps: compared to BillyTP matched hypothetical protein
    /note=
    /note=Start codon: original

CDS 42369 - 44897
    /note=Phagesdb Function: hits with cluster F1 (18 counts) 24% frequency 
    /note=Phagesdb BLAST: top hits with MidnightRain and Richie, e-value significant, call DNA helicase/methylase
    /note=HHPRED: top hit DNA binding protein, e-value significant 
    /note=NCBI BLAST: top hits with MidnightRain and Richie, e-value significant, call DNA helicase/methylase
    /note=Pham Maps: compared to BillyTP matched as DNA helicase/methylase
    /note=
    /note=Start codon: original

CDS 44894 - 45259
    /note=Phagesdb Function: hits with cluster AY (9 counts) 14% frequency 
    /note=Phagesdb BLAST: top hits with Auxilim and Phrank15, e-value significant, call helix-turn-helix DNA-binding (domain) protein
    /note=HHPRED: top hit DNA binding protein, e-value significant
    /note=NCBI BLAST:  top hits with Auxilim and Phrank15, e-value significant, call helix-turn-helix DNA-binding (domain) protein
    /note=Pham Maps: compared to BillyTP matched as call helix-turn-helix DNA-binding domain

CDS 45280 - 46362
    /note=Phagesdb Function: hits with cluster AR (18 counts) 24% frequency 
    /note=Phagesdb BLAST: top hits with RadFad and Auxilium, e-value significant, call DNA polymerase III sliding clamp beta/DNA polymerase III beta subunit
    /note=HHPRED: top hit with DNA polymerase III, beta subunit, e-value significant
    /note=NCBI BLAST: top hits with RadFad and Auxilium, e-value significant, call DNA polymerase III sliding clamp beta/DNA polymerase III beta subunit
    /note=Pham Maps: compared to BillyTP matched as DNA polymerase III sliding clamp beta
    /note=
    /note=Start codon: original

CDS 46359 - 46796
    /note=Phagesdb Function: hits with cluster DV (29 counts) 85% frequency 
    /note=Phagesdb BLAST: top hit with Auxilium, e-value significant, call function unknown 
    /note=HHPRED: hits with DUF, e-value not significant 
    /note=NCBI BLAST: top hit with Auxilium, e-value significant, call hypothetical protein 
    /note=Pham Maps: compared to BIllyTP matched
    /note=
    /note=Start codon: original

CDS 46793 - 48364
    /note=Phagesdb Function: hits with cluster CZ1 (19 counts) 26% frequency 
    /note=Phagesdb BLAST: top hits with Auxilium and Anekin, e-value significant, call DNA methylase/DNA methyltransferase
    /note=HHPRED: top hit with DNA methyltransferase, e-value significant
    /note=NCBI BLAST: top hits with Auxilium and Anekin, e-value significant, call DNA methylase/DNA methyltransferase
    /note=Pham Maps: compared to BillyTP matched as DNA methyltransferase

CDS 48381 - 48815
    /note=Phagesdb Function: hits with cluster AY (3 counts) 33% frequency 
    /note=Phagesdb BLAST: top hits with Auxilium and BenchScraper, e-value significant, call DNA binding protein/ function unknown 
    /note=HHPRED: top hit N-terminal phage replisome organiser, e-value significant 
    /note=NCBI BLAST: top hits with Auxilium and Seahorse, e-value significant, call DNA binding protein/ function unknown
    /note=Pham Maps: compared to BillyTP matched as hypothetical protein
    /note=
    /note=Start codon: original

CDS 49049 - 49420
    /note=Phagesdb Function: hits with cluster J (2 counts) 40% frequency 
    /note=Phagesdb BLAST: top hit with Phrank15, e-value significant, call function unknown 
    /note=HHPRED: top hit DUF, e-value not significant
    /note=NCBI BLAST: top hit with Phrank15, e-value significant, call hypothetical protein 
    /note=Pham Maps: compare to BillyTP matched as hypothetical protein
    /note=
    /note=Start codon: original

CDS complement (49000 - 49584)
    /note=Phagesdb Function: hits with cluster CZ4 (2 counts) 29% frequency 
    /note=Phagesdb BLAST: top hits with Che8 and Harley, e-value not significant, function unknown 
    /note=HHPRED: e-value not significant 
    /note=NCBI BLAST: no data
    /note=
    /note=Start codon: original (called by glimmer nothing to compare with)

CDS 49671 - 50240
    /note=Phagesdb Function: hits with cluster BD1 (6 counts) 46% frequency 
    /note=Phagesdb BLAST: top hits with Kukla and Isolde, e-value significant, call endonuclease/nuclease
    /note=HHPRED: top hit DNA repair protein, e-value significant 
    /note=NCBI BLAST: top hits with Kukla and Isolde, e-value significant, call endonuclease/nuclease
    /note=Pham Maps: compared to BillyTP matched as hypothetical protein

CDS 50237 - 50593
    /note=Phagesdb Function: hits with cluster AAY (1 count) 11% frequency 
    /note=Phagesdb BLAST: top hits with Auxilium and Faja, e-value significant, call function unknown 
    /note=HHPRED: hits with DUF, e-value not significant
    /note=NCBI BLAST: top hits with Auxilium and Isolde, e-value significant, call hypothetical protein

CDS 50590 - 51084
    /note=Phagesdb Function: hits with cluster AY (14 counts) 21% frequency 
    /note=Phagesdb BLAST: top hits with BenchScraper and Hestia, e-value significant, call MazG-like nucleotide pyrophosphohydrolase
    /note=HHPRED: hits with pyrophosphohydrolase, e-value significant
    /note=NCBI BLAST: top hits with BenchScraper and Hestia, e-value significant, call MazG-like nucleotide pyrophosphohydrolase
    /note=Pham Maps: compared to BillyTP matched 
    /note=*CDS notes pyrophosphatase*
    /note=
    /note=Start codon: original

CDS 51087 - 51266
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hit with BillyTP, e-value significant, call function unknown 
    /note=HHPRED: top hit zinc-ribbon domain, e-value significant 
    /note=NCBI BLAST: top hits with Richie and other Arthrobacter phage, e-value significant, call hypothetical protein 
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 51263 - 51688
    /note=Phagesdb Function: hits with cluster AY (18 counts) 30% frequency 
    /note=Phagesdb BLAST: top hits with Hestia and RadFad, e-value significant, call RusA-like resolvase
    /note=HHPRED: top hit Holliday junction resolvase RusA, e-value significant 
    /note=NCBI BLAST: top hits with Hestia and RadFad, e-value significant, call RusA-like resolvase
    /note=Pham Maps: compared to BillyTP matched as RusA-like resolvase
    /note=*CDS notes RusA-like resolvase (endonuclease)*
    /note=
    /note=Start codon: original

CDS 51685 - 51996
    /note=Phagesdb Function: hits with cluster BK1 (2 counts) 67% frequency 
    /note=Phagesdb BLAST: top hits with RadFad and BillyTP, e-value significant, call function unknown 
    /note=HHPRED: e-value not significant 
    /note=NCBI BLAST: top hits with Auxilium and RadFad, e-valu significant, call hypothetical protein 
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 52005 - 52238
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hits with Auxilium and BillyTP, e-value significant, call function unknown 
    /note=HHPRED: hits with DUF, e-value not significant 
    /note=NCBI BLAST: top hits with Auxilium and RadFad, e-value significant, call hypothetical protein 
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 52222 - 52449
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hits with Persistence and Auxilium, e-value significant, call function unknown 
    /note=HHPRED: top hit uncharacterized proteins, e-value not significant 
    /note=NCBI BLAST: top hits with Persistence and Isolde, e-value significant, call hypothetical protein
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: original

CDS 52446 - 52703
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hit with Auxilium, e-value significant, call function unknown 
    /note=HHPRED: top hit RNA-binding domain, e-value not significant
    /note=NCBI BLAST: top hit with Auxilium, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

CDS 52700 - 53323
    /note=Phagesdb Function: hits with cluster AY (7 counts) 24% frequency 
    /note=Phagesdb BLAST: top hits with Auxilium and Isolde, e-value significant, call helix-turn-helix DNA-binding (domain) protein
    /note=HHPRED: top hit DNA-binding domain, e-value significant 
    /note=NCBI BLAST: top hits with Auxilium and Isolde, e-value significant, call helix-turn-helix DNA-binding (domain) protein
    /note=Pham Maps: compared to BillyTP matched as helix-turn-helix DNA binding domain
    /note=
    /note=Start codon: original

CDS complement (53387 - 53668)
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: no data (all hits with drafts) 
    /note=HHPRED: top hit tertiapin, e-value not significant 
    /note=NCBI BLAST: no data
    /note=
    /note=Start codon: original (called by glimmer nothing to compare to)

CDS 53886 - 54041
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: top hits with Auxilium and Zucker, e-value significant, call function unknown 
    /note=HHPRED: e-value not significant 
    /note=NCBI BLAST: top hits with Auxilium and Zucker, e-value significant, call hypothetical protein 
    /note=Pham Maps: compared to BillyTP matched
    /note=
    /note=Start codon: kept start the same 1 to 1 match for target and query for majority of matches

CDS complement (53785 - 53916)
    /note=Phagesdb Function: no data 
    /note=Phagesdb BLAST: no data (only Bhageatrice draft) 
    /note=HHPRED: e-value not significant 
    /note=NCBI BLAST: no data
    /note=
    /note=Start codon: original (called by glimmer nothing to compare to)

CDS 54051 - 54281
    /note=Phagesdb Function: hit with cluster K1 (6 counts) 46% frequency 
    /note=Phagesdb BLAST: top hits with Auxilium and AbbyDaisy, e-value significant, call function unknown 
    /note=HHPRED: top hit signaling protein, e-value not significant 
    /note=NCBI BLAST: top hits with Auxilium and AbbyDaisy, e-value significant, call hypothetical protein
    /note=
    /note=Start codon: original

Genes added in gaps (Only added in the DNA Master Minimal file):
•	SEA_BHAGEATRICE_70 (38238-38393): Hypothetical protein and hits with Phrank15 and other AY phages
•	SEA_BHAGEATRICE_96 (53525-53734): (forward) in the minimal file gets significant hits with cluster AY phages (Auxilium and Windest).

Deleted phams in the DNA MASTER minimal file (Genes not deleted in PECAAN):
•	Phams 41, 57, 87, 98, 100 as annotated in PECAAN; removed the two reverse genes (ORPHAMs) at the end of the genome. Most of the annotated cluster AY phages lack reverse genes at the end of the genome.