New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1207805	6 (0.090)	49 (0.800)	41/252	0.4	87.2%	coding (305/630 nt)	stfP	e14 prophage, uncharacterized protein
	?	NC_000913	= 1209602	9 (0.150)					pseudogene (18/501 nt)	stfE	pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related

AATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGT  >  NC_000913/1209460‑1209741                                                                                                                                               |                                                                                                                                            aaTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                             <  2:1092713/143‑1 (MQ=255)       gCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                                                                                                       >  2:958652‑M2/1‑137 (MQ=255)        caaggtagttGTGGGATAAGCCTGTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCAGTTCAGCCAGCCTGTGGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTtctc                                                                                                                                      <  1:217477/135‑1 (MQ=255)             tatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                                                                                                 >  1:640509‑M2/1‑131 (MQ=255)               tttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGCAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                                                                                               <  1:908161‑M2/143‑15 (MQ=255)                   tCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                                                                                           <  1:692632‑M2/143‑19 (MQ=255)                       aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCATGGAACCGGGACACCAACAGGTAATGCAGCGccttctcccaaaccaagcgttt                                                                                                                       <  2:452431‑M2/143‑23 (MQ=255)                       aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccagaccaaggtttt                                                                                                                       >  1:109014‑M2/1‑121 (MQ=255)                       aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtttt                                                                                                                       >  2:102249‑M2/1‑121 (MQ=255)                        tATGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAATAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                                      >  1:1091734‑M2/1‑120 (MQ=255)                        tAAGCCTCCGCCAGTACCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                                      >  2:103153‑M2/1‑120 (MQ=255)                            ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAAGCGGGACACCAACAGGTAAGGCAGAGccttcgcccaaaccaagggtgtcgaga                                                                                                                  <  1:887738‑M2/143‑28 (MQ=255)                              tttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTGAGCACGCTTGTTTGCGGAGTGGAAAAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggctttcgagagc                                                                                                                <  2:1198106‑M2/143‑30 (MQ=255)                              tttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGACACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagcgc                                                                                                                <  2:1195159‑M2/143‑30 (MQ=255)                                 gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGACGGCCACGGAACCGGGACACCACCAGGTACTGCAGAGccttctcccacaccaaggctcccgagagccgc                                                                                                             <  1:914213‑M2/143‑33 (MQ=255)                                   cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATTCAGAGCCTTCTCCCAACCCAACGTTTATGCAAACgaaga                                                                                                           >  2:1082748/1‑143 (MQ=255)                                   cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccacggaattcgaggtctgttt                                                                                                           <  2:1064142‑M2/143‑35 (MQ=255)                                      ttGCGGGTATTCTTCAGCAGAAATAGCGGCACCATTGCATTTGAGCCAGCCTGTTGGCGGAGTGGCTGAGGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                                                                                        >  1:38262‑M2/1‑106 (MQ=255)                                      ttCTGGGTGTTCTTCTGCAGAAAAAGCGGCTCCATTGCTTTTTAGTCAGCCTGATGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaagccaaggttttcgagagccgtttgca                                                                                                        >  2:104591‑M2/1‑106 (MQ=255)                                      ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                                                                                        <  2:1142977‑M2/143‑38 (MQ=255)                                        ctgggtaTTCTTCAGGAGAAAGAGCGGCACCGTTGCATTTCAGCCAGCCTGTTGGCGGAGTCGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcaca                                                                                                      >  1:998778‑M2/3‑104 (MQ=255)                                        ccGGGTATTCTTCAGCAGAAAAGGCGGCACCATGGCATTTCAGCCAGCATGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                                                                                      <  1:786028/143‑1 (MQ=255)                                          gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccgatgggttcgagagccgtttgcccagt                                                                                                    <  1:779848‑M2/143‑42 (MQ=255)                                           ggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTGGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtattcgagagccggttgcacagtg                                                                                                   <  2:333264‑M2/143‑43 (MQ=255)                                              aTTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgcttgcacattgcca                                                                                                <  1:990166‑M2/143‑46 (MQ=255)                                                   tcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                                                           >  1:440766/1‑143 (MQ=255)                                                    ccgtagaaaaAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaacccaaggttttcgagagccgtttgcacagtgccatccgat                                                                                          >  1:797247‑M2/5‑92 (MQ=255)                                                                       actGCATTTCTGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                       <  1:108267/141‑1 (MQ=255)                                                                         ggcCTTCCAGCCAGCGTGTTGGCGGCGTGCCGGAAGGCCCCGAGTCCATTACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattc                                                                     >  1:133260‑M2/2‑71 (MQ=255)                                                                             tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGctttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                                 >  2:129669‑M2/1‑67 (MQ=255)                                                                                    cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGTAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                          <  2:205743/143‑1 (MQ=255)                                                                                     agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaac                                                         <  2:526964‑M2/143‑85 (MQ=255)                                                                                      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaaca                                                        <  1:497417‑M2/143‑86 (MQ=255)                                                                                         gggtggcGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                                     >  1:826864‑M2/4‑55 (MQ=255)                                                                                         cgtCGGCGGAGTAGCCGAAGGCCACGGACCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtattcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca                                                     >  1:276919‑M2/2‑55 (MQ=255)                                                                                           ttGGCTGAGTGGCTTAAGGCCACGGTACCGGGACTCCAACCGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcactgtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagc                                                   >  1:380659‑M2/1‑53 (MQ=255)                                                                                             ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgc                                                 >  2:569984‑M2/1‑51 (MQ=255)                                                                                              gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcg                                                >  2:844766‑M2/1‑50 (MQ=255)                                                                                                ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAACTAACAAGCAAATGGCATCATTCCTGCTTTCACCAGGGGGATTTAACATGATCACTGGCt                                              >  1:591336/1‑143 (MQ=255)                                                                                                  ggtGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaa                                            >  1:322024‑M2/2‑46 (MQ=255)                                                                                                   gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgagatcgcaacacggcgtccagcagcgtaacagcagcgcgcgaag                                           <  2:749471‑M2/143‑99 (MQ=255)                                                                                                      gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttgtcgagagccgtttgcacagtgccatccgatttgatgtcgccaaacggattgttgcggcttagccgccgcgcgcgacagac                                        <  2:153861‑M2/143‑102 (MQ=255)                                                                                                       cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagacccctttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttcacagcagcgcgccaagggcg                                       <  2:688156‑M2/143‑103 (MQ=255)                                                                                                             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagttcgggaagt                                 <  2:617255‑M2/143‑109 (MQ=255)                                                                                                              ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacaccagcgcgcgaagttcggtaagta                                <  1:685443‑M2/143‑110 (MQ=255)                                                                                                               cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggctttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttcacagcagcgcgccaagcgcggtaagcaa                               <  1:1161073‑M2/143‑111 (MQ=255)                                                                                                                  cgcacctGGACTCCAACAGGTAAGGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagcgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactg                            >  1:9208‑M2/4‑30 (MQ=255)                                                                                                                   gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactgc                           <  1:119265‑M2/143‑115 (MQ=255)                                                                                                                     aCCGGGACAGCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCGTCATTCCTGCTTTTACCAGGGGGCTGCAACATGCTTATTGGCTATGTACGCGGATGACCaaacg                         >  1:917576/1‑141 (MQ=255)                                                                                                                         ggACATCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                     <  1:996992/143‑1 (MQ=255)                                                                                                                            cacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaccggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgc                  >  2:197040‑M2/1‑20 (MQ=255)                                                                                                                             acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcc                 <  1:45965‑M2/143‑125 (MQ=255)                                                                                                                             acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcc                 <  1:956577‑M2/143‑125 (MQ=255)                                                                                                                             acCAACAGGTAATGCAGAGccttctcacaaaccaaggttttcgagaaccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagcaactggtcgcgccgca                 <  1:698345‑M2/143‑125 (MQ=255)                                                                                                                              ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccc                >  2:491060‑M2/1‑18 (MQ=255)                                                                                                                               cAACTGTTAATGCATAGccttctcctaaaccaaggttttctagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcccc               >  1:1004669‑M2/1‑17 (MQ=255)                                                                                                                                  cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGACTATGTACGCGAATCAACAAATGACCAGaatcaaca            >  1:126554/1‑137 (MQ=255)                                                                                                                                         tGCAGAGccttctcccaagccgatgttttcaagaacctgtttcgcagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggc     >  1:526964‑M2/1‑7 (MQ=255)                                                                                                                                          gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgcagtgcggtaagtaactggtcgtgccgccccttctccaggct    <  2:1053407‑M2/143‑138 (MQ=255)                                                                                                                                           ccgcgccgtctcccataccaaggttttcgagtgccgtgtgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctg   >  1:689200‑M2/5‑5 (MQ=255)                                                                                                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcccggctg   <  2:640637‑M2/143‑139 (MQ=255)                                                                                                                                           cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccacgccg   <  2:864793‑M2/143‑139 (MQ=255)                                                                                                                                            agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctgg  <  2:938873‑M2/143‑140 (MQ=255)                                                                                                                                               |                                                                                                                                            AATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGT  >  NC_000913/1209460‑1209741

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.