New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 282512	24 (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

TGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCCGACAG  >  NC_000913/282372‑282666                                                                                                                                                   |                                                                                                                                                     tGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                                  >  2:635991‑M2/1‑147 (MQ=255)  gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                                  <  2:330673‑M2/150‑5 (MQ=255)  gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                                  <  2:506036‑M2/150‑5 (MQ=255)  gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                                  <  2:618839‑M2/150‑5 (MQ=255)  gATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcac                                                                                                                                                  <  2:94080‑M2/150‑5 (MQ=255)    tCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcaccgg                                                                                                                                               <  1:416717‑M2/151‑8 (MQ=255)       acgaTCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGgcaccggca                                                                                                                                             >  1:626363‑M2/1‑141 (MQ=255)                                  gCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAATCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAAc                                                                                                                  <  2:302658/150‑1 (MQ=255)                                    cTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaag                                                                                                                <  1:581498‑M2/150‑39 (MQ=255)                                          ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagccag                                                                                                          >  2:632213/1‑150 (MQ=255)                                          ttcCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTcagccag                                                                                                          >  2:492619/1‑150 (MQ=255)                                                        tGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcg                                                                                           <  2:518829/151‑1 (MQ=255)                                                             aGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAggcgggcg                                                                                       <  1:293762/150‑1 (MQ=255)                                                             aGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcgg                                                                                       <  2:61404‑M2/150‑64 (MQ=255)                                                             aGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCCGATCAgcaccggcacgcgacgatcgacatgatcgatagcaaagcgggcaatggctttacgctcttcgg                                                                                       <  2:556316‑M2/150‑64 (MQ=255)                                                                        cTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCAt                                                                             <  1:408829/149‑1 (MQ=255)                                                                            gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                                                                        <  2:322944/150‑1 (MQ=255)                                                                            gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                                                                        <  2:352449/150‑1 (MQ=255)                                                                            gcgcCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCgtgg                                                                        <  2:531902/150‑1 (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:60014/151‑1 (MQ=255)                                                                                          aaaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctact                                                         <  1:520406/151‑1 (MQ=255)                                                                                          aaaGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCctact                                                         <  1:421134/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:25805/151‑1 (MQ=255)                                                                                                 ttGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGaa                                                  <  2:28524/151‑1 (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)                                                                                                                                                aTCGGCACCGGCGGCACCTACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCCGAc    <  2:519461/150‑1 (MQ=255)                                                                                                                                                aTCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCCGAc    <  2:469555/150‑1 (MQ=255)                                                                                                                                                aTCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCCGAc    <  2:431665/150‑1 (MQ=255)                                                                                                                                                 tCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCCGACAg  <  1:506046/151‑1 (MQ=255)                                                                                                                                                 tCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCCGACAg  <  2:377683/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGATCGACGATCTGATCAAAGCAGGCGTTGACGGCCTGTTCTTCCTGGGCAGCGGTGGCGAGTTCTCCCAGCTCGGCGCCGAAGAGCGTAAAGCCATTGCCCGCTTTGCTATCGATCATGTCGATCGTCGCGTGCCGGTGCTGATCGGCACCGGCGGCACCAACGCCCGGGAAACCATCGAACTCAGCCAGCACGCGCAGCAGGCGGGCGCGGACGGCATCGTGGTGATCAACCCCTACTACTGGAAAGTGTCGGAAGCGAACCTGATCCGCTATTTCGAGCAGGTGGCCGACAG  >  NC_000913/282372‑282666

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.