New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	8 (0.130)	39 (0.690)	36/250	0.5	80.3%	coding (305/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	= 1209602	12 (0.210)					pseudogene (54/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

AATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTA  >  NC_000913/1209460‑1209742                                                                                                                                               |                                                                                                                                             gaTCAGGCGGTTTATTTGTCGGATAAGCCGTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                              <  2:230790/142‑1 (MQ=255) aaTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGCCACCAACAGGCAATGCagag                                                                                                                                              >  1:328203/1‑143 (MQ=255) aaTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                              >  2:37853/1‑143 (MQ=255)  aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                                                                                             >  1:177262‑M2/1‑142 (MQ=255)   tCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTCAAGGCCACGGAACCGGGACACCAACCGGTAATGCAGAGcc                                                                                                                                            <  2:427958‑M2/143‑3 (MQ=255)     aGGCAATTTATTGGTCGGATAAGCCTTGGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCGTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                                                                                                          >  2:905315‑M2/1‑139 (MQ=255)       gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGAAGAGGCtgcg                                                                                                                                        >  2:123001/1‑140 (MQ=255)       gCAATTTATTTGTCGGATAAGCCTTGGCCAGTTCCGGGTAGTCTGCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                                                                                                        >  1:763428‑M2/1‑137 (MQ=255)        cAATTCACTTGCCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                                       >  1:284222‑M2/1‑136 (MQ=255)           tttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGACCTGTGATGGCATCGGGTAATGCAGAGccttctccca                                                                                                                                    <  1:210697‑M2/143‑11 (MQ=255)                 tgacGGATAAGCAAGGGCCAGTTCGGGGTAGTGGCCAGCAGACAAAGCGGGACCATTGCATTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCaa                                                                                                                              <  1:314554/140‑1 (MQ=255)                  gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                                             <  2:384526‑M2/143‑18 (MQ=255)                      gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttt                                                                                                                         <  2:656052‑M2/143‑22 (MQ=255)                      gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttt                                                                                                                         >  1:133930‑M2/1‑122 (MQ=255)                         aaGCCTTTGCCAGTTCCGGGTAGTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTCATGCAGAGccttctcccaaaccaaggcgttcg                                                                                                                      <  2:1016156‑M2/143‑25 (MQ=255)                           gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgag                                                                                                                    >  1:503094‑M2/1‑117 (MQ=255)                           gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGccttctcccaaaccaaggttttcgag                                                                                                                    <  1:1125447‑M2/143‑27 (MQ=255)                            ccTTTGCTCGTTCCGGGTATGCTTCAGCAGAAAAAGCGGCACCATTGCATTTCATCCAGCATGTTGGCGGATTGTCTGAAGGCCACGTAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                                                   <  2:75607/143‑1 (MQ=255)                            ccTCTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCTTTTCAGCCGGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgaga                                                                                                                   >  2:72712‑M2/1‑116 (MQ=255)                                   cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGCAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtcttcgagcgccgttc                                                                                                            <  2:147359‑M2/143‑35 (MQ=255)                                   cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                            <  1:764451‑M2/143‑35 (MQ=255)                                      ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaacccaaggttttcgagagccgtttgca                                                                                                         <  1:514124‑M2/143‑38 (MQ=255)                                      ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                                                                                         >  1:589542‑M2/1‑106 (MQ=255)                                        ccGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                                                                                       <  1:194268/143‑1 (MQ=255)                                                  ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccg                                                                                             <  1:1068251‑M2/143‑50 (MQ=255)                                                    cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTCATGCAGAGccttctcccacaccaaggttttcgcgcgccgtttgcacagtgccatccgat                                                                                           <  1:1070507‑M2/143‑52 (MQ=255)                                                        agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACCGGTAATGCAGAGccttctcccaaaccaaggttttagagaacagtttgcacagtgcaatccgatttga                                                                                       <  1:209150‑M2/143‑56 (MQ=255)                                                        agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                                                                       >  1:628220‑M2/1‑88 (MQ=255)                                                        agAAAAAGCGGCACCATTGCATTGCAGCCAGCCTGGTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACCCCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                                                                       >  1:936606‑M2/1‑88 (MQ=255)                                                         gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCGGAGCCTTCGCCCGATACAACGTTCATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                      >  2:1005222/1‑143 (MQ=255)                                                                   cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                            >  1:285706/1‑143 (MQ=255)                                                                    aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatgcgatttgatgtcgccaaacg                                                                           <  2:505429‑M2/143‑68 (MQ=255)                                                                        ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaacaaagtattagagagccggttgcacagtgaacaccgatttgatagtgcccaacgcgtt                                                                       <  1:476900‑M2/143‑72 (MQ=255)                                                                          gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                     <  1:952504/143‑1 (MQ=255)                                                                            aTTTCAGCCAGCCTGTTGGCAGAGTGGCTGAACGCCACGGCACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                   >  2:604403/1‑143 (MQ=255)                                                                             tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                                  >  1:656052‑M2/1‑67 (MQ=255)                                                                              ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcg                                                                 <  1:228283‑M2/143‑78 (MQ=255)                                                                                      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaca                                                         <  2:519364‑M2/143‑86 (MQ=255)                                                                                         tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                                      <  2:70157‑M2/143‑89 (MQ=255)                                                                                                ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagccgcgcgcg                                               <  2:626841‑M2/143‑96 (MQ=255)                                                                                                  ggTGGCGGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCTAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTGTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                             <  2:315230/142‑1 (MQ=255)                                                                                                  aGTGGCTGAAGGCCACGGAACCGTGACACCAACAGGTAATGCAGGGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                             <  2:875704/143‑1 (MQ=255)                                                                                                     ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtg                                          >  1:183003‑M2/1‑43 (MQ=255)                                                                                                       cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                        >  1:384526‑M2/1‑41 (MQ=255)                                                                                                        tGAAGGCAACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaacaaggttttcgagagccgtgtgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgggcgaagtgcgg                                       <  1:964556‑M2/143‑104 (MQ=255)                                                                                                        cgaAGGCAACGGAAGCGGGACACCAACAGGTAATGCAGAGtcttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcagattagatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcgg                                       >  1:793362‑M2/2‑40 (MQ=255)                                                                                                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaccgcagcgcgcgaagtgcggtacgt                                  <  2:1037710‑M2/143‑109 (MQ=255)                                                                                                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagga                                 <  1:131135‑M2/143‑110 (MQ=255)                                                                                                                       cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTCTGTACGCGTATCAACAAACGAc                        >  2:286784/1‑143 (MQ=255)                                                                                                                         ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgc                      >  1:37853‑M2/1‑23 (MQ=255)                                                                                                                         ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgc                      <  1:704301‑M2/143‑121 (MQ=255)                                                                                                                              ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgctgtgatatagccaacccgattcttgcggcttcacagcgccgcgcgaagtgcgatacgtaactggtcgtgcctgcc                 <  2:1003786‑M2/143‑126 (MQ=255)                                                                                                                              ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc                 <  1:686268‑M2/143‑126 (MQ=255)                                                                                                                                aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtaccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgccgtaactacccactcgtgccgcccgc               <  2:318523‑M2/143‑128 (MQ=255)                                                                                                                                 aCAGGTAATGCAGCGcctcctcccaacccaaggttttcgagagccgtgtgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccctt              >  2:50429‑M2/1‑15 (MQ=255)                                                                                                                                  cAGGTAATGCAGAGCCTGCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGa             <  1:38333/143‑1 (MQ=255)                                                                                                                                     gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgagttgatatcgccaaacggattcttgcggcttcacagcagcgcgcgaagtgcggtaagttactgttcctgccgccccctcccc          <  1:42157‑M2/143‑133 (MQ=255)                                                                                                                                        acgccGAGcctgctccctatccatggttttcgagggccgtttgcccggtgccagccgatttggtatcgccaaccggattcttgcggcttaacagcagcgcgcgacgtgcggtaagtaactggtcgtgccgccccttctccagg       >  2:1087331‑M2/3‑8 (MQ=255)                                                                                                                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACa       >  1:148652/1‑143 (MQ=255)                                                                                                                                         tGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggc      <  1:502262‑M2/143‑137 (MQ=255)                                                                                                                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc     >  1:8529/1‑143 (MQ=255)                                                                                                                                             gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggc  <  2:844314‑M2/143‑141 (MQ=255)                                                                                                                                               |                                                                                                                                             AATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTA  >  NC_000913/1209460‑1209742

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