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.

ATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGC  >  NC_000913/4578449‑4578728                                                                                                                                          |                                                                                                                                               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)                                                                                                                                          gCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGc  <  2:591352/143‑1 (MQ=255)                                                                                                                                          |                                                                                                                                               ATTAACATCATCACTTCGCCAAATACTTTACTGAGATTGGCACGATTGATTTCATCTATAATAAAAATATACTTTTTCTCTGGCTGCTCTTTAGCTTGCTGACAAAAATTGTAAAATATGCCGTCTTTACGTCGGAAGCCGACGCCATTCGGACGATAGCCCTGTATAAAATCCTCATAGCTATAAGATTGATGGAACTGAACCATATTGACGCGTTGCGGAGCCTTTTCTCCTGTCAGCAAGTAAGCCAGACGGCGTGCAACAAAGGTTTTTCCAACGC  >  NC_000913/4578449‑4578728

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.