New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 28251224 (0.590)7 (0.180) 5/282 2.9 21.3% coding (235/909 nt) yagE 2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage
?NC_000913 = 282518 29 (0.750)coding (241/909 nt) yagE 2‑keto‑3‑deoxy gluconate (KDG) aldolase; CP4‑6 prophage

TGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGG  >  NC_000913/282372‑282659
                                                                                                                                            |                                                                                                                                                   
tGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                           >  2:635991‑M2/1‑147 (MQ=255)
 gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                           <  2:94080‑M2/150‑5 (MQ=255)
 gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                           <  2:618839‑M2/150‑5 (MQ=255)
 gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                           <  2:330673‑M2/150‑5 (MQ=255)
 gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                           <  2:506036‑M2/150‑5 (MQ=255)
   tCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcaccgg                                                                                                                                        <  1:416717‑M2/151‑8 (MQ=255)
      acgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcaccggca                                                                                                                                      >  1:626363‑M2/1‑141 (MQ=255)
                     cAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCtgatcggcaccggcacgcgacgatcgaca                                                                                                                        <  2:316050‑M1/149‑30 (MQ=255)
                                 gCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAATCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAAc                                                                                                           <  2:302658/150‑1 (MQ=255)
                                  ccTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCtgatcggcaccggcacgcgacgatcgacatgatcgatagcaaag                                                                                                         <  2:422576‑M1/151‑45 (MQ=255)
                                  ccTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCtgatcggcaccggcacgcgacgatcgacatgatcgatagcaaag                                                                                                         <  2:43546‑M1/151‑45 (MQ=255)
                                   cTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaag                                                                                                         <  1:581498‑M2/150‑39 (MQ=255)
                                         ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagccag                                                                                                   >  2:492619/1‑150 (MQ=255)
                                         ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagccag                                                                                                   >  2:632213/1‑150 (MQ=255)
                                                       tGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcg                                                                                    <  2:518829/151‑1 (MQ=255)
                                                            aGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcgggcg                                                                                <  1:293762/150‑1 (MQ=255)
                                                            aGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcgg                                                                                <  2:556316‑M2/150‑64 (MQ=255)
                                                            aGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcgg                                                                                <  2:61404‑M2/150‑64 (MQ=255)
                                                                       cTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCAt                                                                      <  1:408829/149‑1 (MQ=255)
                                                                           gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                                                                 <  2:531902/150‑1 (MQ=255)
                                                                           gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                                                                 <  2:352449/150‑1 (MQ=255)
                                                                           gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                                                                 <  2:322944/150‑1 (MQ=255)
                                                                              ccGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCtgatcggcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaact                                                              <  2:61404‑M1/150‑88 (MQ=255)
                                                                              ccGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCtgatcggcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaact                                                              <  2:556316‑M1/150‑88 (MQ=255)
                                                                               cGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATc                                                            <  2:560847/151‑1 (MQ=255)
                                                                               cGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATc                                                            <  2:378235/151‑1 (MQ=255)
                                                                                     gCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAccc                                                       >  2:516924/1‑150 (MQ=255)
                                                                                      cGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAAccc                                                       <  1:103300/149‑1 (MQ=255)
                                                                                         aaaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctact                                                  <  1:421134/151‑1 (MQ=255)
                                                                                         aaaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctact                                                  <  1:520406/151‑1 (MQ=255)
                                                                                         aaaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctact                                                  <  1:60014/151‑1 (MQ=255)
                                                                                          aaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctacta                                                 <  1:14505/151‑1 (MQ=255)
                                                                                             ccATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCtactac                                                <  1:614560/149‑1 (MQ=255)
                                                                                                ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaa                                           <  1:174740/151‑1 (MQ=255)
                                                                                                ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaa                                           <  2:28524/151‑1 (MQ=255)
                                                                                                ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaa                                           <  2:25805/151‑1 (MQ=255)
                                                                                                       cTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCtgatcggcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacag                                     <  1:581498‑M1/150‑113 (MQ=255)
                                                                                                       cTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCg                                    >  2:551266/1‑151 (MQ=255)
                                                                                                       cTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacagg                                    <  2:422576‑M2/151‑108 (MQ=255)
                                                                                                       cTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacagg                                    <  2:43546‑M2/151‑108 (MQ=255)
                                                                                                           gCTATCGCTCATGTCGATCGTCGAGTGCCGGAGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCTCAGCCGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAg                                >  2:126487/1‑151 (MQ=255)
                                                                                                                      tGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacaggccgtcaacgcctg                       <  2:316050‑M2/149‑121 (MQ=255)
                                                                                                                               tcgCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTAttt             <  1:324738/150‑1 (MQ=255)
                                                                                                                                    tgccggtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacaggccgtcaacgcctgctttgatcagatcgt        >  1:626363‑M1/1‑9 (MQ=255)
                                                                                                                                      ccggtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacaggccgtcaacgcctgctttgatcagatcgtcga     <  1:416717‑M1/151‑145 (MQ=255)
                                                                                                                                         gtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacaggccgtcaacgcctgctttgatcagatcgtcgatca  >  2:635991‑M1/1‑4 (MQ=255)
                                                                                                                                         gtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacaggccgtcaacgcctgctttgatcagatcgtcgatc   <  2:330673‑M1/150‑147 (MQ=255)
                                                                                                                                         gtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacaggccgtcaacgcctgctttgatcagatcgtcgatc   <  2:618839‑M1/150‑147 (MQ=255)
                                                                                                                                         gtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacaggccgtcaacgcctgctttgatcagatcgtcgatc   <  2:506036‑M1/150‑147 (MQ=255)
                                                                                                                                         gtgccgatcagcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcggcgccgagctgggagaactcgccaccgctgcccaggaagaacaggccgtcaacgcctgctttgatcagatcgtcgatc   <  2:94080‑M1/150‑147 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGG  >  NC_000913/282372‑282659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ 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.