New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 120780523 (0.420)28 (0.570) 27/254 1.1 55.8% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi = 1209602 24 (0.490)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCT  >  NC_000913_3_bme_pgi/1209456‑1209749
                                                                                                                                                  |                                                                                                                                                   
cGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGAGACCAACAGGTAATGCAGAGcc                                                                                                                                                   >  2:887402‑M2/1‑147 (MQ=255)
        aGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCAGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGTGccttctccca                                                                                                                                           <  1:905151‑M2/149‑11 (MQ=255)
               ttatttGTCGGATAAGCCTTTGCAAAAACCGGGTATTCTCCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccacg                                                                                                                                    >  2:714200‑M2/1‑132 (MQ=255)
                         gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                          <  2:392590‑M2/149‑28 (MQ=255)
                           tAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                                                                                        <  1:505387/149‑1 (MQ=255)
                           tAAGCCTTTGCCAGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCTCAAACCAACGTTTATGAAAAt                                                                                                                        <  2:336974/149‑1 (MQ=255)
                              gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                                     >  1:789464‑M2/1‑117 (MQ=255)
                              gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                     <  1:433974/149‑1 (MQ=255)
                              gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAAAGCAGAGCCTTCTCCCAAAACAACGTTTATGAAAAAgaa                                                                                                                     >  1:308564/1‑149 (MQ=255)
                               ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                                                                                    <  1:927728‑M2/149‑34 (MQ=255)
                                cattccCGGTTACGGGTGTGCTTCAGCAGAAAAAACGGAACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                                   <  2:526523‑M2/147‑35 (MQ=255)
                                   tGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                                                                                                <  2:198377‑M2/149‑38 (MQ=255)
                                      cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCTGCCTGTTGGCGGAGTGGCTGAAGGCAACGGAGTCGGGACAGTAACAGGTAATGAAGAGccttctcccaacccagaggtttcgagagacgttgggatag                                                                                                             <  2:270317‑M2/149‑41 (MQ=255)
                                         ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTAGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgc                                                                                                          <  1:788909‑M2/149‑44 (MQ=255)
                                            cGGGTATTTTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCCTCTCCCCACCCAACGTTTATGAAAATGGAGAAATAACAAGCaa                                                                                                       >  2:541553/1‑149 (MQ=255)
                                             gggTATTCTTCAGCAGAAAAAGCGGTACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                                      >  2:783934‑M2/1‑102 (MQ=255)
                                             gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtgttcgagagccgtttgcacagtgccatc                                                                                                      <  1:678503‑M2/149‑48 (MQ=255)
                                                tATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                                                                   <  1:752958/149‑1 (MQ=255)
                                                  tgcttcagCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                                                                                                 >  2:65905‑M2/3‑97 (MQ=255)
                                                   tcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                                >  1:690873‑M2/1‑96 (MQ=255)
                                                     ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                                                                              <  1:823906‑M2/149‑56 (MQ=255)
                                                           agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGTTTAAGGCCACGGAACCGGGCCTCCAACAGGGAATGTAGAGCCTTCTCCCAAACCAATGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                        <  2:253660/149‑1 (MQ=255)
                                                           agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCCGGGccttcttccaggccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc                                                                                        >  2:871046‑M2/1‑88 (MQ=255)
                                                              aaaaGGGGCTCCAGTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCCTTcccactt                                                                                     <  2:430977/149‑6 (MQ=255)
                                                               aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaa                                                                                    >  1:336970‑M2/1‑84 (MQ=255)
                                                                      cACCATGGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                             <  1:184559/149‑1 (MQ=255)
                                                                       aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                                                                            <  2:67375‑M2/149‑74 (MQ=255)
                                                                       aCCATTGCATTTCAGCCAGCCTGTTAGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAGATAACAAGCAAATCGCATCATTCCTGCTTTTACCAggc                                                                            >  2:584072/1‑148 (MQ=255)
                                                                            ttcatcacaggcagcCTGGTGGCGGAGTGTCTTAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                                                                       <  2:298315/142‑1 (MQ=255)
                                                                                tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTGCTCCCAAACCAACGTTTATGAAAATGCAGAAAAAACAAGCAAATGGCAACATTCCGGCCTTTACCAgggggcgtacac                                                                   >  2:324649/1‑142 (MQ=255)
                                                                                      ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTGAACATGCtt                                                             >  1:629925/1‑149 (MQ=255)
                                                                                           cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcg                                                        <  1:802560‑M2/149‑94 (MQ=255)
                                                                                               tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaa                                                    >  1:392590‑M2/1‑52 (MQ=255)
                                                                                                     aGTGGCTGAAGGCTACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAATCAAATGGCATCATTCCTGCTTTTACCAGGGTGATTTAACATGCTTATTGGCTATGTAcgc                                              <  2:131458/149‑1 (MQ=255)
                                                                                                         gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgacgtgcggtaag                                          <  1:553856‑M2/149‑108 (MQ=255)
                                                                                                          cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtacgt                                         <  2:433328‑M2/149‑109 (MQ=255)
                                                                                                             aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgacgtgcggtaagtaac                                      <  1:762098‑M2/149‑112 (MQ=255)
                                                                                                               ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                 <  2:830699/136‑1 (MQ=255)
                                                                                                               ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                 >  1:830699/1‑136 (MQ=255)
                                                                                                                   aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGa                                <  2:273675/149‑1 (MQ=255)
                                                                                                                      gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgc                             <  1:422958‑M2/149‑121 (MQ=255)
                                                                                                                      gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                             >  1:481411/1‑149 (MQ=255)
                                                                                                                                acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                   >  1:251680/1‑149 (MQ=255)
                                                                                                                                  cAACAGGTAATGCAGAGccttctcccaaaccaacaaccaggagagccgtttgcacagtgccatccgatttgatagcgccaaacggattcttgcggcttaacagcagcgcgcgacgtgcggtaagtaactggtcgtgcgcccccgtctcg                 <  2:676919‑M2/149‑133 (MQ=255)
                                                                                                                                   aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTa                <  1:174109/149‑1 (MQ=255)
                                                                                                                                       ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                     >  2:271790/1‑140 (MQ=255)
                                                                                                                                       ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                     <  1:271790/140‑1 (MQ=255)
                                                                                                                                        gcaaTGCAGAGCCTTCTCCCAAACCAATGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCCTATCAACAAATGACCAGAACACAGATCTACAACg           <  1:163391/147‑1 (MQ=255)
                                                                                                                                        gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctg           >  2:823906‑M2/1‑11 (MQ=255)
                                                                                                                                        gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctg           <  2:905382‑M2/149‑139 (MQ=255)
                                                                                                                                            tGCAGAGtcttctctcaagcgaagtttttcgagagccgtttgcacagggccatccgatgtgatatcgccaaacggattcttgcggctgaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcac       >  2:981059‑M2/1‑7 (MQ=255)
                                                                                                                                             gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgggtgcttgcggcttaacagcagctcgcgaagtccggtaagtaactgctcgtgccgcccattctccacgctggcaca      <  2:255586‑M2/149‑144 (MQ=255)
                                                                                                                                               agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaac                                                                                   >  2:634718‑M2/1‑4 (MQ=255)
                                                                                                                                               agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaac                                                                                   <  1:634718‑M2/70‑67 (MQ=255)
                                                                                                                                                 agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTATTACCAGGGGGATTTAACATGCTTATTGGCTATGTAGGCGTATCAACAAATGACCAGAACACACGTCTACAACGTAATGCGCt  >  1:360264/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCT  >  NC_000913_3_bme_pgi/1209456‑1209749

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.