New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	3 (0.050)	48 (0.890)	41/250	0.2	97.3%	coding (305/630 nt)	ycfK	e14 prophage; protein StfP
	?	NC_000913	= 1209602	0 (0.000)					pseudogene (54/537 nt)	stfE	e14 prophage; putative side tail fiber protein fragment

ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAA  >  NC_000913/1209461‑1209738                                                                                                                                              |                                                                                                                                         aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                                                                                         <  1:920407‑M2/143‑2 (MQ=255)    aGGCAAGTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                                                                                                      >  1:1119474‑M2/1‑139 (MQ=255)       cAATTTATTTGTCGGATACGCCTTTGCCTGTTCCGGCTATTCTCCAGCAGAAGAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                                   >  2:321277‑M2/1‑136 (MQ=255)         atttatttGTCGGATAAGCCTTTGCCCGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCGCCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACCGAACCGGGACACCACCAGGTAATGCAGAGccgtctccc                                                                                                                                 <  2:478558‑M2/143‑10 (MQ=255)            tatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccaaaa                                                                                                                              <  2:515639‑M2/143‑13 (MQ=255)                tGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaa                                                                                                                          >  1:915356‑M2/1‑127 (MQ=255)                 gTCGGATAAGCCTTTGCCAGTTCCGGGTATTGTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCGGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                                                                                         >  2:531949‑M2/1‑126 (MQ=255)                  tCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                                                                                        <  1:985715‑M2/143‑19 (MQ=255)                  tCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCGGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCACCAGGTAATGCAGAGccttctccccaaccaagg                                                                                                                        <  2:843779‑M2/143‑19 (MQ=255)                     gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCGGAGccttctcccaaaccaaggttg                                                                                                                     <  1:479277‑M2/143‑22 (MQ=255)                         aGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                 >  1:317452‑M2/1‑118 (MQ=255)                            cTTTGCCAGTTCCGAGTATTCTTCATCAGAAAAAGCTGCATCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTTAAGGCCACGGAACCGGGATACCAACAGGTAATGCAGAGtcttctcccaaaccaaggttttcgagag                                                                                                              >  2:158263‑M2/1‑115 (MQ=255)                                  cggttCCGGGTATTCTTCAGCAGGAAGAGCGGCTCCATTGCACTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgttt                                                                                                        >  2:1066222‑M2/3‑109 (MQ=255)                                  cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACAAACAGGTAATGCAGAGccttctcccaagccaaggttttcgagagccgttt                                                                                                        <  1:152166‑M2/143‑35 (MQ=255)                                         gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                                                                                                 <  1:731989‑M2/143‑42 (MQ=255)                                          ggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccacgggtttcgagagccgtttgcacagtg                                                                                                <  2:661600‑M2/143‑43 (MQ=255)                                                 ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagctccgtttgcaccgggccacgcg                                                                                         <  1:23500‑M2/143‑50 (MQ=255)                                                 ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccg                                                                                         <  1:116055‑M2/143‑50 (MQ=255)                                                     gtagaaAAAGCGGCACCAGTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                                                                     >  1:112949‑M2/3‑90 (MQ=255)                                                      cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagcgccgtttgcccagtgccatccgctttg                                                                                    <  2:176381‑M2/143‑55 (MQ=255)                                                      cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttg                                                                                    <  1:738783‑M2/143‑55 (MQ=255)                                                      cagAAAAAGCGGCACCATTGCACTTCAGCCAGTCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggcttccgagagccgtttgcacagtgccagccgacatg                                                                                    <  2:680476‑M2/143‑55 (MQ=255)                                                       agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                                                                   <  1:549018‑M2/143‑56 (MQ=255)                                                         aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatg                                                                                 <  1:505175‑M2/143‑58 (MQ=255)                                                         aaaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata                                                                                 <  1:792849‑M2/143‑58 (MQ=255)                                                           aaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc                                                                               <  1:535322‑M2/143‑60 (MQ=255)                                                               cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgctttgatatcgcca                                                                           <  1:759041‑M2/143‑64 (MQ=255)                                                                ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgtcatccgatttgatatcgccaa                                                                          >  1:1034535‑M2/1‑79 (MQ=255)                                                                    ccATTGCATTTCAGTCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgg                                                                      >  1:4248‑M2/1‑75 (MQ=255)                                                                         gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggatatcgagagccgtttgcacggtgccatccgatttgatatcgccaaacggattct                                                                 <  2:788673‑M2/143‑74 (MQ=255)                                                                           aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcgtg                                                               <  2:397436‑M2/143‑76 (MQ=255)                                                                            tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaccggattctggc                                                              <  2:246518‑M2/143‑77 (MQ=255)                                                                             ttgcgcaagcctggtGGCGGAGTGGCTGACGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcg                                                             >  2:809729‑M2/5‑66 (MQ=255)                                                                                agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggct                                                          <  1:343223‑M2/143‑81 (MQ=255)                                                                                     gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttctcgagagccgtttgcacagtgccatccgatctgatatcgccaaacggattcttgcggcttaacc                                                     <  2:969383‑M2/143‑86 (MQ=255)                                                                                             gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggctcaacagcagcgcg                                             <  2:195003‑M2/143‑94 (MQ=255)                                                                                                 aGTGGCTGAAGGCCACGGAACCGGGATACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaa                                         >  2:25133‑M2/1‑46 (MQ=255)                                                                                                     gCTGAAGGCCGCGGAACCTGGGAACCAACAGGTAATGCAGTGccttcgcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgc                                     >  2:941019‑M2/1‑42 (MQ=255)                                                                                                     gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttccagagccgtttgcacagtgccatccgatttgatatcgccaaacggcttcttgccggttaacagcagcgcgcgcagggc                                     <  1:352513‑M2/143‑102 (MQ=255)                                                                                                          aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaagcggcttcttgcggcttaacagcagcgcgcgaagtgcggtaa                                >  1:264341‑M2/1‑37 (MQ=255)                                                                                                            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagt                              <  1:618777‑M2/143‑109 (MQ=255)                                                                                                            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggagtcttgcggcttaacagcagcgcgcgaagtgcggtaagt                              <  1:894280‑M2/143‑109 (MQ=255)                                                                                                             ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcctaacggattcttgcggcttaacagcagcgcgcgaagtgcggtacgca                             <  2:176891‑M2/143‑110 (MQ=255)                                                                                                              cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcccaacggattcttgcggcttaaccgccgcgcgcgaagtgcggtaagtaa                            <  1:257918‑M2/143‑111 (MQ=255)                                                                                                               aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttagcagcagagcgcgaagtgcgttaagtgac                           <  1:1063646‑M2/143‑112 (MQ=255)                                                                                                                        cgcaaccaacAGGTAATGCAGATccgtctcccaaaccgaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgc                  >  2:922511‑M2/5‑23 (MQ=255)                                                                                                                          acacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccg                <  1:92458‑M2/143‑123 (MQ=255)                                                                                                                           cacCAGCAGGTAATGCAGAGcctcctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgc               >  2:731989‑M2/1‑20 (MQ=255)                                                                                                                              cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacaccagcgcgcgaagtgcggtaagaaactggtcgtgccgcccc            <  2:178648‑M2/143‑127 (MQ=255)                                                                                                                                aCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtcactggtcgtgccgcacctc          <  2:642223‑M2/143‑129 (MQ=255)                                                                                                                                     tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgcattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcca     <  2:847266‑M2/143‑134 (MQ=255)                                                                                                                                       aTGCAGAGccttctcccaaaccacggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccagg   >  1:76940‑M2/1‑8 (MQ=255)                                                                                                                                       aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccagg   <  2:627256‑M2/143‑136 (MQ=255)                                                                                                                                        ggcAGTGcgttcgccgaatccaaggtttgcgagagccgtgtgcacagtgccatgcgatttgatatcgccagacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggc  >  2:952959‑M2/2‑7 (MQ=255)                                                                                                                                              |                                                                                                                                         ATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAA  >  NC_000913/1209461‑1209738

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.