New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 120780536 (0.910)15 (0.430) 14/254 1.8 35.7% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi = 1209602 22 (0.620)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATA  >  NC_000913_3_bme_pgi/1207657‑1207951
                                                                                                                                                    |                                                                                                                                                  
gCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGAAGCAGG                                                                                                                                                    >  2:105202‑M1/1‑149 (MQ=255)
  aTCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAg                                                                                                                                                  >  1:421279/1‑149 (MQ=255)
        tGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                                                                                                            >  2:336217/1‑149 (MQ=255)
          ccAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatta                                                                                                                                          <  1:570921‑M1/149‑11 (MQ=255)
              ccTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctg                                                                                                                                      >  1:275409‑M1/1‑135 (MQ=255)
               cTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattaattat                                                                                                                                     <  2:432855‑M1/149‑16 (MQ=255)
               cTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                                                                 <  1:199585/137‑1 (MQ=255)
               cTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                                                                                                 >  2:199585/1‑137 (MQ=255)
                tGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                                                    >  1:284606/1‑149 (MQ=255)
                     aaGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTCGCGATATCAAATCCGATGGCTCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGACGGctatgcattacctgttggtcc                                                                                                                               >  2:559379‑M1/1‑128 (MQ=255)
                           cggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGGGAGGGCTAc                                                                                                                         >  2:651414/1‑149 (MQ=255)
                               aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                                     >  1:470721/1‑149 (MQ=255)
                                     aGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                                                                                                               <  2:275409‑M1/149‑38 (MQ=255)
                                           ttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTCGGCGATATCAAATCGGATGGAACTGTGCAAAAGGCTATCGAAAACCTTGGTTTGGGAAAAGGAGCAAAACTCACTGCAGCCACGGCTACATTAGGACGCACCGGt                                                                                                         >  2:28899/1‑149 (MQ=21)
                                                 cACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                                                   <  2:54746/149‑1 (MQ=255)
                                                   ctgcgcgcgctgCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc                                                                                                 <  1:36640‑M1/146‑52 (MQ=255)
                                                       gcgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaaca                                                                                             <  1:63613‑M1/149‑56 (MQ=255)
                                                       gcgcgcTGCGGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGGACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtata                                                                                             <  1:296389/149‑1 (MQ=255)
                                                            gtgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGAAAACGAGCAATACTCCATGCTTCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                                                                        <  1:273713/148‑1 (MQ=25)
                                                             tgctgTTAAACCGCAAGAATCAGATTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCACGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGAACGCACCGGTTTCATAGCTATACCGGtt                                                                                       >  2:152314/1‑149 (MQ=37)
                                                                    aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaat                                                                                <  2:268025‑M1/149‑69 (MQ=255)
                                                                     aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                               <  2:470721/149‑1 (MQ=255)
                                                                         gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAGCGGCTATCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATCACGGTTATGGTTggtggt                                                                           >  1:339798/1‑149 (MQ=38)
                                                                         gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGAAGCAACTGCTACATTAGGACGCAACGGTTTCATAGCTATACCGGTTGTGATtggtgct                                                                           >  2:356516/1‑147 (MQ=255)
                                                                              aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc                                                                      <  1:544400‑M1/149‑79 (MQ=255)
                                                                                   gtatGGCGCTATTAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACGCAATGCAGCACCGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                                 <  2:98155/146‑1 (MQ=38)
                                                                                     ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa                                                               <  2:679435/149‑1 (MQ=255)
                                                                                     tgggcagtatCAAATCGTATGGCACTGGGCAAACGGCTCTCGAAAACCTGGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCa                                                               <  2:27345/142‑1 (MQ=255)
                                                                                          gctctCAAACCGGATGCCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGGGGTATTGAGCAATCAGTAAt                                                          >  2:327591/5‑149 (MQ=255)
                                                                                          gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTCTGATTGGTGGTAATGAGCACTCacacaa                                                          >  2:616550/1‑144 (MQ=255)
                                                                                           atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATAGTTGGTGGTATTGAGCAATCAGTATTc                                                         >  2:303185/1‑149 (MQ=255)
                                                                                                aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgacatgcaatggtgccgctttttctgctgacga                                                    >  2:264093‑M1/1‑53 (MQ=255)
                                                                                                   tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTg                                                 >  1:355828/1‑149 (MQ=255)
                                                                                                      gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGCTTTGGGAGAAGGAGCAAAACTCAATGCATCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                              <  1:614836/149‑1 (MQ=255)
                                                                                                             cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCGATGGAGGTGGAAt                                       >  2:282354/1‑149 (MQ=255)
                                                                                                                 gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCg                                   <  2:518068/149‑1 (MQ=255)
                                                                                                                  cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaagacactggtgccgctttttctcctcaagagaacccgggactgccaaa                                  >  2:23232‑M1/1‑35 (MQ=255)
                                                                                                                   aaaCGGCTCTCGAAAACCTTGGTTCGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATCAGGACGCCCCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCa                                 >  2:337108/1‑149 (MQ=255)
                                                                                                                       ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcgaaggctt                             >  1:126059‑M1/1‑30 (MQ=255)
                                                                                                                          tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatg                                                                          >  1:212912‑M1/1‑27 (MQ=255)
                                                                                                                          tctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatg                                                                          <  2:212912‑M1/101‑75 (MQ=255)
                                                                                                                          tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTACGCCATATTAGTGGTATTTCTATAAAGATTTTAAGTCTAGGTGTGGATCTTGGAGGTATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAt                          >  2:321597/1‑149 (MQ=255)
                                                                                                                            tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcgaaggcttatccg                        >  2:79302‑M1/1‑25 (MQ=255)
                                                                                                                             tgAAAACCTTGGTTTGGGAGAAGGAGCAAAACGCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTg                       <  2:213820/148‑1 (MQ=255)
                                                                                                                                    gttGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                >  1:92202/2‑149 (MQ=255)
                                                                                                                                     ttGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCACCGGTTACGTTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggg               <  1:303185/149‑1 (MQ=255)
                                                                                                                                        gTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattg            >  2:285137‑M1/1‑13 (MQ=255)
                                                                                                                                          ttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcc          <  1:81125‑M1/149‑139 (MQ=255)
                                                                                                                                          ctgGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcc          <  1:168298‑M1/148‑139 (MQ=255)
                                                                                                                                             ggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTg                                                 >  1:272625/1‑107 (MQ=255)
                                                                                                                                             ggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTg                                                 <  2:272625/107‑1 (MQ=255)
                                                                                                                                             ggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATgg       >  1:109130/1‑149 (MQ=255)
                                                                                                                                              gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttt                                                                <  2:271541‑M1/91‑85 (MQ=255)
                                                                                                                                              gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttt                                                                >  1:271541‑M1/1‑7 (MQ=255)
                                                                                                                                               atgtggctctgcattacctgttggggtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgatt     <  2:273584‑M1/145‑144 (MQ=255)
                                                                                                                                                  aGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATa  >  2:369322/1‑149 (MQ=255)
                                                                                                                                                   ggAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                                                                        <  1:286128/62‑1 (MQ=255)
                                                                                                                                                   ggAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                                                                        >  2:286128/1‑62 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATA  >  NC_000913_3_bme_pgi/1207657‑1207951

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.