New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4036819	42 (0.810)	4 (0.080)	4/252	NT	9.4%	coding (1870/2145 nt)	PP_3557	methyl‑accepting chemotaxis transducer
	?	NC_002947	= 4036880	37 (0.750)					coding (1931/2145 nt)	PP_3557	methyl‑accepting chemotaxis transducer

TTTGCCGTGGTCGCCGATGAGGTGCGCAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCGAGGAGCAAAGCGCGGTAGCCGAGGAAATCAACAGCA  >  NC_002947/4036684‑4036966                                                                                                                                          |                                                                                                                                                    tttGCCGTGGTCGCCGATGAGGTGCGCAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTgg                                                                                                                                        >  2:243906/1‑149 (MQ=255)   tGCCGTGGTCGCCGATGAGGTGCGCAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGc                                                                                                                                       <  6:162010/148‑1 (MQ=255)     ccGTGGTCGCCGATGAGGTGCGCAACCTGGCCCAGCGCACTCAGGAGGCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTg                                                                                                                                    <  4:41324/149‑1 (MQ=255)      cGTGGTCGCCGATGAGGTGCGCAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTgg                                                                                                                                   <  8:94596/149‑1 (MQ=255)                gATGAGGTGCGCAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTCTGGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCg                                                                                                                          <  1:80761/149‑1 (MQ=255)                      gTGCGCAACCTGGCCCACCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGATGGCGTGGAGCAGGTCGGCCAAGc                                                                                                                   >  4:138453/1‑149 (MQ=255)                       tGCGCAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCATTGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGc                                                                                                                   <  1:243906/149‑1 (MQ=255)                        gcgcAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGAGTGGAGCAGGTCGGCCAAGCCg                                                                                                                 >  5:324144/1‑149 (MQ=255)                          gcAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAATCCGTc                                                                                                               <  6:293260/149‑1 (MQ=255)                             aCCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCC                                                                                                                                                     >  1:285991‑M1/1‑108 (MQ=255)                             aCCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCC                                                                                                                                                     <  2:285991‑M1/108‑1 (MQ=255)                                 ggCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGc                                                                                                                   <  6:153693/138‑1 (MQ=255)                                 ggCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGc                                                                                                                   >  5:153693/1‑138 (MQ=255)                                   cccAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCtggccaatgcgctgcaatgccgtgacggcttggccgacctgctccac                                                                                                      >  4:41025‑M1/1‑102 (MQ=255)                                   cccAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTg                                                                                                      <  6:18379/149‑1 (MQ=255)                                     cAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTg                                                                                                      >  4:238454/1‑147 (MQ=255)                                       gcgcACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAgc                                                                                                  >  2:327244/1‑149 (MQ=255)                                       ccgcACTCAGGAGGCGGCGGAAGACCCCCCTCACGTCATCGAAGCCCTGCAGAACGCCACCCGCGAAGTCGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAgc                                                                                                  <  8:149551/148‑1 (MQ=255)                                            cTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCAtt                                                                                             >  5:103673/1‑149 (MQ=255)                                               aGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTgg                                                                                                                                   >  2:166402/1‑108 (MQ=255)                                               aGGAGTCGGGGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTg                                                                                            <  5:280914/147‑1 (MQ=255)                                                ggAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCtggccaatgcgctgcaatgccgtgacggcttggccgacctgctccacgccaccttgggcc                                                                                         >  3:12448‑M1/1‑89 (MQ=255)                                                   gTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAgg                                                                                      <  1:303624‑M2/149‑3 (MQ=255)                                                      ggTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCtggccaatgcgctgcaatgccgtgacggcttggccgacctgctccacgccaccttgggcctggcgg                                                                                   <  4:12448‑M1/149‑67 (MQ=255)                                                      ggTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAggcgg                                                                                   <  8:190308‑M2/149‑6 (MQ=255)                                                       gTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCtggccaatgcgctgcaatgccgtgacggcttggccgacctgctccacgccaccttgggcctggcggt                                                                                  <  3:41025‑M1/149‑68 (MQ=255)                                                          gaagaaACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGcggtgac                                                                               >  3:96944/1‑149 (MQ=255)                                                            agaaACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGcggtgacgg                                                                             <  1:327244/149‑1 (MQ=255)                                                            agaaACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGcggtgacgg                                                                             <  3:238454/149‑1 (MQ=255)                                                             gaaACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGcggtgacggt                                                                            <  2:234568/149‑1 (MQ=255)                                                             gaaACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGcggtgacggt                                                                            <  5:244549/149‑1 (MQ=255)                                                                  ccGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCa                                                                       >  5:208272/1‑149 (MQ=255)                                                                       aGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGa                                                                   <  2:303624/148‑1 (MQ=255)                                                                       aGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGAc                                                                  <  1:59669/149‑1 (MQ=255)                                                                          tCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATg                                                               >  4:278946/1‑149 (MQ=255)                                                                                aaGCCCTGCAGAACGGCACCCGCTAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGcc                                                                                                                  <  6:114907/92‑1 (MQ=255)                                                                                aaGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGcc                                                                                                                  >  5:114907/1‑92 (MQ=255)                                                                                aaGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCtt                                                         <  4:96944/149‑2 (MQ=255)                                                                                  gCCCTGCAGACCGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCa                                                       <  4:93648/149‑1 (MQ=255)                                                                                        cAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGc                                                 <  3:278946/149‑1 (MQ=255)                                                                                           aaCGGCACCCGCGAAGTGGGC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCctc                                              <  4:317460/149‑1 (MQ=255)                                                                                                 aCCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGAGCGCCTCTgccgc                                        >  8:64631/1‑149 (MQ=255)                                                                                                   ccGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTgccgcc                                       >  1:217287/1‑148 (MQ=255)                                                                                                      cgAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCgagg                                   <  7:30511/149‑1 (MQ=255)                                                                                                           tGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACCGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGACGAGGCGCa                               >  5:23051/1‑148 (MQ=255)                                                                                                                  gcgcCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCTCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCGAGGAGCAAAGcgcg                        <  3:138453/149‑1 (MQ=255)                                                                                                                     ccA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCGAGGAGCAAAGCGCGGTa                     <  2:217287/149‑1 (MQ=255)                                                                                                                       a‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCATCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACATGCAGATCGCCTCTTCCGCCGAGGAGCAAAGCGCGGTAGc                   <  6:324144/149‑1 (MQ=255)                                                                                                                           gACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCGAGGAGCAAAGCGCGGTAGcc                  >  1:20303/1‑147 (MQ=255)                                                                                                                                 aGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCGAGGAGCAAAGCGCGGTAGCCGAGGAAAt          >  2:117732/1‑149 (MQ=255)                                                                                                                                     aCCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAggcggtggctgtggtccatggcgccgaccacttcgcgggtgccgttctgcagggcttcgatgacctgacgggtttcttccac      <  3:41025‑M2/149‑83 (MQ=255)                                                                                                                                      ccGCCtggccaatgcgctgcaatgccgtgacggcttggccgacctgctccacgccaccttgggcctggcggtggctgtggtccatggcgccgaccacttcgcgggtgccgttctgcagggcttcgatgacctgacgggtttcttccacc     <  8:190308‑M1/149‑145 (MQ=255)                                                                                                                                      ccGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAggcggtggctgtggtccatggcgccgaccacttcgcgggtgccgttctgcagggcttcgatgacctgacgggtttcttccacc     <  4:12448‑M2/149‑84 (MQ=255)                                                                                                                                      ccGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCGAGGAGCAAAGCGCGGTAGCCGAGGAAATCaa       >  3:216725/1‑147 (MQ=255)                                                                                                                                       cGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCGAGGAGCAAAGCGCGGTAGCCGAGGAAATCAACAg    <  6:208272/149‑1 (MQ=255)                                                                                                                                         cctggccaatgcgctgcaatgccgtgacggcttggccgacctgctccacgccaccttgggcctggcggtggctgtggtccatggcgccgaccacttcgcgggtgccgttctgcagggcttcgatgacctgacgggtttcttccaccgac  <  1:303624‑M1/149‑148 (MQ=255)                                                                                                                                          |                                                                                                                                                    TTTGCCGTGGTCGCCGATGAGGTGCGCAACCTGGCCCAGCGCACTCAGGAGTCGGTGGAAGAAACCCGTCAGGTCATCGAAGCCCTGCAGAACGGCACCCGCGAAGTGGTC‑GGCGCCA‑TGGACCACAGCCACCGCCAGGCCCAAGGTGGCGTGGAGCAGGTCGGCCAAGCCGTCACGGCATTGCAGCGCATTGGCCAGGCGGTGACGGTGATCACCGACATGAACCTGCAGATCGCCTCTGCCGCCGAGGAGCAAAGCGCGGTAGCCGAGGAAATCAACAGCA  >  NC_002947/4036684‑4036966

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