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

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_bme_pgi/1209454‑1209748
                                                                                                                                                     |                                                                                                                                                  
cacGTAAATCAGGCAATTTATTTGTGGGATAAGCCTTTGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTACATTATTTCGGTCTTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCcgag                                                                                                                                                    >  1:188536/1‑149 (MQ=255)
 aCGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                                                                                                   >  1:248777‑M2/1‑148 (MQ=255)
    tAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                                                                                                                <  2:669279‑M2/149‑5 (MQ=255)
       aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                                             <  1:59371‑M2/149‑8 (MQ=255)
        tCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                                                                                                            <  1:574990‑M2/149‑9 (MQ=255)
            gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                                                                                                        <  2:426462‑M2/149‑13 (MQ=255)
                 ttatttGTCGGATAAGCCTTCGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATGGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                                                   >  2:427956‑M2/1‑132 (MQ=255)
                      tggcGGATAAGCCTTCGCTAGTTCCGGGTATTCTTCAGCCGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCTGAGTGGCTGAAGGCCACGGAACCGGAACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGtttt                                                                                                                              <  1:615766/146‑2 (MQ=255)
                          ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgag                                                                                                                          >  1:669279‑M2/1‑123 (MQ=255)
                              aaGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCGTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                                                                                      >  2:81930/1‑149 (MQ=255)
                                gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATGTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                                    >  2:187837‑M2/1‑117 (MQ=255)
                                gCCTTTGATAGTTCCGGGTATTCTCCAGAAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                                    >  2:488923‑M2/1‑117 (MQ=255)
                                  cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                                  >  1:478685‑M2/1‑115 (MQ=255)
                                  cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                                  >  2:501090‑M2/1‑115 (MQ=255)
                                       ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacccaggtcttcgagagccgtttgcaca                                                                                                             <  2:471457‑M2/149‑40 (MQ=255)
                                        cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                                                                                            >  1:345924/1‑149 (MQ=255)
                                           ttCCGGCTATTCTTCAGCAGAAAATGCGGCACCATTGCATTTCAGCCAGCCTTTTGGCGGAGTGGCTGAAGGCCACGGAATCGGGACGCCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                                                                                                         >  2:226901‑M2/1‑106 (MQ=255)
                                                 gTATTCTTCAGCAGAAAAAGAGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                   <  2:277370/149‑1 (MQ=255)
                                                  tATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccga                                                                                                  >  1:161894‑M2/1‑99 (MQ=255)
                                                    ttctccAGCAGAAAACGCGGCACCACTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                                                                                                >  2:311170‑M2/1‑97 (MQ=255)
                                                          agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat                                                                                          <  2:607839‑M2/149‑59 (MQ=255)
                                                           gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGGTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACGGGTAATGCCGGGccttctcccaaaccaagggtttcaagagccgttttcacgttgccatccgatttgatatc                                                                                         <  2:719538‑M2/149‑60 (MQ=255)
                                                                 aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaa                                                                                   <  1:805479‑M2/149‑66 (MQ=255)
                                                                      ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttcatatcgccaaacggct                                                                              <  1:128414‑M2/149‑71 (MQ=255)
                                                                         aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                           <  2:59371/149‑1 (MQ=255)
                                                                             ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcg                                                                       >  1:280729‑M2/1‑72 (MQ=255)
                                                                                   tgcagccagcctGTTGGCGGAGCGGTTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaatcgaagtttttcgagagccgtttgaaccgggccatccgatttgatatcgccaaacggattcttgcggcttaa                                                                 >  2:46324‑M2/3‑66 (MQ=255)
                                                                                     cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaca                                                               >  1:417082‑M2/1‑64 (MQ=255)
                                                                                     cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                               <  1:81930/149‑1 (MQ=255)
                                                                                           gccTGTTGGCGAAGAGGCTCGAGGCCACGGAACCGGGACAAATTGATGGTTAGCCGAGCCTTCCCCCCCGCCAACGTGTTTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                         >  2:577372/1‑149 (MQ=255)
                                                                                            ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                        >  1:412380/1‑149 (MQ=255)
                                                                                               gTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                     >  1:318242/1‑149 (MQ=255)
                                                                                                   gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagt                                                 >  1:236782‑M2/1‑50 (MQ=255)
                                                                                                     ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacaccagcgcgcgaagtgc                                               <  1:407202‑M2/149‑102 (MQ=255)
                                                                                                     ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCGATGTAc                                               >  2:64980/1‑149 (MQ=255)
                                                                                                      gggtgGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                              >  2:440607‑M2/3‑47 (MQ=255)
                                                                                                      gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                              >  1:607839‑M2/1‑47 (MQ=255)
                                                                                                        gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggt                                            >  1:401883‑M2/1‑45 (MQ=255)
                                                                                                           gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgcgatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaag                                         <  2:748001‑M2/149‑108 (MQ=255)
                                                                                                              gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaa                                      >  1:381930‑M2/1‑39 (MQ=255)
                                                                                                               aaGGACACGTAATCGGGACACCTACATGTAATGCAGAAccttctcccaaaccaaggtttccgagagccgtttgcatagtaccatccgaattgatatcgccaaacggattcttgccgcttaacagcagcgcgcgaagtgcggtaagtaac                                     >  2:264654‑M2/1‑38 (MQ=255)
                                                                                                                aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcctaacagcagcgcgcgaagtgcggtaagtacct                                    <  1:98136‑M2/149‑113 (MQ=255)
                                                                                                                 ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagagcgcgaactgcggtaagtaacta                                   <  2:38434‑M2/149‑114 (MQ=255)
                                                                                                                 ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTCCCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaacaa                                   >  2:317197/1‑149 (MQ=255)
                                                                                                                  gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgggttcttgcggcttaacagcagcgcgcgaggtgctgtatgtaactgg                                  <  1:169215‑M2/149‑115 (MQ=255)
                                                                                                                      cGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgt                              <  1:483159‑M2/149‑119 (MQ=255)
                                                                                                                        gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgc                            >  2:477408‑M2/1‑29 (MQ=255)
                                                                                                                             gggACACCACAAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc                       <  2:300858‑M2/149‑126 (MQ=255)
                                                                                                                                     aaCAGGTAATGCAGATccttcgcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcca               >  1:501566‑M2/1‑16 (MQ=255)
                                                                                                                                     aaCAGCTAATGCAGCGccttcgcccaaacccaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcca               >  1:294041‑M2/1‑16 (MQ=255)
                                                                                                                                      aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAGCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc              >  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)
                                                                                                                                                     |                                                                                                                                                  
CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  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.