New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4578582	54 (1.300)	4 (0.100)	4/270	5.6	7.2%	coding (756/1380 nt)	mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
	?	NC_000913	= 4578585	52 (1.310)					coding (753/1380 nt)	mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCC  >  NC_000913/4578442‑4578722                                                                                                                                             |                                                                                                                                             ggTTCGATTAACATCATCACTTCGCCAAATACTTGACTGAGGTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTCACGTCCGaa                                                                                                                                            <  2:132782/142‑1 (MQ=255)       aTTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTCCGTCCGCAGCCGAc                                                                                                                                      >  2:222432/1‑143 (MQ=255)        ttAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACg                                                                                                                                     <  1:709784/143‑1 (MQ=255)           ccagcatcaGTTCGGCACATACTTCAATGAGATTGGTACGATGGATTTAATCTATAATAAAAATATACGTTTTCACTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGCTTCCGACgtaa                                                                                                                                  <  2:42411/139‑4 (MQ=255)              tcatcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGACGCCAtgc                                                                                                                               >  1:39935/1‑141 (MQ=255)               catcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGCCGGAAGCCGACGCCATTCg                                                                                                                              >  2:82642/1‑143 (MQ=255)                atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCgg                                                                                                                             <  1:251023/143‑1 (MQ=255)                  cgcTTCGCCAACTACTGTACTGAGATTGGCACGATTTATTTCATCTATAATAACAATATACTTTTTCTATGGATGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGAc                                                                                                                           <  2:367738/141‑1 (MQ=255)                     ttCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGCAGCAGACGCCATTCGGACGAt                                                                                                                        >  1:364694/1‑143 (MQ=255)                      tCGCGAAATACTTTACCGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATa                                                                                                                       <  2:52607/143‑1 (MQ=255)                         ctaaaTATTTGTCTTAGATTGGCACGATTGATTTCAGCTAGAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                                                                    <  2:683035/141‑1 (MQ=255)                         ccAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                                                                    <  2:236815/143‑1 (MQ=255)                         ccAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                                                                    <  1:27688/143‑1 (MQ=255)                                cTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGtata                                                                                                             >  1:374656/1‑143 (MQ=255)                                 tttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaa                                                                                                            >  1:191543/1‑143 (MQ=255)                                  ttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCGCTGGCTGCGCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaa                                                                                                           <  1:337273/143‑1 (MQ=255)                                   tACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAATTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaaa                                                                                                          <  1:269138/143‑1 (MQ=255)                                        agaTTGGGACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTc                                                                                                     <  2:144011/143‑1 (MQ=255)                                        agaTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTc                                                                                                     >  2:145196/1‑143 (MQ=255)                                              gCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCt                                                                                               >  1:611269/1‑143 (MQ=255)                                                aCGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCtat                                                                                             >  2:546888/1‑143 (MQ=255)                                                      gattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGAtt                                                                                       >  1:38993/1‑143 (MQ=255)                                                            aTCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGa                                                                                 >  1:458838/1‑143 (MQ=255)                                                             tCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGaa                                                                                >  2:52055/1‑143 (MQ=255)                                                               tataATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACt                                                                              <  1:636991/143‑1 (MQ=255)                                                                 taataaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAAACCTCATAGCTATAAGATTGATGGAACTGa                                                                            >  1:203339/1‑143 (MQ=255)                                                                 taataaAAATATACTTTTTCTCTAGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGa                                                                            <  2:12027/143‑1 (MQ=255)                                                                  aataaACGTATTCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGaa                                                                           <  2:513660/143‑1 (MQ=255)                                                                        gatatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCatat                                                                     <  2:733266/142‑1 (MQ=255)                                                                         atatACGTTTTCCCTCGCTCCCCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTCTATAAAATCCTCATCGCTATAAGATTGATGGAACTGAACCATAtt                                                                    <  2:306669/143‑1 (MQ=255)                                                                          gataCTTTTTCTCTTGCTGCTCTTTTGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTg                                                                   <  2:344239/142‑1 (MQ=255)                                                                            taCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAc                                                                 <  1:461499/143‑1 (MQ=255)                                                                               tttttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgcg                                                              >  1:32271/1‑143 (MQ=255)                                                                               tttttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgcg                                                              >  1:619724/1‑143 (MQ=255)                                                                                    ctctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACCGCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                                                         <  1:553539/143‑1 (MQ=255)                                                                                      ctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGa                                                       <  2:471690/143‑1 (MQ=255)                                                                                       tGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGCCGCGGAg                                                      >  2:362407/1‑143 (MQ=255)                                                                                       tGGCTGCGCTTTAGCTTGCTGACAAAGATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATGGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAg                                                      <  2:207609/143‑1 (MQ=255)                                                                                          ccgctCTCTAGCGTGCTGACAAAAATTGTATAATATGCCGTCGTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCACAGCTATACGATTGATGGAACTGAACCATATTGACGCGTTGCggtgact                                                   <  1:481646/141‑6 (MQ=255)                                                                                              tctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTtc                                               >  1:716319/1‑143 (MQ=255)                                                                                                  tAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt                                           <  2:619709/143‑1 (MQ=255)                                                                                                  tAGCCTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGACGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt                                           <  1:273593/143‑1 (MQ=255)                                                                                                      ttGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTcc                                       >  1:276452/1‑142 (MQ=255)                                                                                                      ttGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCa                                       >  1:361442/1‑143 (MQ=255)                                                                                                         cTGACAAAAATTTTAAAATATTTTTTGTTAACACCAAAAGCCGACGCCATTGAGACGATAGGCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCa                                    <  1:138078/143‑1 (MQ=255)                                                                                                          tGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAACATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCaa                                   <  2:252580/143‑1 (MQ=255)                                                                                                            aCAAACATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGGACTTAACCATAGTGACGCGTTGCGGAGCCGTTTCTGCTGTCAGCAAGt                                 >  2:527246/1‑143 (MQ=255)                                                                                                              aaaaaTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTaa                               <  1:502012/143‑1 (MQ=255)                                                                                                               aaaaTTGGAAAAGATGCCGGCCGTACGTCGGAAGCCGACGCCCGTCGGACGATAGCCCTTTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAg                              <  2:481974/143‑1 (MQ=255)                                                                                                                    tGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGa                         <  2:95153/143‑1 (MQ=255)                                                                                                                        aaagaggCCGTCTTTGCGTCGGTAGCCGAGGCCTGTCGTATGTTGGTCCTGTATAACATCCTCATAGCTATAAGATTCATGGACCTGAATCATATTGACTCGTTTCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGc                     <  2:398115/137‑1 (MQ=255)                                                                                                                          agactccGTCTTTACGTCGGACGCCGACGCCACTCGGACGATAGCCCTGTATAAAATCCTCAGAGCTATAAGATTGATTGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGt                   <  1:255613/138‑1 (MQ=255)                                                                                                                                  tccctACGTCGGAAGCCGACGCCAGTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACaaa           <  1:130896/139‑1 (MQ=255)                                                                                                                                     ttACGTCGGAAGCCGCCGCCATACGGACGATAGCCCTGTATAAATTCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGt        <  2:260458/143‑1 (MQ=255)                                                                                                                                           cGGAAGCCGACGCCATTCGGACGAGAGCCCTGTATAAATTCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTcc  <  2:691587/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                             TGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCC  >  NC_000913/4578442‑4578722

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