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.

ATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGC  >  NC_000913/4578448‑4578727                                                                                                                                          |                                                                                                                                               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                                                                                                                         <  1:27688/143‑1 (MQ=255)                   ccAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGcc                                                                                                                         <  2:236815/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:619724/1‑143 (MQ=255)                                                                         tttttCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGAcgcg                                                                   >  1:32271/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)                                                                                                                                          gCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGc  <  1:591352/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                               ATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGC  >  NC_000913/4578448‑4578727

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.