New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4578583	54 (1.250)	4 (0.100)	4/270	5.8	7.2%	coding (756/1380 nt)	mcrB	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrB
	?	NC_000913	= 4578586	52 (1.250)					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/4578443‑4578723                                                                                                                                             |                                                                                                                                             ggTTCGATTAACATCATCACTTCGCCAAATACTTGACTGAGGTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTCACGTCCGaa                                                                                                                                            <  1:132782/142‑1 (MQ=255)       aTTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTCCGTCCGCAGCCGAc                                                                                                                                      >  1:222432/1‑143 (MQ=255)        ttAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACg                                                                                                                                     <  2:709784/143‑1 (MQ=255)           ccagcatcaGTTCGGCACATACTTCAATGAGATTGGTACGATGGATTTAATCTATAATAAAAATATACGTTTTCACTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGCTTCCGACgtaa                                                                                                                                  <  1:42411/139‑4 (MQ=255)              tcatcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGACGCCAtgc                                                                                                                               >  2:39935/1‑141 (MQ=255)               catcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGCCGGAAGCCGACGCCATTCg                                                                                                                              >  1:82642/1‑143 (MQ=255)                atcaCTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCgg                                                                                                                             <  2:251023/143‑1 (MQ=255)                  cgcTTCGCCAACTACTGTACTGAGATTGGCACGATTTATTTCATCTATAATAACAATATACTTTTTCTATGGATGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGAc                                                                                                                           <  1:367738/141‑1 (MQ=255)                     ttCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGCAGCAGACGCCATTCGGACGAt                                                                                                                        >  2:364694/1‑143 (MQ=255)                      tCGCGAAATACTTTACCGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATa                                                                                                                       <  1:52607/143‑1 (MQ=255)                         ctaaaTATTTGTCTTAGATTGGCACGATTGATTTCAGCTAGAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                                                                    <  1:683035/141‑1 (MQ=255)                         ccAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                                                                    <  2:27688/143‑1 (MQ=255)                         ccAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                                                                    <  1:236815/143‑1 (MQ=255)                                cTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGtata                                                                                                             >  2:374656/1‑143 (MQ=255)                                 tttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaa                                                                                                            >  2:191543/1‑143 (MQ=255)                                  ttACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCGCTGGCTGCGCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaa                                                                                                           <  2:337273/143‑1 (MQ=255)                                   tACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAATTATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATaaaa                                                                                                          <  2:269138/143‑1 (MQ=255)                                        agaTTGGGACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTc                                                                                                     <  1:144011/143‑1 (MQ=255)                                        agaTTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTc                                                                                                     >  1:145196/1‑143 (MQ=255)                                              gCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCt                                                                                               >  2:611269/1‑143 (MQ=255)                                                aCGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCtat                                                                                             >  1:546888/1‑143 (MQ=255)                                                      gattTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGAtt                                                                                       >  2:38993/1‑143 (MQ=255)                                                            aTCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGa                                                                                 >  2:458838/1‑143 (MQ=255)                                                             tCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGaa                                                                                >  1:52055/1‑143 (MQ=255)                                                               tataATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACt                                                                              <  2:636991/143‑1 (MQ=255)                                                                 taataaAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAAACCTCATAGCTATAAGATTGATGGAACTGa                                                                            >  2:203339/1‑143 (MQ=255)                                                                 taataaAAATATACTTTTTCTCTAGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGa                                                                            <  1:12027/143‑1 (MQ=255)                                                                  aataaACGTATTCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGaa                                                                           <  1:513660/143‑1 (MQ=255)                                                                        gatatACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCatat                                                                     <  1:733266/142‑1 (MQ=255)                                                                         atatACGTTTTCCCTCGCTCCCCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGGCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTCTATAAAATCCTCATCGCTATAAGATTGATGGAACTGAACCATAtt                                                                    <  1:306669/143‑1 (MQ=255)                                                                          gataCTTTTTCTCTTGCTGCTCTTTTGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTg                                                                   <  1:344239/142‑1 (MQ=255)                                                                            taCTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAc                                                                 <  2:461499/143‑1 (MQ=255)                                                                               tttttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgcg                                                              >  2:619724/1‑143 (MQ=255)                                                                               tttttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgcg                                                              >  2:32271/1‑143 (MQ=255)                                                                                    ctctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACCGCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCg                                                         <  2:553539/143‑1 (MQ=255)                                                                                      ctGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGa                                                       <  1:471690/143‑1 (MQ=255)                                                                                       tGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGCCGCGGAg                                                      >  1:362407/1‑143 (MQ=255)                                                                                       tGGCTGCGCTTTAGCTTGCTGACAAAGATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATGGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAg                                                      <  1:207609/143‑1 (MQ=255)                                                                                          ccgctCTCTAGCGTGCTGACAAAAATTGTATAATATGCCGTCGTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCACAGCTATACGATTGATGGAACTGAACCATATTGACGCGTTGCggtgact                                                   <  2:481646/141‑6 (MQ=255)                                                                                              tctTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTtc                                               >  2:716319/1‑143 (MQ=255)                                                                                                  tAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt                                           <  1:619709/143‑1 (MQ=255)                                                                                                  tAGCCTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGACGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCt                                           <  2:273593/143‑1 (MQ=255)                                                                                                      ttGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTcc                                       >  2:276452/1‑142 (MQ=255)                                                                                                      ttGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCa                                       >  2:361442/1‑143 (MQ=255)                                                                                                         cTGACAAAAATTTTAAAATATTTTTTGTTAACACCAAAAGCCGACGCCATTGAGACGATAGGCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCa                                    <  2:138078/143‑1 (MQ=255)                                                                                                          tGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAACATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCaa                                   <  1:252580/143‑1 (MQ=255)                                                                                                            aCAAACATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGGACTTAACCATAGTGACGCGTTGCGGAGCCGTTTCTGCTGTCAGCAAGt                                 >  1:527246/1‑143 (MQ=255)                                                                                                              aaaaaTTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTaa                               <  2:502012/143‑1 (MQ=255)                                                                                                               aaaaTTGGAAAAGATGCCGGCCGTACGTCGGAAGCCGACGCCCGTCGGACGATAGCCCTTTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAg                              <  1:481974/143‑1 (MQ=255)                                                                                                                    tGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGa                         <  1:95153/143‑1 (MQ=255)                                                                                                                        aaagaggCCGTCTTTGCGTCGGTAGCCGAGGCCTGTCGTATGTTGGTCCTGTATAACATCCTCATAGCTATAAGATTCATGGACCTGAATCATATTGACTCGTTTCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGc                     <  1:398115/137‑1 (MQ=255)                                                                                                                          agactccGTCTTTACGTCGGACGCCGACGCCACTCGGACGATAGCCCTGTATAAAATCCTCAGAGCTATAAGATTGATTGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGt                   <  2:255613/138‑1 (MQ=255)                                                                                                                                  tccctACGTCGGAAGCCGACGCCAGTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACaaa           <  2:130896/139‑1 (MQ=255)                                                                                                                                     ttACGTCGGAAGCCGCCGCCATACGGACGATAGCCCTGTATAAATTCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGt        <  1:260458/143‑1 (MQ=255)                                                                                                                                           cGGAAGCCGACGCCATTCGGACGAGAGCCCTGTATAAATTCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTcc  <  1:691587/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                             TGTTCCATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCC  >  NC_000913/4578443‑4578723

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.