New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1207790 =4 (0.130)17 (0.640) 16/254 0.7 81.8% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi 1209619 = 4 (0.150)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGG  >  NC_000913_3_bme_pgi/1207647‑1207938
                                                                                                                                               |                                                                                                                                                    
cgctgtggaggtTTCCGGCGCCAGCCTGGAGACGGGTCGGCACGACCAGTCACTTTCCGCACATCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTT                                                                                                                                                 <  1:71239‑M1/145‑1 (MQ=255)
        attcagcgcattacgttgtagatctgtgttctggtcatttgttcatacgcgtacagagccaataatcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGaa                                                                                                                                         >  2:322708‑M1/136‑149 (MQ=255)
        aGGCGTCAGGTGCCAGCCTGGGGAAGGGGCGGCACGACCAGTTACTTACCGAACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAA                                                                                                                                         <  2:498354‑M1/149‑1 (MQ=255)
         ggCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAG                                                                                                                                        >  2:336727‑M1/1‑149 (MQ=255)
                    ccAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTc                                                                                                                             >  1:425450/1‑149 (MQ=255)
                     cAGCCTGTAGAAGGGGCGGCATGACCAGTTACTTACGGCACTTCGCGCGCTGTTGTTAAGCTGCAAGAATCCGTTTGACGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                                                                                            <  2:386303‑M1/149‑12 (MQ=255)
                     cAGCCTGGGGAAGGGTCGGCACGACCAGTTACTTACCGTACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATTAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                                                                                            <  2:386121‑M1/149‑12 (MQ=255)
                      gttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                           <  1:302997‑M1/28‑1 (MQ=255)
                         gtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                                                                                        <  2:391870‑M1/31‑1 (MQ=255)
                                   ggcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                                                                                              >  1:420698‑M1/1‑124 (MQ=255)
                                       gcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                                                                                          <  2:420698‑M1/149‑30 (MQ=255)
                                               aGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                                                                                                  >  1:166010‑M1/1‑112 (MQ=255)
                                                   acttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCACATCAGCTGGCACGGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag                                                                                              >  1:365912‑M1/1‑108 (MQ=255)
                                                   acttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCCTCGTTTGGGACAAGGctctgcattacctgttgttgtcccggttccgtggccttccg                                                                                              >  2:350303‑M1/1‑108 (MQ=255)
                                                         cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactc                                                                                        >  1:2685‑M1/1‑102 (MQ=255)
                                                          gtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                                                                                       >  1:494068‑M1/86‑149 (MQ=255)
                                                           tacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCGGCAACGGCTACATTTGGACGCGctggtttggtt                                                                                      <  1:55852‑M1/65‑5 (MQ=255)
                                                            aCTTCGCGCGCTGCTGTGAAGCCGCAAGAGTCCGTTTGGCTATATCAAATCGGATGGGACTGTGCAAACGGCTCTCGAAAACCGTGGTTAGGGCGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                                                                                     <  1:151302/149‑1 (MQ=255)
                                                              ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGTAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca                                                                                   <  2:469076‑M1/149‑53 (MQ=255)
                                                              ttCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca                                                                                   >  2:299314‑M1/1‑97 (MQ=255)
                                                                   gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagg                                                                              >  1:132467‑M1/1‑92 (MQ=255)
                                                                     gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct                                                                            <  2:132467‑M1/149‑60 (MQ=255)
                                                                     gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTATCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct                                                                            >  2:417739‑M1/1‑90 (MQ=255)
                                                                              ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtg                                                                   >  2:302997‑M1/66‑149 (MQ=255)
                                                                               taaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtgg                                                                  >  1:74972‑M1/65‑149 (MQ=255)
                                                                               aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatg                                                                  >  2:83612‑M1/1‑80 (MQ=255)
                                                                                gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgc                                                                 >  1:354245‑M1/1‑79 (MQ=255)
                                                                                 ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctccgcattacctgttggtgtcccggttccgtggccttcagccactccgc                                                                                     <  2:439389‑M1/128‑51 (MQ=255)
                                                                                 ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctccgcattacctgttggtgtcccggttccgtggccttcagccactccgc                                                                                     >  1:439389‑M1/1‑78 (MQ=255)
                                                                                     aaGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatgg                                                            >  1:475079‑M1/1‑74 (MQ=255)
                                                                                      ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                           <  1:451182‑M1/92‑1 (MQ=255)
                                                                                        ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGAGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                         >  1:470632‑M1/56‑149 (MQ=255)
                                                                                          tCCGTTTGGCGATATCAAATCGGATGGCGCTGTGCAAACGGATGTCGAAAACCTTGGTTTGGTAGAAGGtgcggcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatgttgccg                                                       >  2:377754‑M1/1‑69 (MQ=255)
                                                                                           gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGATAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCa                                                      <  1:112784‑M1/97‑1 (MQ=255)
                                                                                           gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCa                                                      <  1:382048‑M1/97‑1 (MQ=255)
                                                                                            cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcaaacaggctggctgcaacgcaaaggtgccgct                                                     >  2:148747‑M1/1‑67 (MQ=255)
                                                                                            cGTTTGGCGAAATCAAACCGAACGGCACTGTGCAAACGGCTCCCGAAAACCTTGGTTTGGGAGAACGctctgcattaccttttggtgtcccggttccgtgtccttcagccactccgccaacaggctggctgaaatgcaatggtgccgct                                                     <  2:206485‑M1/149‑83 (MQ=255)
                                                                                             gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                                    <  2:425450/149‑1 (MQ=255)
                                                                                                agcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg                                                 >  1:501289‑M1/48‑149 (MQ=255)
                                                                                                    gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgc                                             >  1:47265‑M1/1‑59 (MQ=255)
                                                                                                               ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacc                                  <  1:372409‑M1/149‑102 (MQ=255)
                                                                                                                gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaataccc                                 >  1:266009‑M1/1‑47 (MQ=255)
                                                                                                                    gaactGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaa                             <  2:447510‑M1/147‑107 (MQ=255)
                                                                                                                      aCTGTGCAAACGGCTCTCGAAACCATTGGTTTGGGAGACGGctctgcgttacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcattggtgccgttttttttggtggagcatgcccggcact                           >  2:224228‑M1/1‑41 (MQ=255)
                                                                                                                       cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactg                          <  2:166010‑M1/149‑110 (MQ=255)
                                                                                                                         gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGc                        <  1:379573/149‑1 (MQ=255)
                                                                                                                          tgCAAACTGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagcatacccggaactggca                       >  1:506750‑M1/1‑37 (MQ=255)
                                                                                                                          tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcc                       >  2:453728‑M1/1‑37 (MQ=255)
                                                                                                                           tttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                                                           >  1:512764‑M1/21‑96 (MQ=255)
                                                                                                                           tttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                                                                           <  2:512764‑M1/76‑1 (MQ=255)
                                                                                                                               ttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAACAACGGCTACATTAAGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaa                  <  2:508203‑M1/133‑1 (MQ=255)
                                                                                                                                tcaatttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATTATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                 >  2:525912‑M1/16‑149 (MQ=255)
                                                                                                                                  attttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGc               <  2:496407‑M1/136‑1 (MQ=255)
                                                                                                                                   ttttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAgcc              <  2:231502‑M1/137‑2 (MQ=255)
                                                                                                                                   ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttat              >  1:316888‑M1/1‑28 (MQ=255)
                                                                                                                                        gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgac         >  1:426628‑M1/1‑23 (MQ=255)
                                                                                                                                             ccTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaata    >  2:301621‑M1/1‑18 (MQ=255)
                                                                                                                                             ccTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacacgctgcctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaatc    >  1:73134‑M1/1‑18 (MQ=255)
                                                                                                                                               TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAATTAGGACGCACCAGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggg  >  1:439552‑M1/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGG  >  NC_000913_3_bme_pgi/1207647‑1207938

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.