New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 12078055 (0.110)36 (0.880) 32/252 0.3 81.4% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi = 1209602 12 (0.290)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_bme_pgi/1209454‑1209748
                                                                                                                                                     |                                                                                                                                                  
cacGTAAATCAGGCAATTTATTTGTGGGATAAGCCTTTGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTACATTATTTCGGTCTTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCcgag                                                                                                                                                    >  1:188536/1‑149 (MQ=255)
 aCGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGc                                                                                                                                                   >  1:248777‑M2/1‑148 (MQ=255)
    tAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGcctt                                                                                                                                                <  2:669279‑M2/149‑5 (MQ=255)
       aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctc                                                                                                                                             <  1:59371‑M2/149‑8 (MQ=255)
        tCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcc                                                                                                                                            <  1:574990‑M2/149‑9 (MQ=255)
            gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaa                                                                                                                                        <  2:426462‑M2/149‑13 (MQ=255)
                 ttatttGTCGGATAAGCCTTCGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATGGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaag                                                                                                                                   >  2:427956‑M2/1‑132 (MQ=255)
                      tggcGGATAAGCCTTCGCTAGTTCCGGGTATTCTTCAGCCGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCTGAGTGGCTGAAGGCCACGGAACCGGAACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAAGGtttt                                                                                                                              <  1:615766/146‑2 (MQ=255)
                          ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgag                                                                                                                          >  1:669279‑M2/1‑123 (MQ=255)
                              aaGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCGTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                                                                                      >  2:81930/1‑149 (MQ=255)
                                gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATGTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                                    >  2:187837‑M2/1‑117 (MQ=255)
                                gCCTTTGATAGTTCCGGGTATTCTCCAGAAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                                    >  2:488923‑M2/1‑117 (MQ=255)
                                  cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                                  >  1:478685‑M2/1‑115 (MQ=255)
                                  cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                                  >  2:501090‑M2/1‑115 (MQ=255)
                                       ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaacccaggtcttcgagagccgtttgcaca                                                                                                             <  2:471457‑M2/149‑40 (MQ=255)
                                        cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                                                                                            >  1:345924/1‑149 (MQ=255)
                                           ttCCGGCTATTCTTCAGCAGAAAATGCGGCACCATTGCATTTCAGCCAGCCTTTTGGCGGAGTGGCTGAAGGCCACGGAATCGGGACGCCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                                                                                                         >  2:226901‑M2/1‑106 (MQ=255)
                                                 gTATTCTTCAGCAGAAAAAGAGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                   <  2:277370/149‑1 (MQ=255)
                                                  tATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                                                                                                  >  1:161894‑M2/1‑99 (MQ=255)
                                                    ttctccAGCAGAAAACGCGGCACCACTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                                                                                                >  2:311170‑M2/1‑97 (MQ=255)
                                                          agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat                                                                                          <  2:607839‑M2/149‑59 (MQ=255)
                                                           gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGGTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACG‑GGTAATGCCGGGccttctcccaaaccaagggtttcaagagccgttttcacgttgccatccgatttgatatc                                                                                         <  2:719538‑M2/149‑60 (MQ=255)
                                                                 aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaa                                                                                   <  1:805479‑M2/149‑66 (MQ=255)
                                                                      ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttcatatcgccaaacggct                                                                              <  1:128414‑M2/149‑71 (MQ=255)
                                                                         aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                           <  2:59371/149‑1 (MQ=255)
                                                                             ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcg                                                                       >  1:280729‑M2/1‑72 (MQ=255)
                                                                                   tgcagccagcctGTTGGCGGAGCGGTTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaatcgaagtttttcgagagccgtttgaaccgggccatccgatttgatatcgccaaacggattcttgcggcttaa                                                                 >  2:46324‑M2/3‑66 (MQ=255)
                                                                                     cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaca                                                               >  1:417082‑M2/1‑64 (MQ=255)
                                                                                     cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                               <  1:81930/149‑1 (MQ=255)
                                                                                           gccTGTTGGCGAAGAGGCTCGAGGCCACGGAACCGGGACAAATTGATGGTTA‑GCCGAGCCTTCCCCCCCGCCAACGTGTTTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                         >  2:577372/1‑149 (MQ=255)
                                                                                            ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                        >  1:412380/1‑149 (MQ=255)
                                                                                               gTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                     >  1:318242/1‑149 (MQ=255)
                                                                                                   gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagt                                                 >  1:236782‑M2/1‑50 (MQ=255)
                                                                                                     ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacaccagcgcgcgaagtgc                                               <  1:407202‑M2/149‑102 (MQ=255)
                                                                                                     ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCGATGTAc                                               >  2:64980/1‑149 (MQ=255)
                                                                                                      gggtgGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                              >  2:440607‑M2/3‑47 (MQ=255)
                                                                                                      gAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                              >  1:607839‑M2/1‑47 (MQ=255)
                                                                                                        gTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggt                                            >  1:401883‑M2/1‑45 (MQ=255)
                                                                                                           gCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcgatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaag                                         <  2:748001‑M2/149‑108 (MQ=255)
                                                                                                              gAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaa                                      >  1:381930‑M2/1‑39 (MQ=255)
                                                                                                               aaGGACACGTAATCGGGACACCTACA‑TGTAATGCAGAAccttctcccaaaccaaggtttccgagagccgtttgcatagtaccatccgaattgatatcgccaaacggattcttgccgcttaacagcagcgcgcgaagtgcggtaagtaac                                     >  2:264654‑M2/1‑38 (MQ=255)
                                                                                                                aGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcctaacagcagcgcgcgaagtgcggtaagtacct                                    <  1:98136‑M2/149‑113 (MQ=255)
                                                                                                                 ggCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagagcgcgaactgcggtaagtaacta                                   <  2:38434‑M2/149‑114 (MQ=255)
                                                                                                                 ggCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTCCCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaacaa                                   >  2:317197/1‑149 (MQ=255)
                                                                                                                  gCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgggttcttgcggcttaacagcagcgcgcgaggtgctgtatgtaactgg                                  <  1:169215‑M2/149‑115 (MQ=255)
                                                                                                                      cGGAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgt                              <  1:483159‑M2/149‑119 (MQ=255)
                                                                                                                        gAACCGGGACACCAACA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgc                            >  2:477408‑M2/1‑29 (MQ=255)
                                                                                                                             gggACACCACAA‑GGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc                       <  2:300858‑M2/149‑126 (MQ=255)
                                                                                                                                     aaCA‑GGTAATGCAGATccttcgcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcca               >  1:501566‑M2/1‑16 (MQ=255)
                                                                                                                                     aaCA‑GCTAATGCAGCGccttcgcccaaacccaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcca               >  1:294041‑M2/1‑16 (MQ=255)
                                                                                                                                      aCA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAGCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc              >  1:65668/1‑149 (MQ=255)
                                                                                                                                            tAATGCAGAGccttctcgcaacccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg         >  2:405782‑M2/1‑10 (MQ=255)
                                                                                                                                            tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactgctcgtgccgccccttctccaggctgg         <  1:244846‑M2/149‑140 (MQ=255)
                                                                                                                                                   gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccgga  >  2:805479‑M2/1‑3 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACA‑GGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_bme_pgi/1209454‑1209748

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.