New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	15 (0.220)	22 (0.360)	16/244	NT	50.4%	coding (305/630 nt)	stfP	e14 prophage; uncharacterized protein
	?	NC_000913	= 1209602	30 (0.500)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCT  >  NC_000913/1209461‑1209749                                                                                                                                              |                                                                                                                                                    aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                                              <  2:127315‑M2/149‑8 (MQ=255) aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTACCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttcgc                                                                                                                                              <  2:1171815‑M2/149‑8 (MQ=255)              tttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTGGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGtt                                                                                                                                 <  2:523405/149‑1 (MQ=255)                 gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                              <  1:377569/149‑1 (MQ=255)                  tgggatAAGCCTGTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTTCATTTCAGCCAGACTGTTGCCGGAGTGGCTGATGGCCACGGAACCGGGCCACCTGCAGGTAATGCAGAGcattctcccaaaccaaggttttag                                                                                                                             <  2:932096‑M2/147‑25 (MQ=255)                      aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaa                                                                                                                         <  1:273693/149‑1 (MQ=255)                          gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                                                                                     <  2:605684‑M2/149‑33 (MQ=255)                           ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCCCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCAAGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                                                                                    <  1:776175‑M2/149‑34 (MQ=255)                           ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACATGTAATGCAGCGccttctcacaacccaaggtttgcgagagccgat                                                                                                                    <  1:219519‑M2/149‑34 (MQ=255)                           ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctactcccaaaccaaggttttcgagagcggtt                                                                                                                    <  1:570373‑M2/149‑34 (MQ=255)                                gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCCGCACCATTGCAAATCAGCCAGCCTGTTGGCGGAGTGGCTGAACGCCAACGACCCGGGCAACCAACAGGAAATCCAGAGCCATCTCCCACCCACACTATTATGAAAATGAAGaaaaa                                                                                                               >  2:792991/1‑147 (MQ=255)                                     ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                                                                                          <  2:1116063/149‑1 (MQ=255)                                               tcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                                <  1:169562‑M2/149‑54 (MQ=255)                                               tcttaAGCCGAAAAAGCGGCACCATTGCATTTCAGCCAGCCGGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGTTAATGCAGAGCCTTCTCCCAACCCAACGTTTATGAAAATGAAGAAAAAACAAGCAAATGGcat                                                                                                <  1:794582/149‑1 (MQ=255)                                                cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGAAGAGCCTTCTCCCAAACCAACGGTTATGAAAATGAAGAAATAACAAGCAAATGGcaac                                                                                               >  2:503627/1‑147 (MQ=255)                                                       agAAAAAGCGGCACCATTGCATTTCAGCCATCCTGTTGGCGGAGTGGCTGAAGGCCAGGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                        <  2:1014362/149‑1 (MQ=255)                                                       agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                        >  2:968316/1‑149 (MQ=255)                                                           aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCACCAGGTAATGCAGAGccttctcccaaaccaagggtttcgagcgccgtttgcacagtgccatccgatttgatatcgccaaa                                                                                    <  2:91246‑M2/149‑66 (MQ=255)                                                                    ccATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                           <  2:445179/149‑1 (MQ=255)                                                                         gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACATGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                      <  2:1014674/149‑1 (MQ=255)                                                                           aTTTCAGCCCGACTTTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACATGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATCACAATCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                    <  2:879353/149‑1 (MQ=255)                                                                           aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCAAGGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                    <  1:642606/149‑1 (MQ=255)                                                                            tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                           <  2:601745/125‑1 (MQ=255)                                                                            tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                           >  1:601745/1‑125 (MQ=255)                                                                                  ccagccTGTTGGCGGAGTGGCTGAAGGCCCCGGAACCGGGACCCCAACAGGTAATGCAGAGCCTTCTCCCAACCCAACGTTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAGTTAACATGCtt                                                             >  2:123007/1‑149 (MQ=255)                                                                                    agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgacatccgatttgatatcgccaaacggattcttgcggcttaacagcagc                                                           >  1:907725‑M2/1‑59 (MQ=255)                                                                                     gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                          <  1:578531/149‑1 (MQ=255)                                                                                       cTGTTGGCGGAGTGGCTGAAGGCCACGGACCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                        >  2:591957/1‑149 (MQ=255)                                                                                       cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                        >  2:596984/1‑149 (MQ=255)                                                                                          ttGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATTCAGAGCCTGCTCACAAACCAACGTTTATGAAAATGAATAAATAAAAAGCAAATTGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                     <  1:776653/149‑1 (MQ=255)                                                                                           tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                    <  2:624797/149‑1 (MQ=255)                                                                                             gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGGGATTTAACATGCTTATTGGCTATGt                                                  >  2:1088102/1‑149 (MQ=255)                                                                                               ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                >  1:969121/1‑149 (MQ=255)                                                                                                     gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                          >  2:153750/1‑149 (MQ=255)                                                                                                      cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagt                                         <  2:582813‑M2/149‑109 (MQ=255)                                                                                                      cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagt                                         <  2:1168624‑M2/149‑109 (MQ=255)                                                                                                      cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                         >  1:136039/1‑149 (MQ=255)                                                                                                             ccACGGATCAGGTACACAAACAGGTAATGCAGAGCCTTCTCCAAAACCAACTTTTATGAAAATGAAGAAATTACAAGCTAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                  <  1:262630/149‑1 (MQ=255)                                                                                                              cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtc                                 <  1:982148‑M2/149‑117 (MQ=255)                                                                                                              cACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATg                                 >  2:851234/1‑149 (MQ=255)                                                                                                                     ccGGGACACCAACAGGTACTGGAGAGcctgctcccaacccaaggttttcgcgagccgtttgcactgtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgc                          >  2:1159200‑M2/1‑26 (MQ=255)                                                                                                                      cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCTAATGTCATCATTCCTGCTTTTACCAGGTGGATTTATCATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                         >  1:786059/1‑149 (MQ=255)                                                                                                                        ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccc                       >  2:169562‑M2/1‑23 (MQ=255)                                                                                                                         gACACCAACAGGTAATGCAGAGccttctcccacaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccct                      <  2:1174191‑M2/149‑128 (MQ=255)                                                                                                                          acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                     <  1:851234/149‑1 (MQ=255)                                                                                                                          acaacaacaGGTAATTAAGAGGCTTCTACCAAACCAAAGTTTATGAAAATGAAGAAATAACAAGCTAATGGCATCATTCCTGCTTTTTCAAGGGGGTTTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCATAAAACAg                     <  1:792991/149‑1 (MQ=255)                                                                                                                             ccAACAGGTAATGCAGAGcctgctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatagcgccaaccggattcttgcggcttaacagcagcgcgcgatgtgcggtaagtaactggtcgtgccgccccttctc                  >  1:1171815‑M2/1‑18 (MQ=255)                                                                                                                              cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcc                 <  2:581776‑M2/149‑133 (MQ=255)                                                                                                                              cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcc                 <  2:38770‑M2/149‑133 (MQ=255)                                                                                                                                    gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctg           <  1:17412‑M2/149‑139 (MQ=255)                                                                                                                                    gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACACGCAAATGGCATCATTCCTGCTTTTACGAGGGGGTTTTAACATGCTTATTGGCTATGTCCGCGTATCAACAAATGACCAGAACACAGATCTACAACg           >  1:116836/1‑149 (MQ=255)                                                                                                                                     tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg          <  2:286083‑M2/149‑140 (MQ=255)                                                                                                                                        tGCAGAGcctgctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcac       >  2:377119‑M2/1‑7 (MQ=255)                                                                                                                                         gCAGAGccttctcccaaacccaggtttacgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcgtaacagccccgcgcgctgttcggtaagtaactggtcgtgccgccccttctccaggctggcacc      <  1:72486‑M2/149‑144 (MQ=255)                                                                                                                                         gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcacc      <  1:332875‑M2/149‑144 (MQ=255)                                                                                                                                             cgccttctcccaacccaaggttttcgagagccgtttgcacagtgccatccgatgtgatatcgccaaacggattcttgcggcgtaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccggat  >  1:1142299‑M2/2‑2 (MQ=255)                                                                                                                                              |                                                                                                                                                    ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCT  >  NC_000913/1209461‑1209749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ 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.